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Perspectives 


Remote Sensing and Human Health: 
New Sensors and New Opportunities 

Louisa R. Beck*t Bradley M. Lobitz.t and Byron L. Woodt 

California State University, Monterey Bay, California, USA; 

NASA Ames Research Center, Moffett Field, California, USA 

Since the launch of Landsat-1 28 years ago, remotely sensed data have been used 
to map features on the earth’s surface. An Increasing number of health studies have 
used remotely sensed data for monitoring, surveillance, or risk mapping, particularly of 
vector-borne diseases. Nearly all studies used data from Landsat, the French Systeme 
Pour rObservatlon de la Terre, and the National Oceanic and Atmospheric 
Administration’s Advanced Very High Resolution Radiometer. New sensor systems are 
In orbit, or soon to be launched, whose data may prove useful for characterizing and 
monitoring the spatial and temporal patterns of Infectious diseases. Increased 
computing power and spatial modeling capabilities of geographic Information systems 
could extend the use of remote sensing beyond the research community Into operational 
disease surveillance and control. This article Illustrates how remotely sensed data have 
been used In health applications and assesses earth-observing satellites that could 
detect and map environmental variables related to the distribution of vector-borne and 
other diseases. 


Remote sensing data enable scientists to 
study the earth's biotic and abiotic components. 
These components and their changes have been 
mapped from space at several temporal and 
spatial scales since 1972. A small number of 
investigators in the health community have 
explored remotely sensed environmental factors 
that might be associated with disease-vector 
habitats and human transmission risk. However, 
most human health studies using remotesensing 
data have focused on data from Landsat's 
Multispectral Scanner (MSS) and Thematic 
Mapper (TM), the National Oceanic and 
Atmospheric Administration (NOAA)'s Advanced 
Very High Resolution Radiometer (AVHRR), and 
France's Systeme Pour I'Observation de la Terre 
(SPOT). In many of these studies (Table 1), 
remotely sensed data were used to derive three 
variables: vegetation cover, landscape structure, 
and water bodies. 

I nternational space agenciesareplanningan 
estimated 80 earth-observing missions in the 
next 15 years (29). During these missions >200 


Address for correspondence: Louisa R. Beck, Earth Systems 
Scienceand Policy, CaliforniaStateUniversity, Monterey Bay, MS 
242-4, NASA Ames Research Center, Moffett Field, CA 94035- 
1000, USA; fax: 650-604-4680; email: Irbeckc^aia.arcnasa.gcv. 


instruments will measure additional environ- 
mental features such as ocean color and other 
currently detectable variables, but at much 
higher spatial and spectral resolutions. The 
commercial sector is also planning to launch 
several systems in the next 5 years that could 
provide complementary data (30). These new 
capabilities will improve spectral, spatial, and 
temporal resolution, allowing exploration of risk 
factors previously beyond the capabilities of 
remote sensing. In addition, advances in 
pathogen, vector, and reservoir and host ecology 
have allowed assessment of a greater range of 
environmental factors that promote disease 
transmission, vector production, and the emer- 
genceand maintenance of disease foci, as well as 
risk for human-vector contact. Advances in 
computer processing and in geographic informa- 
tion system and global positioning system 
technologies facilitate integration of remotely 
sensed environmental parameters with health 
data so that models for disease surveillance and 
control can be developed. 

I n 1998, the N ational Aeronautics and Space 
Administration's (NASA) Center for Health 
Applications of Aerospace Related Technologies 
(CHAART)^ evaluated current and planned 


^CHAART was established at Ames Research Center by NASA's LifeSciences Division, within the Office of Life & Microgravity 
Sciences & Applications, to make remote sensing, geographic information systems, global positioning systems, and computer 
modeling availableto investigators in the human health community. 
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Table 1. Research using remote sensing data to map disease vectors^ 


Disease 

Vector 

Location 

Sensor 

Ref. 

Dracunculiasis 

Cycl ops spp. 

Benin 

TM 

1 


Cycl ops spp. 

Nigeria 

TM 

2 

Eastern equine encephalomyelitis 

Culiseta melanura 

Florida, USA 

TM 

3 

Filariasis 

Culex pi piens 

Egypt 

AVHRR 

4 


Cx. pi piens 

Egypt 

TM 

5,6 

Leishmaniasis 

Phlebotomus papatasi 

SW Asia 

AVHRR 

7 

Lyme disease 

Ixodesscapu laris 

New York, USA 

TM 

8,9 


1 . scapularis 

Wisconsin, USA 

TM 

10 

Malaria 

Anophel es al bi manus 

Mexico 

TM 

11 


An. aibi manus 

Belize 

SPOT 

12 


An. aibi manus 

Belize 

SPOT 

13 


An. aibi manus 

Mexico 

TM 

14 


An. spp. 

Gambia 

AVHRR, Metosat 

15,16 


An. aibi manus 

Mexico 

TM 

17,18 

Rift Valley fever 

Aedes & Cx. spp. 

Kenya 

AVHRR 

19,20 


Cx. spp. 

Kenya 

TM, SAR 

21 


Cx. spp. 

Senegal 

SPOT, AVHRR 

22 

Schistosomiasis 

Biomphalaria spp. 

Egypt 

AVHRR 

23 

Trypanosomiasis 

Glossina spp. 

Kenya, Uganda 

AVHRR 

24 


Glossina spp. 

Kenya 

TM 

25 


Glossina spp. 

West Africa 

AVHRR 

26 


Glossina spp. 

Africa 

AVHRR 

27 


Glossina spp. 

Southern Africa 

AVHRR 

28 


®See Appendix A for explanation of sensor acronyms 


satel I ite sensor systems as a fi rst step i n enabi i ng 
human health scientists to determine data 
relevant for the epidemiologic, entomologic, and 
ecologic aspects of their research, as well as 
developing remote sensing-based models of 
transmission risk. This article discusses the 
results of the evaluation and presents two 
examples of how remotely sensed data have been 
used i n health-rel ated studi es. Thefi rst example, 
a terrestrial application, illustrates how a single 
Landsat TM i mage was used to characterize the 
spatial patterns of key components of the Lyme 
disease transmission cycle in New York. The 
second example, which focuses on the coastal 
environment, shows how remote sensing data 
from different satel I ite systems can be combi ned 
tocharacterizeand map environmental variables 
in the Bay of Bengal that are associated with the 
temporal patterns ofcholera cases in Bangladesh. 
These examples demonstrate how remote sensi ng 
data acquired at various scales and spectral 
resolutions can be used to study infectious 
disease patterns. 

Lyme Disease in the Northeastern 
United States 

During the past 10 years, NASA's Ames 
Research Center has been collaborating with the 


New York Medical Col lege and the Yale School of 
Medicine to develop remote sensing-based 
models for mapping Lyme disease transmission 
risk in the northeastern United States (31,32). 
The first study compared Landsat TM data with 
canine seroprevalence rate (CSR) data summa- 
rized at the municipality level (31). The canine 
data were used as a measure of human exposure 
risk, the assumption being that dogs were more 
likely to acquire tick bites on or near their 
owner's property. The second study usedTM data 
to map relative tick abundance on residential 
properties by using TM-derived indices of 
vegetation greenness and wetness (32). Figure 1 
shows a subset of the TM data used in both 
studies, as well as some of the products (e.g., 
maps) derived from the data. Each product 
illustrated Lymediseasetransmission variables, 
such as vector and reservoir habitats, as well as 
human risk for disease. Figure la shows raw 
Landsat-5 TM data, which are recorded in six 
spectral bands (excluding a seventh thermal 
band) at a spatial resolution of 30 m. These data 
were processed to derive the products shown in 
Figures Ib-d. 

The image in Figure lb was used to explore 
the relationship between forest patch size and 
deer distribution. Because white-tailed deer 
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Figure 1. Landsat Thematic Mapper (TM) satellite data for a 6x6-km area in Westchester County, New York. 
Shown are the raw data (a), as well as products (e.g., maps) derived from the data (b-d) that might be used for 
modeling Lyme disease transmission risk, a) Raw Landsat TM image composed of bands 5, 4, and 3 (mid-infrared, 
near-infrared, and red bands). Vegetation is shown in shades of green, with bare soil and urban areas shown in 
shades of pink and purple. The spatial resolution of these data is 30x30 m. b) Map showing contiguous forest 
patches, derived from a Landsat TM classification. Colors represent discrete patches, with white indicating the 
absence of contiguous forest, c) A 12-class I and cover map derived from the Landsat TM data, d) Composite image 
of three spectral i ndi ces deri ved from the Landsat T M data, showi ng the contri buti ons of scene bri ghtness i n red, 
greenness in green, and wetness in biue. 
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serve as a major host of the adult tick as well as 
its primary mode of transportation, deer 
distribution was a potentially important factor in 
a Lyme disease risk model. 

Figure Ic shows 12 classes used in two 
separate analyses of risk at two different scales 
(31). These classes include water, evergreen 
trees/vegetation, sparse deciduous trees, dense 
deciduous trees, clearings, golf courses (managed 
grass), urban/commercial, miscellaneous urban, 
residential-lawn, residential-sparse vegetation, 
residential -medium vegetation, and residential 
high -vegetation. I n the first scale, the amount of 
remotely sensed deciduous forest was positively 
correlated (r-0.82) with canine exposure to 
Borrelia burgdorferi, as indicated by CSR data 
summarized by municipality. In the second 
analysis, a linear regression of the residential- 
high vegetation pixels (i.e., wood-edge) and CSR 
data resulted in a correlation coefficient of 0.84— 
indicating that human-host contact risk (e.g., 
deer leaving the forest to feed on residential 
ornamental vegetation) might be a good measure 
of human-vector contact risk. 

The i mage i n F igure Id was derived from the 
LandsatTM data by a Tassel ed Cap Transforma- 
tion (33). Tasseled Cap greenness and wetness 
were positively correlated with tick abundance 
on residential properties in this study area (32). 

Cholera in Bangladesh 

The second example of the use of remotely 
sensed data to provide information for health 
research and applications concerns cholera in 
Bangladesh. I n this study, described by Lobitz 
et al. (34), remotely sensed datasets, down- 
loaded from the I nternet at no cost, were used 
to search for temporal patterns in the Bay of 
Bengal associated with cholera outbreaks in 
Bangladesh. 

Figure 2a shows a color-infrared image of 
the Ganges River, where it empties i nto the Bay 
of Bengal. These data, which were acquired by 
NOAA'sAVFI RR sensor, have a spatial resolution 
of 1.1 km. The sediment load, transported to the 
Bay of Bengal by the Ganges and Brahmaputra 
rivers, includes nutrients that could support 
plankton blooms. Plankton is an important 
marine reservoir of Vibrio cholerae, which 
attaches primarily to zooplankton, which, in 
turn, is associated with phytoplankton (35). 

In Figure 2b, the AVFIRR data shown in 
Figure 2A were processed to show sea surface 


temperature (SST) (36). Because these data are 
for large-area studies, they have been processed 
at a spatial resolution of 18 km. Figure 2c 
represents sea surface height (SSFI) anomaly 
data derived from theTOPEX/Poseidon satellite 
(37). These data have a spatial resolution of 1 
degree. I ncreases in SST and SSFI have preceded 
cholera outbreaks in Bangladesh (34). 

I n the next 15 years, new sensors will provi de 
valuable data for studies of infectious diseases 
similar to the ones described here. For Lyme 
disease, new sensors could provide similar 
information about ecotones, human settlement 
patterns, or forests. These sensors include 
ARIES-1, scheduled for launch by Australia; 
CCD and IR/MSS sensors onboard CBERS, 
launched by China and Brazil in late 1999; 
Ikonos, a commercial satellite with 4-m spatial 
resolution; LI SS 1 1 1 , onboard the orbiting I ndian 
I RS-IC and - ID satellites; and ASTER, onboard 
the recently launched Terra satellite. Informa- 
tion from these sensors could also be used to 
address other vector-borne diseases, such as 
malaria, schistosomiasis, trypanosomiasis, and 
hantavirus, whose patterns are likewise influ- 
enced by environmental variables. 

SeaWiFS, the Sea-Viewing Wide Field-of- 
View Sensor, with its increased spectral 
resol ut i on of 1. 1 km, i s al ready provi di ng i magery 
critical to understanding the temporal and 
spatial pattern of cholera risk (35). This sensor 
was specifically designed to gather information 
about ocean color (38) (Figure 2d). 

Sensor Evaluation Project^ 

CFIAART evaluated data from current and 
planned satellite instruments for mapping, 
surveillance, prediction, and control of human 
disease transmission activities, including vector 
ecology, reservoir and host ecology, and human 
settlement patterns. From hundreds of potential 
sensors, 54 were identified that were current (or 
would be launched within the next 5 years), 
operational (not reserved for the scientific 
community), and digital (not photographic). 

Beginningin 1985, NASA has held a series of 
workshops to elicit input from the health 
community on the use of remote sensing in the 
areas of entomology, ecology, epidemiology, 
vector control, and infectious diseases. In 
addition, NASA has participated in sessions on 
remotesensingand health at professional meetings 
sponsored by national and international health 


^The information gathered during the CHAART sensor evaiuation process is avaiiabie at http://geo.arc.nasa.gov/sge/heaith/ 
sensor/sensor.htmi . 
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Figure 2. Datasets used to model the temporal patterns of cholera outbreaks in Bangladesh, a) Advanced Very 
High Resolution Radiometer (AVH RR) satellite image showing the mouth of the Ganges River and the Bay of 
Bengal. Vegetation is shown in shades of red and water in shades of blue. The spatial resolution of these data is 
1.1 km. b) Sea surface temperature data, derived from AVHRR thermal bands. Temperatures range from low 
(purple) to high (red), c) Sea surface height data, derived from TOPEX/Poseidon satellite data. The spatial 
resolution of these data is 1 degree, d) I mage derived from theSea-Viewing Wide Field-of-View Sensor (SeaWiFS) 
showing chlorophyll concentration, ranging from low (blue) to high (red). These satellite data have a nominal spatial 
resolution of 1.1 km. 


organizations. On the basis of this experience as 
well asa review of the scientificliterature(Table 
1), there does not appear to be consensus in the 
health community regarding requirements for a 
remote sensing system. Some investigators use 
remotely sensed data to resolve questions 
regarding the relationship between an aspect of 
disease transmission and an environmental 


variable. Other researchersalready have a model 
of disease transmission and havespecificspatial, 
temporal, or spectral requirements for the 
remotely sensed variables. 

No single spatial, temporal, or spectral 
resolution is universally appropriate for under- 
standing the transmission risk for any disease, 
given the variety of vectors, reservoirs, hosts, 
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geographic locations, and environmental vari- 
ables associated with that disease. Therefore, in 
evaluating the existing sensors, CHAART used 
an approach that allowed individual investiga- 
tors to i dentify satel I i tedata appropri atefor thei r 
own needs. This approach defined 16 groups of 
physical factors that could be used for both 
research and applications. Each factor is 
essentially an envi ron mental variablethat might 
haveadirect or indirect bearing on the survival of 
pathogens, vectors, reservoirs, and hosts. These 
factors may also affect many types of non-vector- 
borne diseases, such as waterborne diseases. The 
factors are vegetation or crop type, vegetation 
green-up, ecotones, deforestation, forest patches, 
flooded forests, general flooding, permanent 
water, wetlands, soil moisture, canals, human 
settlements, urban features, ocean color, SST, 
and SSH. Precipitation, humidity, and surface 
temperature were not included because deriving 
these measurements from raw data requires 
highly specialized processing and calibration, 
routinely performed by qualified groups who 
often make the information availableon Internet 
websites (Appendix B). 

The sensor evaluation project generated a 
series of tables that associated each of the 16 
factors with the 54 sensors according to spatial, 
temporal, and spectral characteristics. For 
example, factors requiring frequent monitoring, 
such as vegetation green-up, are linked with 
sensors with shorter repeat overpasses. Simi- 
larly, factors requiring very high spatial 
resolution, such as mapping urban features, are 
linked with sensors having a spatial resolution of 


10 m or less, regardless of their temporal or 
spectral resolutions. 

Perhaps the broadest use of Landsat and 
SPOT data has been to identify and map 
vegetation or crop types. This factor is important 
because the distribution of vegetation types 
integrates the combined impact of rainfall, 
temperature, humidity, topographic effects, soil, 
water availability, and human activities. Nearly 
all vector-borne diseases are linked to the 
vegetated environment during some aspect of 
their transmission cycle; in many cases, this 
vegetation can be sensed remotely from space. 
The spatial and temporal distribution of vector or 
reservoir and host species may relate to the 
occurrence and distri bution of specific vegetation 
or crop types, not simply to whether an area has 
forests or grasslands. For example, food and 
cover preferences of the white-tailed deer, the 
host for adult ticks that transmit Lyme disease in 
the northeastern United States, might well 
encourage deer to I i ve near certai n types of forest. 
Crop-type information may also be important for 
studying the effects of pesticides (e.g., vector 
resistance; illnesses caused by exposure to toxins). 

The sensor evaluation procedure has identi- 
fied many potentially useful sensors for mapping 
vegetation and crop type beyond the Landsat and 
SPOT systems (Table 2). A ground resolution 
threshold of 30 m was used asthe upper limit for 
exploringtherelationship between vegetation (or 
crop) type and disease vectors, reservoirs, and 
hosts; above 30 m, vegetation and crop type are 
more difficult to ascertain. Many of the sensors 
could also be used for mapping the boundary 


Table 2. Current and proposed sensor systems for identifying and mapping vegetation and crop type^ 


T emporal 
resolution (d) 

1-5 

Spatial resolution'’’ (m) 
6-10 

11-30 

Daily 

2-7 

(QuickBird)^ 

(Orbview-3,4) 

(Almaz-lb MSU-E2) 

(ALOS AVNIR-2) 


(QuickBird) 

(ALOS AVNIR-2) 

(ARIES-1) 



(Orbview-4) 

(SPOT-5a,b3xHRG) 

SPOT-4 2xH RVIR 

8-14 

1 konos 

Priroda/Mir MOMS-2P 

Priroda/Mir MOMS-2P 

15-30 



Terra ASTER 
IRS-1C,D LISS III 
Landsat TM 
Landsat-7 ETM -l- 
SPOT-2 2xH RV 

>30 


(ALOS AVNIR-2) 

(ALOS AVNIR-2) 


®This matrix is the output from an interactive search with the search engine located at http://geo.arc.nasa.gov/sge/health/ 
sensor/senchar.html . 

‘’See Appendix A for explanations of sensor acronyms. 

'’Sensors in parentheses have not yet been launched. 
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between vegetation types, or ecotones, which 
provide habitat for insects and animalscritical to 
the maintenance and transmission of vector- 
borne diseases. These edges may be areas for 
increased risk for vector-human contact, as 
indicated by the relationship between Lyme 
disease transmission and suburban encroach- 
ment into forested areas in the northeastern 
United States. The movement of humans into 
forested edges where potential vectors are 
established could also be important for predict- 


ing malaria or yellow fever transmission. 

The list of 16 factors used in the CHAART 
evaluation includes some that have not yet been 
quantified because available sensors do not 
provide adequate spatial, spectral, or temporal 
resolutions. Two of these factors are briefly 
described below to illustrate how remotely 
sensed data might be used to explore their 
potential links to human health. More links 
between the factors and various diseases are 
listed in Tables. 


Table 3. Potential links between remotely sensed factors and disease 


Factor 

Disease 

Mapping opportunity 

Vegetation/crop type 

Chagas disease 

Flantavirus 

Leishmaniasis 

Lyme disease 

Malaria 

Plague 

Schistosomiasis 
T rypanosomiasis 
Yellow fever 

Palm forest, dry & degraded woodland habitat for triatomines 
Preferred food sources for host/reservoirs 
Thick forests as vector/reservoir habitat in Americas 
Preferred food sources and habitat for host/reservoirs 

Breeding/resting/feeding habitats; crop pesticides vector resistance 

Prairie dog and other reservoir habitat 
Agricultural association with snails, use of human fertilizer 
Glossina habitat (forests, around villages, depending on species) 
Reservoir (monkey) habitat 

Vegetation green-up 

Flantavirus 
Lyme disease 
Malaria 
Plague 

Rift Valley fever 
T rypanosomiasis 

Timing of food sources for rodent reservoirs 

Flabitat formation and movement of reservoirs, hosts, vectors 

Timing of habitat creation 

Locating prairie dog towns 

Rainfall 

Glossina survival 

Ecotones 

Leishmaniasis 
Lyme disease 

Flabitats in and around cities that support reservoir (e.g., foxes) 

Ecotonal habitat for deer, other hosts/reservoi rs; human/vector contact risk. 

Deforestation 

Chagas disease 
Malaria 

Yellow fever 

New settlements in endemic-disease areas 
Flabitat creation (for vectors requiring sunlit pools) 

Flabitat destruction (for vectors requiring shaded pools) 

Migration of infected workers into forests where vectors exist 
Migration of disease reservoirs (monkeys) in search of new habitat 

Forest patches 

Lyme disease 
Yellow fever 

Flabitat requirements of deer and other hosts, reservoirs 
Reservoir (monkey) habitat, migration routes 

F 1 ooded forests 

Malaria 

Mosquito habitat 

Flooding 

Malaria 

Rift Valley fever 
Schistosomiasis 
St. Louis encephalitis 

Mosquito habitat 

Flooding of dambos, breeding habitat for mosquito vector 
Flabitat creation for snails 
Flabitat creation for mosquitoes 

Permanent water 

Filariasis 

Malaria 

Onchocerciasis 

Schistosomiasis 

Breeding habitat for Mansonia mosquitoes 
Breeding habitat for mosquitoes 
Simulium larval habitat 
Snail habitat 

Wetlands 

Cholera 

Encephalitis 

Malaria 

Schistosomiasis 

Vi brio Choi erae associated with inland water 
Mosquito habitat 
Mosquito habitat 
Snail habitat 

Soil moisture 

Flelminthiases 
Lyme disease 
Malaria 

Schistosomiasis 

Worm habitat 
Tick habitat 
Vector breeding habitat 
Snail habitat 

Canals 

Malaria 

Onchocerciasis 

Schistosomiasis 

Dry season mosquito-breeding habitat; ponding: leaking water 
Simulium larval habitat 
Snail habitat 

Fluman settlements 

Diseases 

Source of infected humans: populations at risk for transmission in general 

Urban features 

Chagas disease 
Dengue fever 
Filariasis 
Leishmaniasis 

Dwellings that provide habitat for triatomines 
Urban mosquito habitats 
Urban mosquito habitats 
Flousing quality 

Ocean color (Red tides) 

Cholera 

Phytoplankton blooms; nutrients, sediments 

Sea surface temp. 

Cholera 

Plankton blooms (cold water upwelling in marine environment) 

Sea surface height 

Cholera 

Inland movement of Vibrio-contaminated tidal water 
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Urban Features 

The detection of urban features requires 
higher spatial resolution systems than needed for 
detecting the presence of human settlements. 
Some disease vectors are associated with specific 
urban features such as housing type, which can 
only be detected by sensors with very high spatial 
resolution. I nthefuture, new sensors may be able 
to provide information on the urban environment 
(Table 4). 


Table 4. Current and proposed sensor systems for 
identifying and mapping urban features^ 


Temporal 

resolution 

(d) 

Spatial resolution‘s (m) 

1-5 6-10 

Daily 

(QuickBird)^ 


2-7 

(ALOS AVNIR-2) 

(Almaz-lb MSU-E2) 


(Orbview-3,4) 

(ALOS AVNIR-2) 


(QuickBird) 

(ARIES-1) 


(SPOT-5a,b3xHRVIR) 

IRS-1C,D PAN 
(Orbview-4) 

SPOT-4 2xHRVIR 
(SPOT-5a,b3xHRVIR) 

8-15 

1 konos 

IRS-1C,D PAN 
Priroda/Mir M0MS-2P 

15-30 


IRS-1C,D PAN 
SPOT-2 2xH RV 

>30 

(ALOS AVNIR-2) 

(ALOS AVNIR-2) 


®This matrix is the output from an interactivesearch with the 
search engine located at http://geo.arc.nasa.gov/sge/health/ 
sensor/senchar.html 

*^See Appendix A for explanations of sensor acronyms. 
■^Sensors in parentheses have not yet been launched. 


Soil Moisture 

Wet soils indicate a suitable habitat for 
species of snails, mosquito larvae, ticks, and 
worms. Several types of sensors can detect soil 
moisture, including synthetic aperture radars 
(SARs), shortwave-infrared, and thermal-infra- 
red sensors (Table 5). SARs are particularly 
i mportant for sensi ng ground conditions i n areas 
of cloud cover or vegetation canopy cover, two 
factors commonly found in the tropics. 

Conclusions 

The extent to which remotely sensed data 
are used for studying the spatial and temporal 
patterns of disease depends on a number of 
obstacles and opportunities. Many of the 
obstacles— including cost, inadequate spatial, 
spectral, or temporal resolutions, and long 
turnaround times for products— have restricted 
the use of remote sensing within the user 
community as a whole. Many of these barriers 
will be addressed by new sensor systems in the 
next 5 years. The recently launched Landsat-7 
ETM-h sensor, for example, is now providing 30- 
m multispectral data, a 15-m panchromatic 
band, and an improved 60-m thermal infrared 
band, all at a cost that is an order of magnitude 
less than current Landsat-5TM data. 

With the higher spatial and spectral 
resolutions, morefrequent coverage, lower price. 


Table 5. Current and proposed sensor systems for identifying and mapping soil moisture^ 


T emporal 





resolution 


Spatial Resolution's (m) 


(d) 

11-30 

101-500 

501-1,000 

1,001-4,000 

Daily 

2-7 

(Almaz-lb SAR-70)^ 


(ADEOS II GLI) 

NOAA AVHRR 


(ARIES-1) 


(Almaz-lb MSU-SK) 



(E N VI SAT-1 ASAR) 


(Almaz-lb SROSM) 



Radarsat 


(E N VI SAT-1 AATSR) 



SPOT-4 2XHRVIR 


Terra MODIS 
(EOS PM-1 MODIS) 
Resurs-01 N2,3 
MSU-SK 


8-14 

(LightSAR) 

Priroda/Mir MOMS-2P 


Priroda/Mir MSU-SK 


15-30 

Terra ASTER 
ERS-1,2 AMI-SAR 
Landsat TM 
Landsat-7 ETM-l- 

Landsat TM TIR 



>30 

ERS-1,2 AMI-SAR 





®This matrix is the output from an interactive search with the search engine located at http://geo.arc.nasa.gov/sge/health/ 
sensor/senchar.html . 

*^See Appendix A for explanations of sensor acronyms. 

■^Sensors in parentheses have not yet been launched. 


Emerging Infectious Diseases 


224 


Vol. 6, No. 3, May-]une 2000 



Perspectives 


and i ncreased avai I abi I ity offered by the range of 
new sensors, human health investigators should 
be able to extract many more environmental 
variables than previously realized. These 
improvements will provide new opportunities to 
extend the uses of remote sensing technology 
beyond a few vector-borne diseases to studies of 
water- and soil-borne diseases (for example, 
cholera and schistosomiasis [waterborne] and 
the helminthiases) and the mapping of human 
settlements at risk. The next generation of earth- 
observing remote sensing systems will also allow 
investigators in the human health community to 
characterize an increasing range of variables key 
to understanding the spatial and temporal 
patterns of disease transmission risk. These 
improved capabilities, when combined with the 
i ncreased computi ng power and spati al model i ng 
capabilities of geographic information systems, 
should extend remote sensing into operational 
disease surveillance and control. 
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Appendix A. Acronyms used in the text and tables 


Acronym 

Mission 

1 nstru ments 

Country 

ADEOS II 

Advanced Earth Observation Satellite 

GLI 

J apan 

ALOS 

Advanced Land Observing Satellite 

AVNIR 

J apan 

ARIES 

Australian Resource 1 nformation & Environment Satellite 

ARIES 

Australia 

CBERS 

China-Brazil Earth Resources Satellite 

CCD, IR/MSS 

China/Brazil 

ENVISAT 

Environmental Satellite 

AATSR, ASAR 

Europe 

EOS 

Earth Observation System 

ASTER, MODIS 

USA 

ERS-2 

ESA (European SpaceAgency) Remote Sensing 

AMI-SAR 

Europe 

IRS 

Indian Remote Sensing Satellite 

PAN, LISS 

1 ndia 

NOAA 

National Oceanographic & Atmospheric Administration 

AVHRR 

USA 

SPOT 

Systeme Pour TObservation de la Terre 

HRV, HRVIR 

F ranee 

Acronym 

1 nstrument 

Mission 

Country 

A AT SR 

Advanced Along Track Scanning Radiometer 

ENVISAT 1 

ESA 

AMI-SAR 

Active M icrowave 1 nstru mentation Synthetic Aperture Radar 

ERS-1, 2 

ESA 

ASAR 

Advanced Synthetic Aperture Radar 

ENVISAT 1 

ESA 

ASTER 

Advanced Spaceborne Thermal Emission & Reflection Radiometer 

Terra 

USA 

AVHRR 

Advanced Very Fligh Resolution Radiometer 

NOAA 

USA 

AVNIR 

Advanced Visible & Near Infrared Radiometer 

ALOS 

J apan 

CCD 

Charged Couple Device Camera 

CBERS 

China/Brazil 

ETM + 

Enhanced Thematic Mapper Plus 

Landsat-7 

USA 

GLI 

Global Land 1 mager 

ADEOS II 

J apan 

HRV 

Fligh Resolution Visible 

SPOT 1, 2 

F ranee 

HRVIR 

Fligh Resolution Visible& Infrared 

SPOT 4, 5 

F ranee 

IR-MSS 

Infrared-Multispectral Scanner 

CBERS 

China/Brazil 

LISS III 

Linear 1 maging Self-Scanning System 

IRS-IC, D 

1 ndia 

MODIS 

Moderate Resolution Imaging Spectro Radiometer 

Terra, EOS PM 1-3 

USA 

MOIVIS-2P 

Modular Optoelectronic Multispectral Scanner 

Priroda/Mir 

Russia 

MSU-E2 

Multizone Eli gh-Resoluti on Electronic Scanner 

Almaz-IB 

Russia 

MSU-SK 

Multizone Middle-Resolution Optomechanical Scanner 

Almaz-IB 

Russia 



Priroda 

Russia 



Resurs-01, 02 

Russia 

PAN 

Panchromatic 

IRS-IC, D 

India 

PAN 

Panchromatic 

1 konos-2 

Space 1 maging 

SAR-70 

Synthetic Aperture Radar (70 cm) 

Almaz-IB 

Russia 

SeaWIFS 

Sea-Viewing Wide Field-of-View Sensor 

TOPEX/Poseidon 

F rance/USA 

SROSM 

Spectroradiometer for Ocean Satellite Monitoring 

Almaz-IB 

Russia 

TM 

Thematic Mapper 

Landsat 

USA 

Acronym 

Miscellaneous 



ESA 

European SpaceAgency 



TIR 

Thermal Infrared 




Appendix B. Partial list of Internet locations that contain (or will contain) products derived from remotely sensed 
precipitation, moisture, relative humidity, and surface temperature data 


Parameter 

Mission/ sensor^ 

Spatial 

Temporal 

Web site 

Precipitation 

TRMM/FMI 

10 km 

D 

daac.gsfc.nasa.gov/CAMPAIGN_DOCS/hydrology/ 

hd_main.html 

Precipitation 

TRMM/FMI 

S" 

M 

daac.gsfc.nasa. gov/C A MPAIGN_DOCS/hydrology/ 
hd_main.html 

Precipitation 

Terra7DAS‘‘ 

1° 

8-day 

eos-am.gsfc.nasa.gov 

Precipitation 

GOES/S ounder 

10 km 

hr 

www.nndc.noaa.gov/phase3/productaccm.htm 

Precipitation rate 

Terra7DAS‘‘ 

2° 

8/day 

eos-am.gsfc.nasa.gov 

Precipitation amount 

TOVS/MSU 

0.5° 

D, M, Y 

ghrc.msfc.nasa.gov/ghrc/list.html 

Moisture 

GOES/I mager 

8 km 

hr 

www.nndc.noaa.gov/phase3/productaccm.htm 

Relative humidity 

Terra7DAS‘‘ 

2° 

4/day 

eos-am.gsfc.nasa.gov 

Surface temperature 

TerraVMODIS 

1 km 

D 

eos-a m. gsfe. n asa . gov 

Surface temperature 

TerraVMODIS 

1° 

D, 8-day, M 

eos-am.gsfc.nasa.gov 

Surface temperature 

Terra7DAS‘‘ 

2° 

8/day 

eos-am.gsfc.nasa.gov 

Surface temperature 

Envisat/AATSR 

17 km 

D 

envi sat.estec.esa.nl/envisat-welcom.html 

Surface temperature 

Envisat/AATSR 

50 km 

D 

envi sat.estec.esa.nl/envisat-welcom.html 

Surface temperature 

ERS/ATSR 

1 km 

W 

earthl.esrin.esa.it/ERS 

Surface temperature 

Meteosat/VISSR 

5 km 

48/day 

www.eumetsat.de/en 

Surface temperature 

GOES/I mager 

4 km 

hr 

www.nndc.noaa.gov/phase3/productaccm.htm 


®See Appendix A for explanation of sensor acronyms, 
‘’hr, hourly: D, daily; W, weekly: monthly: Y, yearly. 

'’F uture launch. 

‘‘Data Assimilation System. 
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A Dynamic Transmission Modei 
for Predicting Trends in 
Heiicobacter pyiori and Associated 
Diseases in the United States 

Marcia FT. Rupnow* Ross D. Shachter* Douglas K. Owens, *t and 

Julie Parsonnet* 

^Stanford University, Stanford, California, USA; tDepartment of Veterans 
Affairs Health Care System, Palo Alto, California, USA 


To assess the benefits of intervention programs against Helicobacter pylori 
infection, we estimated the baseline curves of its incidence and prevalence. We 
developed a mathematical (compartmental) model of the intrinsic dynamics of H. pylori, 
which represents the natural history of infection and disease progression. Our model 
divided the population according to age, infection status, and clinical state. Case- 
patients were followed from birth to death. A proportion of the population acquired H. 
py/or/ infection and became ill with gastritis, duodenal ulcer, chronic atrophic gastritis, or 
gastric cancer. We simulated the change in transmissibility consistent with the incidence 
of gastric cancer and duodenal ulcer over time, as well as current H. pylori prevalence. 
In the United States, transmissibility of H. pylori has decreased to values so low that, 
should this trend continue, the organism will disappear from the population without 
targeted intervention; this process, however, will take more than a century. 


Helicobacter pylori, a common human 
bacterial pathogen (1), causes peptic ulcer 
disease, gastric cancer, and gastric mucosa- 
associated lymphoid tissue lymphoma. Current 
U .5. guidelines from the National I nstitutes of 
Health recommend antimicrobial treatment 
only for H. pylori patients with peptic ulcer 
disease (2). Because asymptomatic infection is 
very common, treatment of all asymptomatic 
persons would be expensive and might promote 
antibiotic resistance. Prophylactic and therapeu- 
tic vaccines against H. pylori are being developed 
(3,4); however, if vaccines are to be cost-effective, 
companies must take into account the changing 
epidemiology of H. pylori and related diseases. 
Knowledgeof these trends can also allow health 
agencies to predict resource allocations needed 
for these diseases. 

In industrialized countries, H. pylori incidence 
has been decreasing in successive generations, 
without any targeted intervention (5-7). 


Address for correspondence: Marcia F.T. Rupnow, Health 
Research and Policy T221, Stanford, CA 94305-5405, USA; fax: 
650-725-6951; e-mail: marcia.rupnowc^tanfordal umni.org. 


Quantifying the benefits of intervention in 
reducing disease incidence or cost requires an 
analytical model that estimates the natural 
course of H. pylori and associated diseases. The 
model could then estimate the decrease in 
H. pylori incidence that would result from an 
intervention strategy, relative to the natural 
course of infection. 

Wepresent an analytical frameworkto model 
H. pylori transmission dynamics and its 
subsequent disease progression. Our goal is to 
estimatefuturetrendsof H. pylori and associated 
diseases in the United States based solely on its 
natural history (i.e., without intervention). 

Methods 

Model Description and 
Epidemiologic Assumptions 

Our analysis is based on a dynamic 
compartmental model, a technique in which the 
population is divided into compartments and 
mathematical equations describe the transfer of 
persons from one compartment to another (8). I n 
thistechnique,theincidenceof H. pylori infection 
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is calculated as a function of the number of 
susceptible and infected persons in the popula- 
tion and a constant called the transmission 
parameter (15). We sought transmission param- 
eters that would best explain historical trends 
and allow estimates of future trends of H. pylori 
prevalence and the incidence of H. pylori- 


associated gastric cancer (GC) and duodenal 
ulcer (DU). 

We developed a compartmental model that 
captures the age-dependence of H. pylori 
infection and disease progression in infected 
persons (Figure 1)(6,9-17). The population was 
compartmentalized according to three factors: 


bhrth 





15 + 

yrs 

(Adult) 


I: not-suKceplible; 

S: sunrcptibEc; 

AG: antnim-pTiefJomLn^TU 
g^LHlntLH.; 

CG: cuq^uii-pfrdDrniiiant 
ga.siiiti.&^ 

DU: duodenal ulcer; 
CAG: ctirOnic utruphic: 

GC; gastric cancer 


X|T prcihabtlily of a Kuitcqnible acquiring Up Jind devclopiri^ AG; 
X^\ probabilily of a suscepiible acquiriug Up and developing CO; 
5 i; progression rate from AG to CG; 

Sii progreKion rate from AG to DU; 
prugrcssion rate from DU to CAG ; 
progression rare from CG lo CAG; 

65; progression rare from CAG to GC; 

Hixi i nwitality rale due to DU; 

mojlAlity mlc due to GU; 
pLij.),': murtaJity rate due lo GC 

ugc-Hpecifiu murtaJity rate from al] causes 


Figure 1. Compartmental model of Helicobacter pylori transmission and disease progression. I n this model, the 
population is divided into compartments according to age, infection status, and clinical state. Boxes represent 
population subgroups and arrows indicate transitions between subgroups, as well as flow into and out of the 
population (birth and death). 
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age (child [4 years old], youth [5 to 14 years old], 
and adult [>15 years old]); infection state 
(uninfected and infected); and clinical state 
(normal stomach, antrum-predominant gastritis, 
corpus-predominant gastritis, DU, chronic atro- 
phic gastritis [CAG], and GC) (Appendix). DU 
and GC states corresponded to persons whose 
illnesses were caused by H. pylori infection 
(H. pylori -associated DU and GC). We excluded 
from our model DU and GC not related to 
H . pylori. Computer simulation was used to solve 
the system of equations numerically and 
calculate the number of persons in each 
compartment over time. 

We modeled the possibility of being born 
susceptible or not susceptible to H. pylori 
(F igure 1). The nonsuscepti bleor "isol at ed" group 
comprised persons who could not become i nfected 
for physiologic, physical, or immunologic rea- 
sons. This conforms with epidemiologic studies 
indicating that, even within populations at high 
risk, a small percentage (10% to 20%) does not 
become infected with H. pylori (9,10). 


A susceptible person could become infected 
with H. pylori at any age and either antrum- or 
corpus-predominant gastritis could develop. The 
model represents the net flow from one clinical 
state to the other. Thus, we represented a 
positive flow from antrum gastritis to corpus 
gastritis becausethis progression is much larger 
than the reverse. From the gastritis states, we 
modeled the possibility of further clinical 
progression to DU, CAG, and GC. The model is 
consistent with previously reported outcomes in 
that persons with antrum-predominant gastritis 
are at higher risk for DU , whi le those who have 
corpus-predominant gastritis are at higher risk 
for CAG, GU, and GC (11,12-17). We also 
incorporated the possibility of CAG and GC 
developing in DU patients. 

Input Data and Sources 

Input assumptions are divided into disease 
and population parameters (T able). Weadopted a 
constant demographic characteristic of birth and 
death rates, to eliminate complexities associated 


Table. Input variables and sources 


Symbol 

Description 

Value 

Explanation 

Source 

Population parameters 




p (a) 

Death rate (per person per year) 

Age-specific 


18 

Z 

II 

Population size (persons) 
Birth rate (persons per year) 

200,000 

2,962 

Derived from disease-free 
equilibrium simulation 

Assumption 

Pi 

Disease 

Proportion (%) nonsuscepti ble 

parameters 

20 

Based on data from 
developing countries 

Assumption 

Pc 

Proportion (%) of children with 
antrum (vs. corpus) gastritis 

5 


Assumption 

Py 

Proportion (%) of youths with 
antrum (vs. corpus) gastritis 

75 

Based on simulation 

Assumption 

Pa 

Proportion (%) of adults with 
antrum (vs. corpus) gastritis 

95 


Assumption 

-5ic 

AG(- to CG(- (per person per year) 

0 


Assumption 

V 

AGy toCGy (per person per year) 

0 


Assumption 

^lA 

AG^ to CG^ (per person per year) 

0.0010 

1% progression in 10 yrs 

Assumption 

-52 

AG^ to DU (per person per year) 

0.0046 

14% progression in 32 yrs 

22 

_^3 

DU to CAG (per person per year) 

0.0020 

2% progression in 10 yrs 

15 


CG^ to CAG (per person per year) 

0.0112 

30% progression in 32 yrs 

22 


CAG to GC (per person per year) 

0.0030 

3% progress! on i n 10 yrs 

13 

Pdu 

Mortality rate due to DU 
(per person per year) 

0.0050 


23 

— Pgu 

Mortality rate due to GU 
(per person per year) 

0.0100 


24 

— Pgc 

M ortal i ty rate du e to G C 
(per person per year) 

0.3219 

80% death i n 5 yrs 

25 


GC = Gastric cancer; DU = Duodenal ulcer; CG = Corpus-predominant gastritis; AG = Antrum-predominant gastritis; 
CAG =Chronic atrophic gastritis. 
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with changing demographics (8). Baseline 
demographics were characteristic of the general 
U.S. population in 1950. The age-specific death 
rates from all causes were derived from the 
survival curve of the Vital Statistics Report (18). 
We assumed an initial population of 200,000, 
which corresponds to the population of a medium- 
sized city i n the U nited States (19). The birth rate 
was set so that the size of the population 
remained constant when H. pylori -associated 
diseases were not present. This birth rate was 
deri ved from a di sease-free si mul ation exerd se i n 
which we assumed zero infection with H . pylori. 

The proportion of isolated (nonsusceptible) 
persons in the population, p,, is an assumption 
based on data from developing countries, in 
which the prevalence of H. pylori among adults 
achieves an asymptote of approximately 80% by 
age 30 (1,20). We assumed that earlier 
acquisition of H. pylori is associated with 
development of corpus-predominant gastritis 
and subsequent CAG, GU, and GC, and that 
acquisition at older ages is more likely to lead to 
antrum-predominant gastritis and subsequent 
DU (11,12,21). Therefore, we adopted 5% and 
95% for the proportions of gastritis initially 
restricted to the antrum in children and adults, 
respectively. We chose 75% as the proportion of 
antrum gastritis i n youths; this esti mation had to 
be closer to that of adults than to that of children 
for the delay of Hp acquisition to result in more 
persons with antrum-predominant gastritis and 
subsequent DU. 

We derived the rates of clinical progression 
from gastritis to cancer from published informa- 
tion (13,15,22-25). The model also incorporated 
deaths from DU, GU, and GC. GU is not an 
explicit compartment in our model; instead, we 
modeled the possibility that a person with CAG 
would die of GU, in accordance with previous 
studies that documented an increased risk for 
GU in persons with CAG (12,16). 

Transmission Parameters, Vaiidation 
of the Modei, and Future Trends 

We calculated incidence as a function of the 
numberof susceptible and infected persons in the 
population and the transmission parameter, 15, 
which is a constant that characterizes infectivity 
of the pathogen. The number of infected persons 
at each period was determined by adding the 
number of persons with antrum- and corpus- 
predominant gastritis, DU, and CAG. We 


multiplied the number of CAG patients by a 
factor (<L), to reflect the fact that the infection 
may subside (or H . pylori density may be lower) 
as a result of atrophy (26), and therefore would 
contribute less to the transmission of H. pylori. 
Furthermore, we assumed that the role of GC 
patients as a source of infection was negligible. 

(5 may differ accordingto the population, age 
of i nfecti ve and suscepti ble persons, envi ron men- 
tal conditions, household sanitary and hygiene 
practices, and genetic characteristics of the 
organism itself (e.g., strain differences). Direct 
measurement of (5 is not possible for most 
infections (8). However, the value of 15 and its 
change over ti me must be esti mated to assess the 
impact of public health programs against H. pylori. 
Since the incidence and prevalence of H. pylori 
have been decreasi ng i n the absence of a vacci ne 
(which would decreasethe number of susceptible 
persons) or widespread eradication therapies 
(which would decrease the number of infected 
persons), the historical decrease in H. pylori 
incidence must represent a decrease in (5. Thus, 
weestimatedthevalueof Bin the 19th century by 
simulating an endemic equilibrium condition of 
H. pylori intheU.S. population. Wethen estimated 
the change in 15 that would account for the 
observed patterns of DU and GC. Finally, we 
extrapol ated (5 for the 21st century to esti mate the 
futuretrends in H. pylori -associated DU and GC. 

We set up the model with nine transmission 
parameters, according to the age groups of 
infected and susceptible persons. To find (5s in 
mid- 1800s, we assumed that the pattern of 
infection in the United States at that time was 
similar tothat in developing countries today, i.e., 
rapid acquisition of H . pylori at younger ages (up 
to 5 years old) and subsequent slower acquisition 
rate, achieving a maximum prevalence of 
approximately 80%. I n personal interviews, four 
experts in H. pylori epidemiology provided an 
assessment of the transmissibility among 
different age groups; they estimated transmis- 
sion among children to be 5 to 10 times higher 
than among adults (Hazel I SL, Megraud F, B laser 
M , Correa P, unpub. comm.). No assessment could 
be obtained for the different inter-age group 
transmissions ((5 (-y- (^ca- (^•c' ^a’ ^ac- ^ay^- 
thebaselineanalysis, weassum^that inter-age 
group transmissions are negligible compared 
with intra-age group transmissions (l5(-(~, (5yY, 
15^^) because on average, interaction among 
persons of the same age group is the highest. 
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Once we established transmission param- 
eters for the mid-1800s, we estimated changes in 
transmissibility that could explain the historical 
patterns of GC and DU. GC incidence and deaths 
have been decreasing si nee the beginning of the 
20th century (25,27,28). DU, however, is 
characterized by a rise and fall within the 19th 
century (29-31). For DU, long-term estimates 
were based on other statistics, such as deaths, 
physician visits, and hospitalized patients. 

Investigators have postulated that the 
different patterns of GC and DU could be 
explained by a more rapid decline in childhood 
than in adult incidence, thereby decreasing the 
proportion of new infections acquired in 
childhood (10,32). Weevaluated whether changes 
in 15 consistent with this hypothesis could 
reasonably explain disease patterns. We then 
extrapolated the transmissibility values intothe 
year 2100 on the basis of our projections for 
future socioeconomic changes that would affect 
the dynamics of H. pylori transmission, and 
estimated the future trends in H. pylori 
prevalence, as well as the incidence of H. pylori- 
associated GC and DU. 

Model Implementation 

We implemented the model using the 
software package Powers! m Constructor Version 
2.5 (Powers! m Corp., Herndon, VA), integrated 
with Microsoft Excel 97 (Microsoft Corp., 
Redmond, VA). We set Powers! m to use Runge- 
Kutta 4th-order integration method and time- 
step of 1 year. 

To obtain correct flows of age groups, at the 
implementation level we further subdivided the 
child and youth groups by 1-year age increments 
(i.e., birth through first birthday, 1 year of age 
through second birthday ... 14 years of age 
through 15th birthday). We subdivided the adult 
compartments by 5-year age increments up to 85 
years of age (15 years of age through 20th 
birthday ... 80 years of age through 85th 
birthday, 85 -f), which allowed us to program the 
age-specific death rate and track the aging of 
cohorts with relative accuracy. 

Results 

Transmissibility of H. pylori in 1850 

The values transmissibility most consistent 
with the endemic pattern of infection were: l5(-(- = 
0.000055, I5 yy = 0.000011, and 15^^ = 0.0000012. 


This translated into Potential Transmission 
Rates (PTR, defined as the average number of 
secondary infections per unit of time that an 
infected person would produce in an infection- 
free population) of PTR^-^. =0.6 for children (one 
infected child would be expected to transmit the 
pathogen to 0.6 susceptible child per year), 
PTRyy = 0.25 for youths, and PTR^^ = 0.15 for 
adults. In real populations, these transmission 
parameters would result in lower incidence rates, 
since no population is infection free. 

Change in Transmissibility from 1850 to 1995 

We estimated the change in transmissibility 
from 1850 to 1995 that was most consistent with 
the historical trends of DU and GC and with 
current H. pylori prevalence in the United States 
(Figure 2). The fastest decrease occurred in the 
latter half of the 19th century. Transmissibility 
then leveled off at =0.000042, (Syy =0.0000044, 
and 15^^ =0.0000009. Only a pronounced decrease 
in (5 at the end of the 19th century could explain 
the decline in GC incidence and deaths observed. 
Such a plateau in transmissibility was necessary 



Figure 2. Temporal change in transmission param- 
eters. The transmissibility values in 1850 yielded a 
pattern of infection similar to that observed in 
developing countries today (rapid acquisition in 
younger ages and lower acquisition in older ages). The 
rapid decline of transmissibility in the latter half of 
the 19th century is consistent with GC and DU 
patterns. Although the graph does not show the 
decrease in (5^^ because of the scale ((5^^ is much 
smaller than (5^^ and (5 yy), (5^^ also decreased over time. 
Constant extrapolation of transmissibility assumes 
no change in standards of living that would affect 
H. pylori transmission. 
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fortheprevalenceof H. pylori infect ion to coincide 
with current rates. An S-shaped Gompertz curve 
(33)— in which the decrease is slow initially, 
accelerates, and slows down again, trending to a 
stable level— best illustrates the past trends in 
H. pylori and associated DU and GC. 

H. pylori Prevalence, DU Incidence, and GC 
Incidence 

Wesimulated the temporal trend of H . pylori 
prevalence in the general U.S. population by age 
category (Figure 3a). Among children, prevalence 
decreased from 30% in 1850 to 1% by the end of 
the 20th century; among youths, prevalence 
decreased from approximately 70% to 5%, and 
among adults, from approximately 80% to less 
than 20% in the same period. We compared the 
simulation outputs with available U.S. data on 
incidenceof GC and DU (Figure 3b). According to 
the two sources on GC incidence and deaths 
(25,28), total GC incidence declined from approxi- 
mate! y 34 per 100,000 i n 1930 to 6.6 per 100,000 i n 
1995. In comparison, FI. pylori -associated GC, 
estimated from the simulation, decreased from 24 
per 100,000 in 1930 to 5.8 per 100,000 in 1995. 
Kurata and colleagues, based on data from a large 
health-maintenance organization in southern 
California from 1977 to 1980, estimated the 
incidence of DU to be 58 per 100,000 (34). In 
comparison, our model estimated FI. pylori-asso- 
dated DU incidencetobe47 per 100,000 in 1980. 

To extend the simulation through 2100, we 
assumed that the transmissibility of FI. pylori 
would remain stable at the level in 1995 
(Figure 2) and kept the input parameters 
constant. FI. pylori prevalence, as well as the 
associated DU and GC, would continue to 
decrease in the 21st century. By 2100, the model 
predicted incidences of 1.3 and 12.2 per 100,000 
population for FI. pylori -associated DU and GC, 
respectively, with overall FI. pylori prevalence 
decreasing to 4.2%. 

Sensitivity Analyses 

We tested the model by using different input 
assumptions. For each change in the input 
assumptions, we generated curves of FI. pylori 
prevalence and incidences of FI. pylori -associated 
GC and DU.^ The most sensitive variables were 
the proportion of persons with antrum- (vs. 



A Children Youths Adults 



Year 

— Hp GC (simulated) -m— Hp DU (simulated) 

-*-SEERGC -*-NCHSGC 

g -at- Kurata et al. DU (Ref. 38) 


Figure 3. Simulated temporal change in Helicobacter 
pylori and associated diseases, a) Prevalence of 
H. pylori in the U.S. population by age category 
(simulation output), b) Incidence of DU and GC per 
100,000 U.S. popuiation. "Hp GC" and "Hp DU" are 
simulation outputs. "SEER GC" represents incidence 
of GC from all cases, from the Surveillance, 
Epidemiology, and End Results report (25). "NCHS 
GC" represents incidence of GC from all cases, derived 
from the mortality statistics published by the 
National Center for Health Statistics and the Bureau 
of the Census (28). "Kurata et al. DU" represents 
incidence of DU from all cases, based on data from a 
large health maintenance organization in southern 
California (34). The simulation model predicted that 
the prevalence of H. pylori, as well as incidence of 
H. pylori -associated DU and GC, would continue to 
decrease in the 21st century. 

GC = gastric cancer; DU =duodenal ulcer; NCHS =National 
Center for Health Statistics; SEER = Surveillance, 
Epidemiology, and End Results system. 


^Specific graphs can be made available to readers upon request. 
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corpus-) predominant gastritis, rate of progres- 
sion from AG to DU , and rate of progression from 
CAG to GC. For example, when the rate of 
progression from CAG to GC was 0.0035, the 
incidence of H. pylori -associated GC decreased 
from 6.7 per 100,000 i n 1995 to 1.5 per 100,000 i n 
2100; when the rate was assumed to be 0.0025, 
the incidence in 1995 was 4.9 per 100,000, 
decreasing to 1.1 per 100,000 by the end of the 
21st century. Overall, we found that H. pylori 
preval ence wou I d decrease to <10% i n t he general 
U.S. population and the incidence of H. pylori- 
associated DU andGC would decreaseto <20 and 
1.5 per 100,000 population, respectively, by the 
end of the 21st century. 

Discussion 

This work provides insights into the intrinsic 
dynamics of H. pylori and associated diseases, 
applying and extending the basic principles of 
mathematical modeling used extensively in 
infectious disease analysis. First, we found that 
i n the U nited States, transmissi bi I ity of FI . pyl ori 
has a I ready become so low that the i ncidence and 
prevalence of the organism will continuously 
decrease in the foreseeable future, without any 
targeted intervention. The disappearance of 
FI. pylori, however, would take more than a 
century without intervention. Second, we found 
that only a rapid decline in FI. pylori 
transmissi bi I ity during the second half of the 
19th century could explain the rapid decrease in 
GC incidence observed in the 20th century. We 
tested different curves of decrease in FI. pylori 
transmi ssi bi I ity and found that a Gompertz curve 
best fits the past trends of FI. pylori and 
associated GC and DU. Gompertz curves have 
been used to descr i be deci i ne of I i vi ng organ i sms, 
as well as mechanical devices (33). Our findings 
suggest that those curves could also be used to 
represent waning of transmissi bility (infectivity 
or transmission potential) of an organism in a 
given host population. 

The decI i ne i n transmissi bi I i ty paral I el ed the 
rapid improvement in sanitation and hygiene 
that occurred i n the 19th and early 20th centuries 
in theU nited States. I n a historyof hygiene in the 
United States, Floy describes how the United 
States, once a country with poor sanitary 
conditions, became one with exemplary cleanli- 
ness (35). The author states that "in 1850, 
cleanlinessintheUnitedStates, north and south, 
rural and urban, stood at Third World levels." 


During the Civil War era, appreciation for 
cleanliness expanded. Frequent washing of 
clothes and bathing with soap, expanded 
availability of clean running water, and 
development of a sewage network i n the late 19th 
century were among improvements in personal 
hygiene. Although we cannot rule out other 
causes, this massive improvement in sanitary 
and hygienic conditions (better household 
hygiene and community sanitation projects, 
including water purification systems), which was 
responsible for reducing the incidence of 
infectious diseases such as shigellosis and 
typhoid, could also explain the decrease in 
transmissi bility of FI . pylori. 

This model attempted to explain quantita- 
tively how theshift in age of FI. pylori acquisition 
could have produced the different patterns of DU 
and GC outcomes. We cannot rule out other 
factors (such as change in FI. pylori strains), not 
included in our model, which could have been 
responsible for the rise and fall of DU and the 
parallel decrease in GC. 

We can indirectly estimate FI. pylori- 
associated GC from National Center for Flealth 
Statistics (NCFIS) and Surveillance, Epidemiol- 
ogy, and E nd Results (SE E R) curves showi ng GC 
incidence from all cases, based on attributable 
risk (Figure 3b). Since attributable risk for 
FI. pylori is believed to have decreased because 
the proportion of proximal GC is higher today 
than it was at the beginning of the 20th century, 
the total GC curve underestimates the decline in 
FI. pyl ori -associated GC. Thus, the incidence of 
FI. pyl ori -associated GC from our simulation is 
higher than the one estimated by the method of 
attributable risk applied to NCFI S/SEER data. 
This discrepancy can be explained by the 
existence of other factors acting in conjunction 
with FI. pylori topromotethedevelopmentof GC. 
The discrepancy, however, does not affect our 
conclusions about future decreasing trends in 
FI . pylori and associated diseases. 

In this study, we did not model the 
demographic changes— such as increase in birth 
rate and life expectancy— that occurred in the 
United States from 1850 to the present. 
Demographic changes affect the epidemiology 
and disease statistics in many ways. For 
example, an increase in birth rate would affect 
population density and crowding, depending on 
availability of land. Change in the size of a family 
unit could also affect thetransmission dynamics. 
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The model quickly becomes too complex and 
intractable. We therefore fixed the population 
parameters as in other infectious disease 
transmission models (8,36,37). Because our 
si mul ati on covers a I ong ti me span ( 1850 to 2100), 
to compromise between the past, present, and 
future, we derived the population characteristics 
fromthedemographicdataof 1950. 1 n sensitivity 
analyses, we also tested different birth and death 
rates, but the predictions for future H. pylori 
prevalence and incidence of H. pylori -associ- 
ated DU and GC did not differ significantly 
from the base case. 

I n our baseline model, we assumed that 20% 
of the population would not become infected with 
H. pylori. This is based on studies from 
developing countries, where prevalence in older 
adults is typically 70% to 90%. Given the high 
rates of transmission in childhood required to 
account for this high endemicity, the prevalence 
of nonsusceptible persons (because they are 
either immuneor physically isolated) iscritical in 
explaining why 100% prevalence has not been 
observed in any population. The model, however, 
can accommodate a range of values, including 
100% prevalence. 

Wesimplified thequantifi cation of transmis- 
sion parameters by considering only transmis- 
sion between persons in the same age group. I n 
reality, transmission can also occur between 
adults and children or adults and youths. 
However, no data are avail able on the magnitude 
of transmission between or among different age 
groups. Therefore, we also tested the model by 
using positive numbers for inter-age group 
transmission, but the predictions for future 
H. pylori prevalence and incidence of associated 
GC and DU did not differ significantly from our 
base case results. 

Given the model estimation that H. pylori 
transmissibility has remained relatively con- 
stant sincethebeginningof 20th century, despite 
many recent developments and changes in life- 
style, we assumed that maintenance of the 
current H. pylori transmissibility would be a 
reasonable extrapolation for the next century. 
Despite this assumption of constant transmissi- 
bility, the model predicted that H. pylori 
prevalence, as well as H. pylori -associated DU 
and GC, would continue to decrease without 
intervention. Thus, the decreasing trends 
indicate that transmissibility of H. pylori has 


already declined below the threshold value 
needed to maintain the organism endemic in the 
population. Without any targeted intervention, 
however, the disappearance of H. pylori from the 
U.S. population will take more than a century. 
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Appendix 


The mathematical equations underlying our compartmental model of H. pylori is a 
system of partial differential equations; 

^ ^ = -U (a, <)+ 'll (o. ') + M")] ■ •?(<>.') 

at da 

{‘^) + /) 

dt da 

^ + ^^^X,{a,tys(aj)-[SAa) + M{a)]-CG{a,t) 
at da 

^^ + -^^ = '52(o)/IC(a,()-['5j(a)+ft,i' +//(a)]D{/(o,;) 
dt da 

= S,{a)DU{a,l)+ S, (a) • CC(a, <) - fc («)+//<,„ + M] ■ C^C(fl, () 

dt da 

^ + ^ = S,(a)- CAG(a, l)-Luac + M] ' GC{a. t) 
dt da 

where; 

(a, /) — p{a)- Jy5(a',fl) \AG{a' j] + CG{a' ,A + DU{a' ,l) + a ■ CAG{a\t)-\-GC{a' ,ty da' 

0 

X 2 {a,t) - (1 - p{a))- Jy5(<j',a)- \AG{a' ,t) + CC{a' d) + XXU{a' ,t] + a ■ CAG{a' ,t) + GC{a' d'^da' 


The boundary conditions are; 

i{Q.t)=p, n 

x(o,()=(i-ft)n 

AG{0, t) = CC(0, l) = DU{0,l) = CAGiO, t) = GC(0,t) = 0 
Notation; 

a, a ’ age index 

t time index 

77 birth rate per unit time 

l(a,t) number of isolated (not-susceptible) individuals of age a, at time t 
S(a,t) number of susceptible individuals of age cj, at time t 

AG(a,t) number of infected individuals of age a with antrum-predominant gastritis, at 
time t 

CG(a,t) number of infected individuals of age a with corpus-predominant gastritis, at 
time t 

DU(a,t) number of individuals of age a with duodenal ulcer, at time t 
CfGffl.tlnumber of individuals of age a with chronic atrophic gastritis, at time t 
GC(a,t) number of individuals of age a with gastric cancer, at time t 
Pi proportion of population that is not-susceptible at birth 

Xi(a.t) rate at which one susceptible of age a acquire infection and develop antrum- 
predominant gastritis 

X 2 (a,t) rate at which one susceptible of age a acquire infection and develop corpus- 
predominant gastritis 

P(a \a) transmission parameter; probability that an infective of age a ' will infect a 
susceptible of age a 

p(a) proportion of newly infected individuals of age a developing antrum (vs. 
corpus) predominant gastritis 

Si(a) transition rate from antrum- to corpus-predominant gastritis in age group a 
52 (a) progression rate from antrum-predominant gastritis to duodenal ulcer in age 
group a 

Ss(a) transition rate from duodenal ulcer to chronic atrophic gastritis in age group a 

84 (a) progression rate from corpus-predominant gastritis to chronic atrophic gastritis 
in age group a 

Ss{a) progression rate from chronic atrophic gastritis to gastric cancer in age group a 

p(a) age-specific background mortality rate due to all cases 

Pov mortality rate due to duodenal ulcer 

Pqu mortality rate due to gastric ulcer 

Pqc mortality rate due to gastric cancer 
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Using Remotely Sensed Data 
To Identify Areas at Risk 
For Hantavirus Pulmonary Syndrome 
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The 1993 U.S. hantavirus pulmonary syndrome (HPS) outbreak was attributed to 
environmental conditions and increased rodent populations caused by unusual weather 
in 1 991 -92. In a case-control study to test this hypothesis, we estimated precipitation at 
28 HPS and 1 70 control sites during the springs of 1 992 and 1 993 and compared it with 
precipitation during the previous 6 years by using rainfall patterns at 196 weather 
stations. We also used elevation data and Landsat Thematic Mapper satellite imagery 
collected the year before the outbreak to estimate HPS risk by logistic regression 
analysis. Rainfall at case sites was not higher during 1992-93 than in previous years. 
However, elevation, as well as satellite data, showed association between 
environmental conditions and HPS risk the following year. Repeated analysis using 
satellite imagery from 1 995 showed substantial decrease in medium- to high-risk areas. 
Only one case of HPS was identified in 1996. 


In 1993, a disease characterized by acute 
respiratory distress with a high death rate 
(>50%) among previously healthy persons was 
identified in the southwestern United States. 
This disease, hantavirus pulmonary syndrome 
(HPS), was traced to infection with an 
unrecognized, directly transmissible virus— 
Sin Nombre virus (SNV) (Bunyaviridae; 
Hantavirus) (1). The virus was maintained and 
transmitted primarily within populations of a 
common native rodent, the deer mouse 
(Peromyscus maniculatus) (2), and transmission 
to humans occurred through contact with 
secretions and excretions of infected mice (3). 

It has been hypothesized that the El Nino 
Southern Oscillation (ENSO) of 1991-92 was the 


Address for correspondence: Gregory E. Gurri Glass, 
Molecular Microbiology & Immunology, TheJ ohns Hopkins 
School of Public Health, 615 N. Wolfe Street, Baltimore, MD 
21205, USA; fax: 410-955-0105; e-mall: gurrlglasscgmsn.com. 


major climatic factor producing environmental 
conditions leading to the outbreak of HPS in 
1993. Unseasonable rains in 1991 and 1992 
during the usually dry spring and summer and 
the mild winter of 1992 are thought to have 
created favorable conditions for an increase in 
local rodent populations (4,5). 

This hypothesis is based primarily on the 
following observations: 1) ENSO tends to 
influence the timing and abundance of precipita- 
tion in the southwestern United States. 2) Some 
Peromyscus populations increased dramatically 
in areas where precipitation was above average 
but remained near normal levels where 
precipitation did not increase— this observation 
is based on comparison of data from only two 
study areas: the University of New Mexico's 
Sevilleta Long-Term Ecological Research station 
(90 km south of Albuquerque, NM), where 
precipitation was 2 to 3 times the previous 
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20-year average, and Moab, Utah, where 
precipitation was at or below normal in the 
summer of 1992. Before and during the HPS 
outbreak, populations of P. maniculatus did not 
increase in Moab, while at sites on the Long- 
Term Ecological Research station they were 10 to 
15 times higher than normal. Moreover, both 
Moab and the research station were 200 km to 
300 km from the epicenter of the 1993 HPS 
outbreak, making comparison with conditions 
where disease occurred uncertain. 3) 1993 case 
studies found that rodent abundance varied 
dramatically over short distances. Rodent 
populations were higher at HPS households than 
at neighboring households without disease or 
randomly selected households at least 25 km 
away (6); however, these one-time case studies 
provide little information on the responses of 
rodent populations to changes in environmental 
conditions. Although associating weather with 
H PS outbreaks is consistent with these observa- 
tions, supporting data are limited. This reflects 
the situation for many emerging diseases, as 
active surveillance for unidentified diseases is rare. 

Current attempts to understand the factors 
leading to HPS outbreaks focus on detailing the 
chai n of events from weather, through changes i n 
vegetation, to virus maintenance and transmis- 
sion within rodent populations, culminating in 
changes in human disease risk (trophic cascade 
hypothesis) (4,5,7). An impediment to this 
approach, however, is the focal nature of SNV 
within local populations of P. maniculatus. 
Among local rodent populations the rates of 
infection vary, and some populations appear 
uninfected (8). The reason for this is uncertain 
but may be related to the stochastic loss of the 
horizontally transmitted viruses within local 
populations of the reservoir (9). Alternatively, 
the dynamics of local populations may be such 
that SNV cannot be maintained at very low 
population numbers. Under either circumstance, 
responses of local populations to environmental 
fluctuationscouldsubstantially alter human risk. 
Given the sporadic natureof HPS outbreaks, ongoing, 
longitudinal monitoringof rodent populations in the 
vicinity of subsequent cases is unlikely (6). 
Consequently, inferring human disease risk from 
rodent-SNV population dynamics will again 
require extrapolating from studies in regions 
other than the site of the outbreak. 

I nthisarticle, weexaminetherelationshipof 
theenvironmenttoH PS risk by using locations of 


HPS cases as sites where people were associated 
with infected rodents. This approach avoids the 
immediate need to establish the conditions that 
lead some reservoir populations to be uninfected 
by SNV. We compare the environmental 
characteristics of sites where people were 
i nf ected w i t h t h ose atsiteswhere peop I e were not 
infected. Differences in environmental condi- 
tions could indicate factors that influence either 
the abundance of rodents or the occurrence of 
virus, creating testable hypotheses of environ- 
mental conditions that influence SNV infection 
patterns in reservoir populations. As a partial 
test of our method, the analysis was repeated 
when cases of H PS were uncommon. U nder these 
conditions, the identified environmental factors 
should indicate low levels of disease risk. 

Two sources of information were used as 
measures of the local environmental conditions 
precedi ngtheH PS outbreak: 1) Monthly patterns 
of precipitation from March toj une (generated 
from archived weather station data in the region 
to estimate local precipitation patterns at both 
case and control sites) and 2) satellite imagery 
(obtained before the H PS outbreak and used as a 
measure of variable, local environmental condi- 
tions and HPS risk evaluated by epidemiologic 
analysis). 

Methods 

Study Population and Region 

The epidemiologic analysis was performed as 
a case-control study. Twenty-eight (93.3%) of 30 
sites with confirmed cases of HPS identified in 
the region between November 1992 (identified 
retrospectively) and November 1994 were 
selected. I nclusion criteria were based on clinical 
disease consistent with the Centers for Disease 
Control and Prevention case definition that was 
confirmed by serologic, nucleic acid, or immuno- 
histochemical tests (1). One case was excluded 
because the likely site of exposure could not be 
established, and a second case because we could 
not confi rm that the proper geographic location of 
the site was recorded during data collection. 

S i tes of exposure were est abl i shed previ ousi y 
by investigation of each case-patient's activities 
and, for most fatal cases, demonstration of 
sequence homology of polymerase chai n reaction 
(PCR)-amplified regions of SNV nucleic acids 
obtained from case-patients and rodents col- 
lected at the imputed sites of exposure (1,2). Sites 
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of exposure for the HPS cases were at or in the 
immediate vicinity of households, and there 
usually was a history of activities that could have 
generated exposure to contaminated aerosols (1). 

T 0 control for issues related to access to care 
and socioeconomic conditions, controls were 
selected randomly from all households that used 
the same health clinics as the HPS patients 
during the same period as the HPS patients. 
Controls were randomly selected among persons 
without HPS. A total of 170 persons with 
different residential addresses were identified 
from a listing of visits to all the clinics during the 
time of the HPS outbreak. This represented an 
approxi matel y 2% random sampi e of addresses of 
the patient population. A previous study showed 
that subclinical infection with SNV was not 
observed among controls (10). Geographic 
locations of case and control sites were 
established by using global positioning system 
(GPS) receivers to record latitude and longitude 
at each site. 

The study area of 105,200 km^ was located i n 
the southwestern United States, incorporating 
the region of the original HPS outbreak (1,2). 
Epidemiologic surveillance was part of the H PS 
outbreak investigation. 

Environmental Characterization of Sites 

Monthly precipitation data recorded at 196 
weather stations throughout the region from 
1986tol993werefromtheU .S. National Oceanic 
and Atmospheric Administration's National 
Climatic Data Center. Spring precipitation at the 
weather stations was calculated by aggregating 
monthly precipitation for March through J une. 
Spring precipitation at individual case and 
control sites was estimated by interpolation 
among the nearest weather stations to each case 
or control site. I nterpolation procedures used the 
eight nearest weather stations to estimate 
precipitation by applying trend surface algo- 
rithms from the GIS software (IDRISI, 11). 
Spring precipitation during 1992 to 1993 at each 
site was compared with spring precipitation 
during previ ous years by pai red differences tests, 
with correction for multiple tests. The null 
hypothesis for each case or control site was that 
there was no increase in precipitation duringthe 
spring of 1992 to 1993 compared with the 
previous years. We al so used satel I ite i magery to 
develop detailed characterization of local envi- 
ronmental conditions. Weselectedthreearchived 


Landsat Thematic Mapper (TM) images origi- 
nally recorded in mid-J une 1992 for analysis. TM 
records digital numbers (DN ) from reflected I ight 
in six bands, threein visibleand threein infrared 
(IR) portions of the electromagnetic (EM) 
spectrum. (An additional band of thermal I R 
energy also is recorded but was not used.) 

The TM images were merged to form the 
study area of 105,200 km.^ Nominal pixel 
resolution was 30 m. I mages were geometrically 
and radiometrically corrected by the U.S. 
Geological Survey (USGS) Earth Resources 
Observation Systems (EROS) Data Center. The 
images were imported into a raster-based 
geographic information system (GIS) for geo- 
graphic registration by using control points 
obtai ned from USGS quadrangle maps (11). The 
images were resampled by a bilinear resampling 
procedure with a quadratic mapping function 
(11). Corrections for atmospheric scattering of 
bands one to three were performed by the method 
of Chavez (12). Latitude and longitude positions 
of case and control households were imported as 
a data layer in IDRISI GIS. Additional 
environmental variables incorporated in GIS 
included elevation, slope, and aspect of the case 
and control sites. Elevation was derived from the 
USGS digital elevation models (1:250,000 scale), 
while si ope and aspect were generated from these 
data by using software in Gl S. 

Three satellite images from mid-J une 1995 
were selected to further validate the analysis of 
the 1992 i magery. There was no E N SO from 1994 
to 1996, and we predicted that any relationship 
between the satellite imagery in 1992 and the 
case and control sites would indicate reduced risk 
in 1995. DNs for selected locations on the 1992 
and 1995 images were compared to determine 
any significant changes in sensor calibration. 

Epidemiologic Analysis of 
HPS by Satellite Imagery 

The spatial distribution of HPS sites (cases) 
was compared with that of control sites to 
determine if cases were spatially aggregated 
within the study region (13). Then, the 
relationship between HPS and environmental 
factors measured byTM imagery was examined. 
Becausethe analysis used threeti led images, the 
strategy for model development involved model 
fitting by using a portion of the study region, 
followed by external validation (14). A sample of 
HPS sites and control sites was selected for 
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analysis, by using logistic regression to examine 
the relationship between the odds of a site being 
an HPS site (outcome variable) and the DN in 
each TM band, elevation, slope, and aspect 
(predictor variables). The analysis was then 
repeated by using the remaining HPS and control 
sites to validate the model. Identifying the same 
model in thetwoanalyses would indicatethatthe 
HPS risk model was robust for that period. 

The initial model used a test area of 12,279 
km^ from the east-central region of the study 
area. This area included 14 case and 36 control 
sites and was entirely within oneTM image. The 
validation analysis used the remaining sites (14 
case and 134 control sites) and i ncorporated al I or 
parts of the three TM images covering 92,921 
km.^TheaverageDN for each of the si xTM bands 
used in the analysis was calculated with a 3 x 3 
pixel filter centered on the location for each case 
or control site. This sampled a local region of 
approximately 8,100 m^ around each case or 
control site. 

We used logistic regression analysis to 
identify the best combination of TM bands and 
environmental variables associated with HPS 
status (14). Elevation was dichotomized at the 
median elevation for case and control sites, 
combined. I nclusion of remotely sensed variables 
was based on the statistical significance of their 
coefficients in the model. Elevation was retained 
in all models because of its observed association 
with P. maniculatus abundance (8). The logistic 
model was evaluated with the Deviance and the 
Hosmer-Lemeshow goodness of fit statistics (C) 
for deciles of risk. To examine the accuracy of the 
model predictions, we evaluated the sensitivity 
and specificity by creating a Receiver Operator 
Characteristic function (15). The function 
compares the true-positive rate (sensitivity) 
agai nst the false-positive rate (1 - specificity) of a 
model by using various predicted values as 
thresholds identifying case and control sites. I n 
addition to examining spring precipitation, we 
evaluated the trophic cascade hypothesis by 
examining the relationship between HPS risk 
and vegetation growth before the H PS outbreak 
(5), regardless of habitat type. We used data from 
the near infrared (band 4) and red (band 3) 
portions of the spectrum to generate a 
normalized difference vegetation index of the 
region. The index is a standard algorithm used as 
a measure of vegetative growth. We compared 
the index and the TM bands identified in the 


initial epidemiologicanalysisfor estimating H PS 
risk by comparing the Receiver Operator 
Characteristic's generated by each analysis. 

Results 

Cases and controls occupied the same 
general, geographic area with the greatest 
difference of most extreme sites of 20 km in the 
north-south and 13 km in east-west directions. 
Despite the broad geographic overlap, cases of 
HPS were not randomly distributed within the 
area. HPS cases were spatially clustered (Table 
1). Despitethisclustering, H PS sites were widely 
separated geographically. The average distance 
between case sites was 50.3 km (SD =23.8 km), 
and the nearest neighbor! ng sites (k =l;Tablel) 
were not themselves likely to be case sites. 


Table 1 . Spatial clustering of households with hantavirus 
pulmonary syndrome 


k 

T[k] 

z 

P 

1 

6 

0.91 

0.181 

2 

14 

1.88 

0.030 

3 

24 

3.05 

0.001 

4 

33 

3.71 

<0.001 

5 

40 

3.87 

<0.001 

6 

50 

4.57 

<0.001 

7 

64 

5.85 

<0.001 

8 

75 

6.50 

<0.001 

9 

79 

6.14 

<0.001 

10 

86 

6.21 

<0.001 


Cuzick and Edwards' test (13), T[k], its z score under the null 
hypothesis and the p value using Simes correction for 
multiple tests (16). Results show significant spatial 
clustering among case households for the second through 
tenth nearest neighbors (k=2 10). 

Spring precipitation patterns showed sub- 
stantial interannual variation at case and control 
sites (Figure 1). From 1986 through 1993, 
precipitation was 4.5 mm (1989) to 110 mm 
(1992). Spring precipitation decreased markedly 
between 1992 and 1993 (Figure 1). Overall, 
precipitation at control sites tended to be lower 
(65 mm) than at case sites (72 mm) during the 
spring each year, but there was broad overlap 
among sites, and yearly variation at control sites 
tracked that of the case sites (Figure 1). None of 
the case sites had higher precipitation during 
1992 to 1993 than during the preceding 6 years 
(p >0.05). 

There was a significant relationship between 
local environmental conditions and HPS risk as 
measured by the statistical association between 


Vol. 6, No. 3, May-]une 2000 


241 


Emerging Infectious Diseases 



Research 


March - June Precipitation 



Time (Years) 


Figure 1. March-J une precipitation patterns at case 
sites (solid symbols) and control sites (open symbols) 
from 1986 through 1993. Vertical bars are 1 standard 
deviation in precipitation values. 


theDN recorded by the satellite in 1992 and H PS 
risk the following year. The logistic regression 
analysis developed for the training area fit the 
observed data well (Deviance =45.45; p =0.49, 
df =46; Hosmer Lemeshow C =6.62; p =0.58, 
df = 8), and none of the sites were obvious 
outliers. Higher level interactions between the 
i ndependent vari abl es did not change the results. 
Three of the six bands from the TM images (1, 5, 
and 7) were associated with the odds of HPS 


(Table 2). Sites abovethe median elevation (2,094 
m) were marginally associated with risk in the 
training area, but el evati on was retai ned because 
of the relatively small sample sizes used to 
estimate the parameters, as well as the biologic 
rationale outlined elsewhere (8). 

H igh DN values in the blue (band 1) and mid- 
infrared (band 7) portions of the EM spectrum 
were associated with decreased risk for HPS, 
while high values in the mid-infrared (band 5) 
portion of thespectrum werea risk factor for H PS 
(T able 2). The DN values were approximately 58 
to 233 units. Each unit change in the average DN 
around sites altered the odds of H PS risk by 6% 
(band 1), 15% (band 7). Sites above 2,094 m i n the 
test area were >4 times as likely to be H PS sites 
as sites below 2,094 m. Slope and aspect at the 
sites were not associated with risk. 

The Receiver Operator Characteristic graph 
of sensitivity and specificity of the predictor 
function showed that at least 95% of the case sites 
were correctly identified until the proportion of 
control sites correctly identified exceeded 56% 
(20 of 36 control sites) (Figure 2). The same 
predictor variables from the TM imagery were 
identified when the analysis was repeated for the 
validation area. The coefficients of the validation 
analysis did not differ significantly from those for 
thetrainingarea (Table2). This model alsofitthe 
data well (Deviance =74.49; p = 1.00, df = 144; 
Hosmer Lemeshow C =6.78; p =0.56, df =8). 

Because the logistic models did not differ 
between thetrainingand validation areas, all the 


Table 2. Logistic regression analysis of LANDSAT Thematic Mapper data and elevation 


Predictor 


Training area 



Validation area 


Odds 

ratio 

95% Cl 

P 

Odds 

ratio 

95% Cl 

P 

t 

Band 1 

0.94 

(0.89, 0.99) 

0.027 

0.95 

(0.90, 0.99) 

0.019 

ns 

Band 5 

1.14 

(1.03, 1.26) 

0.013 

1.08 

(1.01, 1.17) 

0.034 

ns 

Band 7 

0.85 

(0.73, 0.99) 

0.039 

0.90 

(0.80, 1.00) 

0.046 

ns 

Elevation 

4.63 

(0.88, 24.29) 

0.070 

2.56 

(0.77, 8.43) 

0.123 

ns 



Overall 


Controls 

Cases 


Odds 







Predictor 

ratio 

95% Cl 

P 

Mean 

Range 

Mean 

(Range) 

Band 1 

0.96 

(0.93, 0.99) 

0.012 

111.9 

(70.6, 164.0) 

97.8 

(72.0, 134.3) 

Band 5 

1.07 

(1.01, 1.13) 

0.013 

170.6 

(110.7, 253.3) 

160.7 

(110.2,226.2) 

Band 7 

0.91 

(0.84, 0.99) 

0.015 

102.6 

(60.6, 163.3) 

93.1 

(58.2, 140.4) 

Elevation 

4.60 

(1.96, 10.79) 

0.005 

1935 

(1465, 2500) 

2081 

(1825, 2341) 


Odds ratios, 95% confidence intervals, and statistical significance of HPS predictor variables. Odds ratios of the Thematic 
Mapper bands indicate thechange in risk for each unit change in digital number recorded by the satellite sensor. Elevation was 
divided as either above or below 2,094 m for the analysis. There was no difference between theT raining and Validation areas, 
as tested by t-test (t). The Controls and Cases columns show the mean and range of values at the control and case sites, 
respectively. 
ns=not shown. 


Emerging Infectious Diseases 


242 


Vol. 6, No. 3, May-]une 2000 




Research 



Figure 2. Receiver Operator Characteristic (ROC) 
function of the logistic model for HPS risk in the 
training area (closed squares) and overall for 1992 
(open diamonds) as the threshold for predicted case 
households was varied from p =0.00 - 0.85 in 0.05 
increments. The ROC for normalized difference 
vegetation index model (triangles) had p values of 0.00 
to 0.40 and multiple points occurred together. The 
early rapid loss of sensitivity for the N DVI model was 
the result of poor model specification. 


sites were combined to give an overall model 
(Table 2). The overall Receiver Operator 
Characteristic (Figure 2) had a sensitivity and 
specificity similar to those of the test area (95% 
sensitivity, 62% specificity). This threshold 
corresponded to a predicted value for HPS of 
approximately 0.10. Thus, using a predicted log 
odds ratio of at least 0.10 as a threshold for 
increased risk included 95%of the case sites and 
excluded 62% of the control sites. 

The logisticfunction was applied toeach pixel 
in the study area to produce a map of predicted 
risk (Figures) (17). The analysis was repeated by 
using satellite images from J une 1995 to predict 
HPS risk for 1996 (Figure 3). There was a near 
elimination of predicted high-risk areas in the 
1995 imagery and a broad expansion of low-risk 
areas compared with the 1992 images. The single 
case of HPS reported from the region in 1996 
occurred at a site with a predicted risk of 0.16 
(i.e., above the HPS threshold). 

Areas of high vegetation growth in 1992, as 
measured by the normalized difference vegeta- 



Figure 3. Comparison of predicted HPS risk for 1993 
(top) and 1996 (bottom) by satellite imagery taken in 
1992 and 1995, respectively, in the study area. Low- 
risk areas are in dark blue and high-risk areas are in 
red and yellow. There was a significant reduction in 
predicted high-risk areas in 1996 compared with 1993. 


tion index (Figure4) incorporated only a portion 
of the HPS high-risk areas (Figure 3). We 
evaluated vegetation growth, as measured by the 
index, asa predictor of HPS risk by modeling the 
case-control data, using the index and elevation 
as predictor variables. The vegetation growth 
model that best fit the observed data included an 
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Figure 4. The normalized difference vegetation index 
(N DVI ) scores of the study area by Thematic Mapping 
bands 3 and 4. Vegetation growth increased from 
brown through yellow to green. There was a 
substantial portion of high-risk area (especially the 
eastern portion of the image) where the N DVI image 
pixels did not obviously correspond to high-risk areas 
(see Figure 3 for comparison). 


exponential transformation of the normalized 
difference vegetation index and sites above the 
med i a n el evat i on . T h i s i n dex model accou n ted for 
a significant part of the variation in the H PS risk 
(deviance = 147.62; p =0.99, df =196) but did not 
accurately model the odds of HPS. In this 
analysis 11 (39.3%) of 28 case sites had 
standardized residuals exceed! ng three standard 
deviations, suggesting a poor fit to the data— an 
interpretation supported by the Hosmer 
Lemeshow statistic (C =20.09; p =0.01, df =8), 
which indicatedthattheformof risk model fitthe 
data poorly. The receiver operator characteristic 
of the vegetation index analysis (Figure 2) also 
lost sensitivity more rapidly than the analysis 
using TM bands 1, 5, and 7, especially over the 
range of values near the threshold of increasing 
HPS risk. 

Conclusions 

Satellite imagery, combined with epidemio- 
logic surveillance, retrospectively identified 
areas at high risk for HPS associated with 


Peromyscus populations over broad geographic 
regions during the 1993 outbreak. TM data 
identified environmental conditions near HPS 
sites that were measurably different from 
conditions in rural, populated sites wheredisease 
did not occur for nearly 1 year before the 
outbreak. These envi ron mental conditions varied 
with the presence of ENSO (Figure 3). The 
geographic extent and general level of predicted 
H PS risk were higher during ENSO, supporting 
theviewthat El N i ho may increasethe likelihood 
of HPS outbreaks. The hypothesized pathway 
between ENSO, increased spring precipitation 
leading to increased vegetation growth, and 
subsequent H PS risk, however, was not strongly 
supported by the data. Possible reasons for this 
lack of support are discussed below. 

In this study, we used a retrospective 
epidemiologic approach to risk assessment (15). 
Therefore, odds ratios of the environmental 
characteristics were used to estimate the 
population's relative risk for H PS. This approach 
is valid when the cases used in the study are 
representative of all cases, the controls are 
representative of the general population, and the 
di sease i s rel at i vel y rare. U nder these condi ti ons, 
odds ratios approximate relative risk (15). 

HPS is rare-fewer than 1,000 cases have 
been identified in North America, although most 
occurred in the southwestern region of the 
United States. In this study, the cases included 
nearly all (28 of 30) sites i n the region where HPS 
occurred during the outbreak. Although bias 
induced from this factor isunlikely, we cannot be 
certain that environmental factors identified 
with outbreaks of H PS are similar to those with 
the sporadic, single cases of HPS reported each 
year. However, the accurate identification of the 
sitewherethesinglecaseof H PS was observed in 
1996 suggests that the classification may also 
identify risk characteristics for this group. I ssues 
related to personal privacy make accessing these 
data difficult, however, because geographic 
locations of residences, for example, are 
considered personal identifiers. 

The selection of controls focused on a 
population from the same socioeconomic and 
geographic region as the HPS cases. Although 
HPS cases were clustered, the maximal 
geographic extent of both cases and controls was 
similar, suggesting that the enrollment proce- 
dure adequately fulfilled this goal. Random 
selection of controls also was intended to control 
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for access to care in a region where travel may be 
difficult. Restricting controls to those using the 
same health-care facilities again raises the issue 
of the applicability of the results to areas with 
different socioeconomic and cultural conditions 
and probably excluded much of the population 
within urban areas of the study site. 

These are relatively minor potential prob- 
lems. H PS cases are rare i n urban areas because 
the primary reservoir species in North America 
rarely occur within urban settings. Moreover, 
surveys of rural housing show that infestations 
by Peromyscus are nearly ubiquitous i n theabsence 
of focused rodent exclusion methods (6,18). 

Absence of a significant difference in spring 
preci pitati on at case sites during 1992-93 and the 
previous 6 years (Figure 1) may reflect either the 
absence of an effect (contradicting the trophic 
cascade hypothesis) or practical problems with 
estimating precipitation and incorporating con- 
ditions associated with past HPS outbreaks. 
Although nearly 200 weather stations were used 
in estimating spring precipitation at case and 
control sites, this still represents a relatively 
sparse network of sampling locations; therefore, 
localized precipitation could have been attoofine 
a spatial scale to detect. However, when we 
estimated precipitation at weather stations by 
using the surrounding stations and comparing 
the results with the observed precipitation, we 
found no difference between observed and 
predicted results in 1992-93. A more likely 
possibility is that the relatively short time series 
of precipitation data used makes demonstrating 
astatistical effect difficult. Additionally, if ENSO 
isatriggeringevent, outbreaksof H PS must have 
occurred in the past. Therefore, previous ENSO 
events may "contaminate" comparisons with past 
precipitation data because they include unrecog- 
nized H PS outbreaks. 

The trophic cascade hypothesis predicts that 
ENSO leads to increased precipitation that 
affects vegetation growth, subsequently influenc- 
ing HPS risk. The association between HPS risk 
and vegetation growth, as measured by the 
normalized difference vegetation index, was 
inconsistent. Areas with a high index (Figure 4) 
did correspond to areas at highest risk for HPS 
(Figure 3) in 1992. Similarly, areas at low risk 
were generally those with low normalized 
difference vegetation indexes. However, broad 
regions of moderate- to high-risk areas did not 


relate to the vegetation index, and the logistic 
regression model did not perform well (Figure 2). 

The failure of the normalized difference 
vegetation index to predict HPS risk may 
indicate that the ecologic connections hypoth- 
esized by the trophic cascade hypothesis are 
complex and modulated by intervening ecologic 
variables. Alternatively, the normalized differ- 
ence vegetation index, which is the normalized 
difference of DN in red and near-infrared portion 
of the EM spectrum, may have difficulty 
accurately characterizing vegetation growth in 
semi-arid areas that contain a complex mixtureof 
vegetation and bare ground (19). Further, 
detailed studies incorporating "ground truthing" 
to establish the relationship between local 
ecological dynamics of plant and rodent 
populations and satellite sensor readings will be 
needed to determine which of these alternatives 
may apply (19,20). 

Field validation of interpretations, which is 
critical to testing hypotheses derived from 
satellite data, should also be applied to the 
epidemiologic analyses of HPS risk. Our 
approach associates three TM bands and 
elevation with human risk. Interpretation of 
what these bands detect in the environment 
varies (soil moisture, soil type, and vegetation 
structure) (19). Our classification is being used to 
identify other sites with similar reflectance 
patterns in the same bands and characterize the 
structure and dynamics of rodent reservoir 
populations. Preliminary analyses show a good 
relationship between HPS risk predicted from 
satellite imagery and P. maniculatus population 
abundance (r =0.92). 

The case-control model using high-resolution 
spatial data from satellite imagery supports 
previous epidemiologic observations indicating 
that changes in rodent population densities and 
HPS risk could occur dramatically over relatively 
short distances (6). Although extensive areas of 
high and low risk were evident (Figure 3), 
substantial interdigitation of these zones at 
higher resolution created a mosaic of high- and 
low-risk areas. This suggests possibly wide- 
spread "environmental islands" of suitable 
reservoir habitat imbedded within less suitable 
habitat and may account for the apparently focal 
nature of HPS outbreaks and the near random 
distribution of cases relative to their nearest 
neighbors observed in this study (Table 1, k =1). 
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The results also support epidemiologic I nvestiga- 
tions indicating that the only measurable risk 
factor around HPS sites during the 1993 epidemic 
was the abundance of P. maniculatus (6). 

Satellite imagery analysis provides an 
efficient survey of large geographic regions for 
environmental indicators of disease risk affect- 
ing human populations and has the potential to 
make surveillance of disease risk for rare 
zoonotic and vectorborne diseases practical for 
public health applications (18,20-26). For many 
diseases, the basis for the supposition that 
remotely sensed data will be useful for 
anticipating disease risk is that pathogen 
transmission is facilitated by arthropods, whose 
survival and reproduction are influenced by 
variations in temperature and humidity. The 
effect of climatic variability, however, on directly 
transmissible zoonotic agents maintained in 
vertebrate, especially mammalian reservoirs, is 
less certain and has received little attention. 

Additionally, although the reason to assume 
a relationship between climate variability and 
infectious disease outbreaks is clear (27), few 
studies have evaluated whether this presumed 
relationship actually exists. This study indicates 
that if these relationships do occur, they are 
modulated by a number of poorly understood 
ecologic and social conditions that will require 
substantial detailed studies of the pathways 
influencing disease risk. 
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Remote Sensing and Geographic 
information Systems: Charting 
Sin Nombre Virus infections in Deer Mice 
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Robert S. DeBaca.t Edward A. Kuhn,* Pascal Villard,* 

Peter F. Brussard,* and Stephen C. St. Jeor* 
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We tested environmental data from remote sensing and geographic information 
system maps as indicators of Sin Nombre virus (SNV) infections in deer mouse 
{Peromyscus maniculatus) populations in the Walker River Basin, Nevada and 
California. We determined by serologic testing the presence of SNV infections in deer 
mice from 1 44 field sites. We used remote sensing and geographic information systems 
data to characterize the vegetation type and density, elevation, slope, and hydrologic 
features of each site. The data retroactively predicted infection status of deer mice with 
up to 80% accuracy. If models of SNV temporal dynamics can be integrated with 
baseline spatial models, human risk for infection may be assessed with reasonable 
accuracy. 


Remote sensing (RS) and geographic infor- 
mation systems (CIS) are map-based tools that 
can be used to study the distribution, dynamics, 
and environmental correlates of diseases (1,2). 
RS is gathering digital images of the earth's 
surface from airborne or satellite platforms and 
transforming them into maps. GIS is a data 
management system that organizes and displays 
digital map data from RS or other sources and 
facilitates the analysis of relationships between 
mapped features. Statistical relationships often 
exist between mapped features and diseases in 
natural host or human populations (1). Examples 
include malaria in southern Mexico and in 
Asia (3,4), Rift Valley fever in Kenya (5), Lyme 
disease in Illinois (6), African trypanosomiasis 
(7), and schistosomiasis in both humans (8) and 
livestock in the southeastern United States (9). 
RS and GIS may also permit assessment of 
human risk from pathogens such as Sin Nombre 
virus (SNV; family Bunyaviridae), the agent 
primarily associated with hantavirus pulmonary 
syndrome (HPS) in North America (10,11). RS 
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USA; fax: 775-784-1620; e-mail: boonecgscsr. nevada.edu. 


and GIS are most useful if disease dynamics and 
distributions are clearly related to mapped 
environmental variables. For example, if a 
disease is associated with certain vegetation 
types or physical characteristics (elevation, 
average precipitation),RSandGIS cou I d i dent i f y 
regions where risk is relatively high. 

We examined whether RS and Gl S data were 
useful indicators of the spatial pattern of SNV 
infections in populations of the primary rodent 
host, the deer mouse (Peromyscus maniculatus) 
(12-15). Our approach involved determining the 
infection status of rodents at 144 field sites, 
collecting RS and Gl S data for each site, testing 
for statistical relationships between these data 
and infection, using the statistical relationships 
to retroactively classify infection status of 
rodents at these sites, and using the classifica- 
tions to estimate prediction accuracy. Predictions 
derived from RS and GIS data could identify the 
ecologic settings where human exposure to SNV 
is most likely to occur. 

SNV and its Host 

Since thefirst recognized outbreak of H PS in 
the southwestern United States in 1993, 
approximately 240 cases have occurred, with a 
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death rate of approximately 40% (J . Mills, pers. 
comm.) (16). Information about SNV host-virus- 
environment relationships is limited (16,17). No 
simple relationships have been found between 
host density and antibody seroprevalence (16-18), 
but morecomplex nonlinear relationships appear 
toexist (17). SNV infectionsalsoappearto beless 
frequent in relatively high- or low-elevation 
habitats (16,17). 

Study Design 

Types of Data 

RS data are commonly used to generate 
maps of vegetation types. Vegetation types can 
be useful indicators of environmental charac- 
teristics, including moisture, soil type, and 
elevation. However, transforming RS images 
into vegetation maps can be subjective and 
imprecise (19,20); therefore, we supplemented 
our vegetation maps with other RS/GIS data, 
including elevation, slope, vegetation density, 
and hydrology. 

We sampled rodents over four field seasons 
(j une to October during 1995 to 1998). However, 
in 1997, population densities of deer mice in our 
study area averaged approxi mately 25% of 1995- 
96 and 1998 levels (unpub. data, Boone et al.). 
Most of the 47 sites sampled in 1997 had three or 
fewer deer mice, and 14 sites had none. Simplet- 
tests (SAS ver. 6.10) showed that the mean 
number of animals per site was statistically 
equivalent in 1995, 1996, and 1998 (11.1, 10.3, 
and 8.9 animals per site, respectively; p>0.10for 
all comparisons), but differed significantly in 
1997 from all other years (2.4 animals per site; 
p <0.0001 for all comparisons). I n 1997, antibody- 
positive animals were significantly less likely 
to be positive by reverse transcription- 
polymerase chain reaction (RT-PCR) for viral 
RNA in the blood than in any other year, 
suggesting unusual infection patterns (25% 
were PCR positive in 1997 and ?60% to 70% in 
other years; chi-square test, p <0.002 for all 
comparisons of 1997 to other years; p >0.20 for all 
comparisons of years excluding 1997). On the 
basis of these tests, we pooled data from 1995, 
1996, and 1998 and excluded 1997 data from all 
analyses because host density and infection 
dynamics appeared atypical and likely to obscure 
the baseline spatial infection patterns we sought 
to identify (21). 


Infection Status of Sites 

Presence of SNV infections is commonly 
i nferred by determi ni ng anti body seroprevalence 
in a host population (14,16-18,21-23). However, 
antibody prevalence at the same site may vary 
considerably (<5%to >60%) over relatively brief 
periodsof<l year (17,18,22,23), probably because 
of rapid turnover of rodent populations through 
death, reproduction, dispersal, and migration. 
We focused on the presence or absence of SNV 
infections inferred from antibody data, a more 
stable measure than antibody prevalence. 
However, determining infection status is 
complicated by several factors: animals may 
remain antibody-positive well after the trans- 
missible phase of an infection (17); noninfectious 
but anti body-positive deer mice may migrateto 
a site where no active SNV infection is present; 
and detectable antibody response requires at 
least 1 to 2 weeks to develop in newly infected 
animals (17). 

Because of these uncertainties, we used two 
criteria to demonstrate the effect of cl assifi cation 
on analytical outcome. "Status 1" classified sites 
with one or more antibody-positive animals as 
positive (active infection present). This criterion 
may havefalsely assigned positive status to some 
sites where no active, transmissible infections 
were present. "Status 2" required two or more 
antibody-positive deer mice or an overall 
antibody seroprevalence of at least 10% for a site 
to be classified positive. This criterion may have 
falsely assigned negative (active infection 
absent) status to some sites that had a single 
infectious animal. 

Site Selection 

Our study area was the Walker River Basin, 
a 10,200-km2 region in western Nevada and east- 
central California northeast of Yosemite Na- 
tional Park (Figure 1). At least nine cases of H PS 
have occurred in the area since 1993. Major 
vegetation types in the river basin along an 
increasing elevational gradient (1,200 m to 
3,760 m) are salt desert scrub, sagebrush-grass 
scrub, pihon-juniper woodland, coniferous forest, 
montane shrubland, and alpine tundra, with 
riparian habitat and meadows at a wide range of 
elevations (24). 

We compiled a GIS database for the study 
area, including a second-generation map of 
vegetation types (Figure 1) (25). The vegetation 
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Figure 1. Location of Walker River Basin (17) and its eight major vegetation types, as well as developed areas. 
Pihon-juniper woodland and montane shrubland tend to be highly interspersed and were combined for visual 
clarity. Because meadows occurred in very small patches, they could not be represented on this map. Map 
generated at Utah State U niversity as part of the GAP conservation mapping project. 
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map, which was generated from Landsat 
Thematic Mapper images and digital elevation 
data, had a 100-hectare mapping unit. We 
aggregated the 36 vegetation subtypes on the 
GAP map into the eight general vegetation types 
described above. To estimate vegetation density, 
we used the normalized difference vegetation 
index (NDVI), a transformation of near infrared 
(TM band 4) and red wavelengths (TM band 3) 
correlated with the amount and productivity, or 
rate of plant growth, of vegetation (5,26,27). The 
standard deviation of NDVI within a local area 
was calculated to estimate the uniformity of 
vegetation density at each field site. Elevation 
and slope (i.e., steepness) data were derived from 
the 2-arcsecond digital elevation model of the 
U.S. Geological Survey. Because riparian zones 
could influence rodent population densities and 
facilitate rodent dispersal across arid regions, we 
calculated proximity to streams and bodies of 
water on the U.S. Geological Survey's 1:100,000- 
scaledigital line graph datasets. 

In 1995, we sampled rodents at 42 sites 
before the GAP map became available. These 
sites were selected as representative of the five 
most common vegetation types in the Walker 
River Basin. I n 1996 and 1998, full Gl S datasets 
and the GAP map were used to distribute 102 
new field sampling sites systematically across 
the widest possible range of environmental 
conditions. We categorized each GIS variable 
according to its relevance to each of the eight 
vegetation types. For example, 'distance to 
streams' was a meaningful distinction within salt 
desert but not within riparian habitat; elevation 
varied substantially within sagebrush scrub but 
not within alpine tundra. For each vegetation 
type, the relevant variables were divided into 
high and low ranges. Theresulting binary cl asses 
for each variable were then intersected in Gl S to 
produce distinct environmental "combinations," 
or strata, for each vegetation type (Figure 2). 
Randomly located sample sites were selected 
within each stratum so that they were within 0.5 
km of a passable road and at least 1 km from any 
other sample site (Figures 1,3). The number of 
replicates within each stratum (including 1995 
sites, which were included retroactively) was 
proportional to its spatial extent, with a minimum 
sample size of two. This Gl S-based stratification 
is a more objective and randomized variation of 
thegradsect sampling method (28, 29). 


All samples were collected from early] uneto 
early October to minimize seasonal effects on 
host density and antibody prevalence (17). 
Seasonal influences were mini mi zed by sampling 
the replicate sites within each environmental 
stratum at different times throughout the field 
season. 

Field and Laboratory Procedures 

Deer mice were live-trapped at all field sites 
according to a fixed protocol (17). Each site had 48 
live-traps in pi ace for 3 days. A blood sample was 
collected from each deer mouse by retroorbital 
puncture with a heparinized capillary tube or 
Pasteur pipette. Blood samples were placed on 
dry ice and returned to the laboratory for 
enzyme-linked immunosorbent assay testing for 
immunoglobulin G antibody to SNV, which 
indicates current or past infections (14). Relative 
population density was estimated by counting 
the number of animals captured during a 
trapping session. 

Analytical Methods 

Of the 144 sites sampled in 1995, 1996, and 
1998, 25 were excluded from analysis because no 
deer micewerecaptured. Status 1 classified 38 of 
the remaining 119 sites as negative and 81 as 
positive. Status 2 classified 70 sites as negative 
and 49 as positive (i.e., 32 sites had differing 
infection status under the two criteria). We 
tested (by chi-square, SAS ver. 6.10, PROG 
FREQ) for differences among the proportion of 
positive sites for each vegetation type. Then, with 
a canonical linear discriminant function analysis 
[DFA] [SAS ver. 6.10, PROC Dl SCRIM], we 
examined relationships between infection status 
and the alternate set of RS and GIS variables 
with slope, elevation, density and uniformity of 
vegetation, and distance from streams as 
indicators of SNV infection status (3,17). Prior 
probabilities were adjusted to reflect actual 
proportions of positive and negative sites. 

The relationships derived from these two 
analyses were then used to classify sites 
retroactively according to their expected infec- 
tion status. Because error rates were not 
distributed evenly among sites classified as 
positive and negative, we present these results 
separately. Classification accuracy is a general 
estimate of the prediction accuracy of each 
method if it were applied to new sites in a similar 
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Figure 2. Environmental strata within the mapped extent of the sagebrush-grass scrub vegetation type. Each 
color represents a unique combination of high or low vegetation density index, standard deviation of vegetation 
density index, slope, elevation, and distance from stream. Whiteareas represent the other seven vegetation types 
(strata not shown). 
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Figure 3. Seroprevalence of IgG antibody to Sin Nombre virus at sample sites in the Walker River Basin. 
(According to Status 1, sites with 0% prevalence were negative and all other sites were positive.) 
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landscape. For each prediction accuracy esti- 
mate, we calculated 95% confidence limits 
derived from the binomial distribution. Standard 
deviation for the binomial distribution is 


n 

where p = accuracy estimate and n = number of 
samples. A normal approximation of confidence 
limits was obtained by multiplying the standard 
deviation of each estimate by the t-table value 
associated with 95% confidence and the 
appropriate number of samples. These confi- 
dence intervals also allowed us to determine 
whether classification accuracy differed signifi- 
cantly between methods. 

Results 

Vegetation Types 

The proportion of positive sites in salt desert 
scrub (34% of 29 sites by Status 1, 14% by Status 
2) was significantly lower than in any other 
vegetation type(p =0.05 criteria for significance). 
Nosignificant differences were found among any 
of the other seven vegetation types, where 
positive sites were more common by both Status 1 
(50% to 100%) and Status 2 (50% to 83%) 
(Figure 3). By assigning the predominant 


infection status to all sites within a given 
vegetation type, overall classification accuracies 
of 76% (Status 1) and 59% (Status 2) could be 
achieved (Table 1). The Status 1 criterion 
resulted in better classification accuracy (for 
negative sites and for all sites combined) than 
Status 2. F or both Status 1 and Status 2, positive 
classification was more accurate (> 88%) than 
negative classification ( 50%). 

DFAsfor both Status 1 and Status 2 produced 
significant canonical correlations showing that 
negative sites were associated with low eleva- 
tions and sparse vegetation (Table 2). These 
qualities most often occur in salt desert scrub 
(24). In contrast, positive sites were higher and 
generally had more dense but less uniform 
vegetation. Slope and distance from streams 
were relatively unimportant factors. For both 
Status 1 and Status 2, negativeclassifi cation was 
more accurate than positive classification 
(Table 1). Positive classification was more 
accurate i n Status 2 than i n Status 1. 

Discussion 

RS and G I S data were usefu I i ndi cators of the 
SNV infection status of deer mice in our study 
area. Sites with typical salt desert scrub 
characteristics were less likely to have infected 
mice than other sites. If the 25 sites where no 
deer mice were captured (primarily salt desert 


Table 1 . Classification and prediction accuracies for Sin Nombre virus infection status, by vegetation type method and 
canonical discriminant function analysis 


Site in Walker River Basin 

% classified negative 

% classified positive 

% overall accuracy 

Vegetation type Status 1 

50 ± 123 

88 ± 8 

76 ±10 

Vegetation type Status 2 

36 ± 9 

92 ± 6 

59 ±12 

DFA^ Status 1 

84 ± 9 

69 ±11 

74 ±10 

DFA Status 2 

84 ± 9 

80 ±10 

82 ± 9 


®The error terms following each estimate are 95% confidence intervals derived from the binomial distribution. 
‘’DFA =discriminant function analysis. 


Table 2. Sin Nombre virus infection status of sites in the Walker River Basin: canonical discriminant function analyses 
for Status 1 and Status 2 


Canonical Canonical loadings^’ 



correlation^ (p-value) 

Elevation 

NDVI^ 

NDVI Std'' 

Slope 

Streams® 

Status 1 

0.41 (0.0003) 

0.77 

0.52 

0.69 

0.06 

0.04 

Status 2 

0.53 (0.0001) 

0.79 

0.46 

0.68 

0.09 

0.05 


®The canonical correlation and its significance level define the overall association between infection status and the indicator 
variables. 

‘^Canonical loadings for each indicator variable indicatetheir relative importance in producing a significant overall canonical 
correlation. 

"^NDVI =vegetation density. 

‘^NDVI Std. = uniformity of vegetation density. 

^Streams = distance to nearest mapped stream or body of water. 
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scrub sites) had been incorporated into our 
analyses as negative sites, this relationship 
would have been more pronounced. The 
relationship may be explained by the level of 
connectivity (i.e., biological interchange) among 
host popu I at i on s. S a 1 1 desert scruborsimilararid 
habitats in the western United States are 
frequently dominated by heteromyid rodents 
(kangaroo rats, pocket mice) rather than by deer 
miceand other potential hosts for SNV. Although 
deer mice were found in salt desert scrub in the 
Walker River Basin and were sometimes locally 
abundant, their overall population density was 
somewhat lower than in other vegetation types, 
and they were more likely to be locally absent (17). 
Wesuspect that SNV infections are less likely in 
deer mouse populations that inhabit such regions 
because of their relative isolation from neighbor- 
ing populations (30,31). Such fragmentation of 
host populations may reduce the rate of disease 
propagation across space and the frequency of 
infection recurrences within local sites. This 
hypothesis is supported by the clustering of 
negative sites in landscapes dominated by salt 
desert scrub (F igure 3), despitethefact that some 
of these sites had relatively dense deer mouse 
populations. 

Spatial Versus Temporal Disease Patterns 

Because the RS and GIS maps summarize 
relatively fixed spatial properties of the 
environment, we focused on investigating the 
corresponding spatial patternsof SNV infections. 
SNV infections also exhibit temporal dynamics 
(13,16-18,22,23) superimposed on the baseline 
spatial pattern. However, a robust temporal 
study would require many years of replicated, 
longitudinal field data, as well as real-time RS 
data describing temporally variable environ- 
mental characteristics (such as climatic 
variables) for the corresponding period. Wedid 
not incorporate weather or climate data into 
GIS because weather monitoring stations are 
widely scattered throughout most of the study 
area, preventing meaningful extrapolations to 
most of the field sites. 

Sampling Design 

Because characterizing large-scale spatial 
disease patterns requires a large sample size, we 
maximized the number of sites sampled rather 
than visiting fewer sites on multiple occasions. 
This cross-sectional approach captured substan- 


tial ecologic diversity and provided statistical 
replicates of sites with similar characteristics. 
The disadvantage of the approach was a degree of 
uncertainty in determining the actual infection 
status at each site. However, when generaliza- 
tion of results is an important goal, a large, 
replicated, and diverse dataset that has a modest 
degree of measurement error is statistically 
preferable to a smaller, more precisely measured 
but poorly replicated dataset (32). 

Comparison of Methods (Table 1) 

The vegetation type approach was based on 
possible relationships between infection status 
and a preexisting vegetation classification that 
might or might not be relevant to deer mice and 
SNV infections. DFA, in contrast, generated a 
linear function that best distinguished the 
properties of positive and negative sites. Our 
results suggest that DFA yields a better balance 
between classification accuracies for positive and 
negative status (especially for Status 2). 

Thevegetation type method could not classify 
negative status as effectively as the DFA, and 
balance between error rates for positive and 
negative classifications was poor. This could be a 
result of using predefined vegetation types 
(rather than making environmental distinctions 
from actual infection patterns) or inaccuracies in 
identifying and mapping vegetation types. Site 
visits suggested that the DFA identified sites 
with pronounced salt desert features more 
effectively than the vegetation map. The 
substantial environmental variability within the 
mapped extent of salt desert scrub was easily 
captured by the set of RS and GIS variables but 
was analytically "invisible" to our aggregated 
GAP map. Some variability might have been 
captured bytheGAP map's 36 original vegetation 
subclasses, but using all these in our analysis 
would have presented serious statistical problems. 

Other analytical approaches are possiblethat 
were not presented here. For example, decision 
tree analysis (33,34) offers advantages if nonlinear 
relationships exist; hierarchical information on 
the effects of each predictor variableisdesired; or 
ease of interpretation is important (29). 

Classification and Prediction Accuracy 

Classification accuracy varied significantly 
between the Status 1 and Status 2 criteria 
(Table 1), with Status 2 giving better classifica- 
tion balance for DFA and Status 1 producing 
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better results for the vegetation type analysis. 
Unfortunately, the biological significance of 
these analytical differences is difficult to 
determine. However, the infection status of 73% 
of the sites was classified similarly by the two 
criteria. The remaining 32 sites of ambiguous 
infection status might represent regions where 
infection status changes with relatively high 
frequency. If so, this produces an intrinsic 
limitation in the capabilities of the methods we 
present. The choice of technique might be based 
on the relative risks and costs of false- negative 
versus false-positive predictions. 

Both methods may have occasionally been 
unable to detect positive sites because of failure 
to capture positive deer mice. The likelihood of 
this error would be proportional tothe number of 
resident animals not captured at a site. Our 
longitudinal data (17; unpub. data, Boone et al.) 
suggest that the 3-day sampling sessions 
captured most of the animals present at a site. 
Duringfour 7-day trapping sessions, 86%±9%of 
animals were captured during the first 3 days of 
trapping. Additionally, examination of 123 3-day 
sessions within the context of their extended 
longitudinal infection timelines (17; unpub. 
data, Boone et al.) suggested that infection 
status was classified with 85% (Status 1) and 
81% (Status 2) accuracy. 

A cross-sectional, replicated, and randomized 
sampling approach should capture most sites 
while they exhibit their most typical infection 
status. H owever, a 'background 'rate of classifica- 
tion errors is to be expected regardless of 
analytical method, given thetemporally dynamic 
nature of SNV infections (17,18,22,23). For 
instance, even where infection status is 
predominantly positive, some sites may be 
sampled during atypical periods when infection 
is temporarily absent; the reverse could also 
occur. Additionally, a subset of sites might 
frequently change their infection status and not 
exhibit primary infection status.Thusit might be 
difficult to improve upon the highest overall 
classification success we achieved (with DFA and 
Status 2 criterion), unless temporal infection 
dynamics are incorporated into the predictive 
model. Another option would be to omit sites from 
analysis if they fail to meet unambiguouscriteria 
for positive or negative status; however, this 
might result in the loss of biological insight. 


Future Directions 

Weexplored the ability of RS and Gl S data to 
predict the baseline spatial patterns of SNV 
infections across an ecologically variable land- 
scape. Our findings should be at least somewhat 
relevant to a number of other regions in the arid 
western United States, especially if infection 
dynamics are ultimately driven by host connec- 
tivity patterns. To expand these findings, we 
developed methods to filter environmental data 
to remove statistical noise and a computer 
simulation model to explore infection dynamics 
on a variety of virtual landscapes. Further work 
will focus on the role of landscape structure in 
producing spatial patterns of disease (35). For 
instance, deer micein small patches of salt desert 
scrub within a matrix of more desirable habitat 
types might be more likely to be infected than 
mice living in large contiguous regions of salt 
desert scrub. Finally, it would be useful to test 
other types of RS and GIS data as possible 
indicators of SNV infections. 

Further work is needed to identify possible 
climatic correlates of periodic outbreaks and the 
degree to which useful indicators of these 
outbreaks can be derived from RS and GIS data 
sources. I n contrast to predictions in large-scale 
outbreaks, specific a priori predictions of 
temporal SNV infection dynamics in local sites 
may remain difficult. Once infections are 
initiated at a site (presumably by random 
dispersal events), changes in antibody and virus 
prevalence cannot be easily explained by changes 
in host density or environmental factors (17,18). 
However, it should be possible to estimate the 
frequency (if not the specific timing) of new 
infections as a function of a site's local 
environment. Additionally, extended longitudi- 
nal studies could identify typical infection 
trajectories of sites based on their environmental 
characteristics or demographic profiles of their 
host populations. When combined, these ap- 
proaches should advance our ability to quantify 
and predict disease dynamics and human risk. 
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Two human deaths caused by Australian bat lyssavirus (ABL) infection have been 
reported since 1996. Information was obtained from 205 persons (mostly adults from 
south Brisbane and the South Coast of Queensland), who reported potential ABL 
exposure to the Brisbane Southside Public Health Unit from November 1,1996, to 
January 31, 1999. Volunteer animal handlers accounted for 39% of potential exposures, 
their family members for 12%, professional animal handlers for 14%, community 
members who intentionally handled bats for 31%, and community members with 
contacts initiated by bats for 4%. The prevalence of Lyssavirus detected by fluorescent 
antibody test in 366 sick, injured, or orphaned bats from the area was 6%. Sequelae of 
exposure, including the requirement for expensive postexposure prophylaxis, may be 
reduced by educating bat handlers and the public of the risks involved in handling 
Australian bats. 


Australian bat lyssavirus (ABL) was first 
reported in J uly 1996 in a black flying fox 
(Pteropus alecto) from Ballina, New South Wales 
(1,2). ABL has been confirmed in five species of 
Australian bat: four species of flying fox 
(suborder Megachiroptera, genus Pteropus) and 
one species of insectivorous bat (suborder 
Microchi roptera, Saccolaimus flaviventris). Two 
cases of human ABL infection have been 
reported. The first case occurred in a 39-year-old 
female animal handler from Rockhampton, 
Queensland, in November 1996, within 5 weeks 
of her being scratched and possibly bitten by a 
yellow-bellied sheath-tailed bat (S. flaviventris) 
(R. Taylor, pers. comm.). The second case 
occurred in a 27-year-old woman from Mackay, 
Queensland, in December 1998, >2 years after a 
bite from a flying fox. Both patients died (3,4). 


Address for correspondence: Bradley J . McCall, Brisbane 
Souths! dePublic Health Unit, P.O. Box 333, 39 Kessels Road, 
Cooper's Plains, Old 4108, Australia; fax: 61-7-3000-9130; e- 
mail: mccalibcghealth. qld.gov.au. 


ABL is a member of the family Rhabdoviri- 
dae. Although ABL possesses marked serotypic, 
antigenic, and molecular sequence similarities to 
classic rabies virus, it represents a distinct, new 
genotype, genotype 7 of the Lyssavi rus genus (5). 
The clinical signs of ABL infection in the two 
human cases were consi stent with those of cl assi c 
rabies infection and included a diffuse, nonsuppu- 
rativeencephalitisthatledtodeath (3,4). Bats with 
ABL infection are frequently reported to have 
had hind limb paresis. While most infected bats 
are depressed when found, some exhibit 
uncharacteristic aggression toward humans or 
other bats. Frequently, a nonspecific, nonsuppu- 
rative meningoencephalitis is seen in brains of 
infected animals (6,7). Vaccine protection trials 
in mice conducted at the Centers for Disease 
Control and Prevention (CDC), Atlanta, Georgia, 
supported the decision to use human diploid cell 
vaccine (HDCV) for human ABL prophylaxis 
(7-9). Historically, Australia has been considered 
freeof rabies and rabieslike viruses. Thus, before 
thefirst human case of ABL infection in 1996, no 
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measures existed to prevent rabies or rabieslike 
disease acquired as a result of contact with 
Australian domestic animals or wildlife. Since 
the first human ABL case, the Queensland 
Health Department, in accordance with the 
recommendations of the national Lyssavirus 
Expert Group, has provided postexposure 
prophylaxis (PEP) to persons who report 
potential exposure to ABL through bites, 
scratches, and permucosal or percutaneous 
exposure to bat saliva or neural tissue (9,10). 
Preexposure prophylaxis is recommended for 
persons who report frequent contact with bats. 

Colonies of flying foxes are common in 
suburban areas of southeast Queensland. The 
black flying fox (Pteropus alecto) and the grey- 
headed flying fox (P. poliocephalus) live there 
throughout the year, and the little red flying fox 
(P. scapulatus) occurs seasonally. While the 
population of flying foxes may be decreasing in 
southeast Queensland, fragmentation of colonies 
has resulted in a wider distribution of smaller 
colonies (L. Hall, pers. comm.). Direct contact 
with bats by the general public and animal 
handlers is not uncommon (11). Volunteer animal 
handlers rehabilitatesick, injured, and orphaned 
bats and are frequently bitten, scratched, or 
exposed to bat saliva. Since November 1996, the 
BrisbaneSouthside Public Health Unit (BSPH U) 
and other state public health units have been 
involved in coordinating lyssavirus PEP. This 
article describes the pattern of potential human 
exposure to ABL reported to the Communicable 
Disease Control Section of BSPHU between 
November 1996 and J anuary 1999 and subse- 
quent PEP. Disease prevalence findings are 
presented for bats surveyed in southeast 
Queensland by the Animal and Plant Health 
Service of the Queensland Department of 
Primary I ndustries. 

The Study 

During the study (November 1, 1996, to 
J anuary 31, 1999), the Communicable Disease 
Control Section of BSPHU served a population of 
approximately 1.1 million persons in several local 
government areas, south Brisbane (part of the 
Brisbane City Council Area), Logan, Redlands, 
Beaudesert, and Gold Coast (12). Al I persons who 
reported a potential ABL exposure (bat bite, 
scratch, percutaneous, or permucosal exposureto 
bat saliva or neural tissue) were asked to 
completea standard questionnaire, which sought 


demographic information (including occupation, 
history of professional or volunteer bat handling, 
history of rabies vaccination, potential rabies 
exposure [bite, scratch, provoked, unprovoked], 
circumstances that led to the exposure, treat- 
ment received, and any laboratory investigation 
of the bat). 

A separate questionnaire was completed for 
each occasion a person contacted BSPHU to 
report potential ABL exposure. Potential expo- 
sures were reported retrospectively, and the 
dates of notification and potential exposure for 
each case were included. All information was 
recorded and analyzed by usingan Epi-I nfo 6.04b 
database (13). Age and gender-specific notifica- 
tion rates were calculated by using estimated 
resident population data for 1997 (12). 

During the same period, healthy bats, sick 
and injured bats, and bats involved in a potential 
human exposureto ABL weretestedfor infection 
with a fluorescein-labeled antirabies monoclonal 
globulin (CENTQCQR) in a direct fluorescent 
antibody test (DFAT) on fresh brain impression 
smears at the Queensland Department of 
Primary I ndustries Animal Research I nstituteor 
at the CSIRQ Australian Animal Health 
Laboratories. Material from most bats that 
tested positive for ABL infection and from bats 
associated with a potential human exposureto 
ABL were sent to either the Australian Animal 
Health Laboratory or Queensland Health 
Scientific Services for confirmation by DFAT, 
virus isolation, and polymerase chain reaction. 

Results 

Atotal of 205 notifications toBSPHU met the 
criteria for potential ABL exposure during the 
study period, an average annual notification rate 
of 8.1/100,000. Complete information was 
obtained from 202 persons. Total notifications 
included 86 males and 119 females (M :F ratio of 
1:1.38). The age- and gender-specific average 
annual notification rates are presented in Figure 
1. Most reported potential exposures (116 of 204) 
were among persons 19 to 49 years of age. The 
months of potential ABL exposure and notifica- 
tion arepresented in Figure2. Most notifications 
(131of 205) were made with in 2 months of each of 
the two fatal human cases. Nine (11%) of 80 
notifications made in the 2 months following the 
first reported case were related to exposures that 
occurred >2 months before the fi rst human case 
was publicized. The median interval between 


Emerging Infectious Diseases 


260 


Vol. 6, No. 3, May-]une 2000 



Research 


Mo(tftc«tton«i1 00 jOOO 



Figure 1. Age and gender-specific average annual 
notification rates of potential human exposure to 
Australian bat lyssavirus (n =204) south Brisbane 
and South Coast, Queensland, 1996-1999. 



Figure 2. Dates of potential Australian bat lyssavirus 
exposures and notifications to the Brisbane Souths! de 
Public Flealth Unit, south Brisbane and South Coast, 
Queensland, 1996-1999. 


exposure and notification of these 80 potential 
exposures was 17 days (0 to 1,080 days). I n the 2 
months fol I owi ng the second case, 22 (43%) of 51 
notifications were related to potential exposures 
that had occurred before the reporting of the first 
human case. The median interval between 
exposure and notification of the 51 potential 
exposures was 728 days (0 to 2,907 days). A 
further 14 (27%) of the 51 notifications were 
related to potential exposures that had occurred 
since the first human case but had not been 
reported to BSPH U at the time of exposure. 

Season of Exposures 

Potential exposures to ABL were reported to 
have occurred from 1991 to 1999, most during 
spri ng and summer (September to F ebruary) (n = 
151, 74%). Whiletheoccurrenceofthetwohuman 
ABL cases in spring and summer influenced the 
reporting of potential exposures at these times, 
this trend of increased spring-summer potential 
exposures persisted in the period between the 
two reported human ABL cases. The highest 
number of potential exposures (105) was reported 
in theyear of thefirst human case; 99 occurred in 
the spring or summer of 1996-97. 

Groups at Risk 

Notifications were categorized according to 
the person's life-style and occupational potential 
for exposuretoABL (Table). Thegroup at highest 
risk, volunteer bat handlers, reported 79 (39%) 
potential exposures; 8 of these handlers reported 
a second potential exposure during the study. 


Table. Groups at risk for exposure to Australian bat lyssavirus 


Groups at risk 

No. of 
potent! al 
exposures 
(n=203) 

M ean age 
and age 
range (yrs) 

Gender 

(m/f) 

Median 
interval 
between 
exposure and 
notification (d) 

Bite/non bite 
injury 
(n=202) 

P revoked 
(%) 

(n=202) 

Volunteer bat handlers 

79 

40.5 

15/64 

19 

56/23 

79/79 



(16-83) 


(0-2,105) 


(100) 

Flousehold or family member 

24 

17.5 

12/12 

27 

18/6 

24/24 

of volunteer bat handlers 


( 5-51) 


(0-1,809) 


(100) 

Professional animal worker 

28 

34 

15/13 

4 

13/14 

27/27 



(17-69) 


(0-1,818) 


(100) 

Community-intentional 

63 

39 

40/23 

10 

41/22 

62/63 

potential exposure 


( 6-85) 


(0-2,907) 


( 98) 

Community-uni ntential 

9 

32 

4/5 

2 

3/6 

4/9 

exposure 


(16-49) 


(0-32) 


( 44) 
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Twenty-four (12%) notifications of potential 
exposure were among household or family 
members of volunteer bat handlers. Professional 
animal handlers (e.g., veterinarians, wildlife 
biologists, park rangers) reported 28 (14%) 
exposures. Community members who handled 
bats (usually in thecourse of freeing them from a 
fence or entanglement) reported 63 (31%) 
potential exposures. Community members re- 
ported 9 unintentional potential exposures in 
which contact was initiated by the bat. 

The pattern of notification varied within 
groups at risk during the study. The number of 
potential exposures reported by volunteer and 
professional animal handlers declined. Notifica- 
tions by all groups were highest in the months 
after the reported fatal human cases of ABL 
infection. The number of potential exposures 
reported in the 2 months after the first human 
case (n = 80) was higher than after the second 
case (n =51), particularly among volunteer bat 
handlers, who reported the highest number of 
potential ABL exposures in the 2 months after 
the first human case (43 [53%] of 80), decreasing 
to 12 (24%) of 51 i n the 2 months after the second 
case. Notifications of potential exposures among 
community members who had intentionally 
handled bats rose from 11 (14%) of 80 in the 2 
months after the first human ABL case to 23 
(45%) of 51 in the 2 months after the second 
human ABL case. 

Nature of Exposure 

Potential exposures were classified as biteor 
nonbite exposures in accordance with interna- 
tional recommendations (14). Most potential 
exposures were bites (n = 132, 64%). The ratio of 
biteto non bite potential exposures within groups 
at risk was highest among volunteer bat handlers 
(56:23) (Table). Potential exposures associated 
with unintentional contact with bats by 
community members were predominantly 
scratches (3 bite: 6 non bite), whereas potential 
exposures from intentional contact with bats by 
all other risk groups were predominantly by 
bites (128 bites: 65 nonbites). Potential 
exposures were categorized as provoked 
(arising from intentional contact with a bat) or 
unprovoked (a contact initiated by the bat). 
Most potential exposures (97%) were described 
as provoked (Table). 


Treatment 

PEP was offered to all persons who reported 
potential ABL exposures, in accordance with 
international and Australian recommendations 
(8,14). Standard PEP for unvaccinated persons 
consisted of human rabies immune globulin 
(H Rl G, 20 I U/kg) on day 0 and 5 doses of H DCV 
administered on days 0, 3, 7, 14, and 28. PEP for 
immunized persons consisted of 2 booster doses 
of HDCV administered on days 0 and 3. A 
national shortage of HRIG required modifica- 
tions to the standard PEP regimen. Sixty-two 
potentially exposed persons received standard 
PEP, 100 received 5 doses of HDCV only, 16 
vaccinated persons received 2 booster doses of 
H DCV, and 25 persons did not receive treatment 
when the bat tested negative. Two persons 
refused vaccination because of concerns about 
potential vaccine side effects. Sixteen persons 
ceased treatment when the bat tested negative. 

The estimated cost of providing PEP during 
this study was A$137,368, which included 
A$30,930for medical services funded through the 
Commonwealth Medicare system (calculated on 
the cost of six visits to a medical practitioner for 
each person requiringa 5-dose course of PEP and 
three visits for each person requiring 2 doses of 
PEP); A$8,200 for public health officers who 
interviewed potentially exposed persons; 
A$10,600 for laboratory testing of the bats; and 
A$87,638 for HDCV and HRIG. The cost of all 
vaccines was met by the Queensland Health 
Department. 

ABL Test Results in Bats 

All bats retrieved from a human exposure 
incident underwent postmortem examination 
and testing for ABL infection. Thirty-six bats 
were tested; two were positive on DFAT and 
polymerase chain reaction testing for 
Lyssavirus— a black flying fox and a little red 
flying fox. The tested bats included 13 black 
flying foxes, 11 grey-headed flying foxes, 5 little 
r^ flying foxes, and 7 insectivorous bats. 

In a separate investigation, the Queensland 
Animal Research I nstitute tested bats by DFAT 
on brain impression smears for evidence of ABL 
infection since J une 1996. From November 1, 
1996, to J anuary 31, 1999, some 153 healthy 
wild-caught flyingfoxes; 181 healthy wild-caught 
insectivorous bats; 366 sick, injured, or orphaned 
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flying foxes; and 39 sick or injured insectivorous 
bats from the area served by BSPH U and greater 
Brisbaneweretested. Of these, 21 (6%)of the366 
sick, injured, or orphaned flying foxes tested 
positive for ABL infection, including the 2 
involvedin human exposures in theBSPHU area. 
All other bats tested negative. 

Discussion 

This is the first description of PEP provided 
to an Austral i an community after the recognition 
of human risk for ABL infection. The first human 
case triggered a large national public health 
campaign and considerable public awareness 
about the risks from bats, particularly in 
communities such as south Brisbane and the 
South Coast of Queensland, where large colonies 
of bats live close to human urban populations and 
bat/human interaction is not uncommon. In- 
creased concern was demonstrated by the large 
number of notifications of potenti al exposure that 
followed reports of the two human cases 
(Figure 2). Most potential human exposures 
were among adults (ages 25 to 49). The i ncr eased 
proportion of women reflects the high proportion 
of female volunteer bat handlers in the study 
population. 

The first 2 months of notifications represents 
a catch-up period in which PEP was provided to 
persons with exposures dating back several 
years. Relatively few notifications occurred after 
the first 2 months of the study, and it was only 
after the second human case, which had an 
assumed incubation period of approximately 2 
years, that another cluster of notifications 
occurred (4). The median interval between 
potential exposure and notification increased 
from 17 days for those notified in the 2 months 
after the first human ABL case to 728 days for 
those notified in the 2 months after the second 
case. The potential exposures reported in the 2 
months after the second case i ncl uded 22 persons 
who were potenti al I y exposed beforethefi rst case 
and 14 with >1 month between potential 
exposure and notification. The second human 
case with its prolonged incubation period 
rei nforced the pubi ic perception of the severity of 
this disease and prompted more notifications. 

Potential exposures occurred most commonly 
in spring and summer, coinciding with the 
birthing season (October to December) of the black 


and grey-headed flying foxes in southeastern 
Queensland (15). During each birthing season in 
southeastern Queensland, 100 to 300 neonatal 
and j u veni I e bl ack or grey-headed f I yi ng foxes are 
reared by volunteer bat handlers (H. Luckhoff, 
pers. comm.). These orphans are commonly 
assumed to have been abandoned or separated 
from their dams. Frequently, orphans are found 
still clinging to the body of their dam. Further 
research is required to identify any association 
between orphaned bats and theABL statusof the 
dam. A case of clinical disease in an in-care 
juvenile black flying fox and the associated 
exposure of eight humans has been described (6). 

Most potential exposures (107 [52%] of 205) 
were reported from groups who handled bats. 
These groups were the target of initial public 
health information campaignsto raise awareness 
of the risks for ABL infection. PEP was provided 
to members of these groups after the first human 
case, and a recommendation was issued that all 
workers in these fields be vaccinated with FI DCV 
and that unvaccinated persons, including family 
members of volunteer bat handlers, not handle 
bats. Seventy-two (35%) of 205 potential 
exposures occurred among members of the 
community. Most of these (63 [88%] of 72) had 
rescued a trapped or fallen bat. The test results 
from bats indicate that sick, injured, or orphaned 
bats have a significantly higher crude prevalence 
of ABL infection (p <0.001) than healthy wild- 
caught or captive bats. Consequently, the risk for 
ABL exposure among volunteer and professional 
bat handlers and persons who rescue bats may be 
relatively increased because these groups 
primarily handlesick, injured, or orphaned bats. 

Reporting of potential exposures among 
groups at risk changed with time during the 
study. Qne important factor in the management 
of PEP wasthe requirement (introduced in 1997) 
that all bats involved in a potential human 
exposure be surrendered for postmortem exami- 
nation and laboratory testingfor ABL.Thosewho 
care for bats are often reluctant to surrender 
them for ABL testing. Notifications from 
volunteer bat handlers declined during the study 
period. Whilethis may reflect a decline either in 
the number of bat handlers or in potential 
exposures among volunteer bat handlers, 
underreporting may be occurring in this group. 
Anecdotal evidence suggests that this reduction 
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in notifications may reflect handlers' concern for 
the bats. Such underreporting could be associ- 
ated with future human cases. Most potential 
exposures resulted from intentional handling of 
bats. The few potential exposures from unpro- 
voked encounters suggest that bats rarely 
initiate contact with humans. 

The recognition of ABL infection has resulted 
in a large public health program to provide 
education, counseling, and prophylaxis to 
volunteer and professional bat handlers and 
members of the community who may be exposed 
to ABL. The focus of the program has been to 
encourage preexposure vaccination of bat 
handlers, prevention of potential exposures by 
avoidance of bat handling by nonvacci nated 
persons, and prompt medical carewhen potential 
exposures occur. The cost of PEP for all those 
potentiallyexposedtoABL in south Brisbaneand 
the South Coast of Queensland during the study 
was considerable. Future public health interven- 
tions should continue to emphasize the risks 
associated with interaction with bats to reduce 
the requirement for PEP and the likelihood of 
human cases of ABL infection. 
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Genetic Variation in Pneumocystis carinii 
isoiates from Different Geographic 
Regions: impiications for Transmission 
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To study transmission patterns of Pneumocystis carinii pneumonia (PGP) in 
persons with AIDS, we evaluated P. carinii isolates from patients in five U.S. cities for 
variation at two independent genetic loci, the mitochondrial large subunit rRNA and 
dihydropteroate synthase. Fourteen unique multilocus genotypes were observed in 1 91 
isolates that were examined at both loci. Mixed infections, accounting for 17.8% of 
cases, were associated with primary PGP. Genotype frequency distribution patterns 
varied by patients’ place of diagnosis but not by place of birth. Genetic variation at the 
two loci suggests three probable characteristics of transmission: that most cases of PGP 
do not result from infections acquired early in life, that infections are actively acquired 
from a relatively common source (humans or the environment), and that humans, while 
not necessarily involved in direct infection of other humans, are nevertheless important 
in the transmission cycle of P. carinii sp. hominis. 


Pneumocystis carinii pneumonia (PCP) in 
humans, caused by the opportunistic fungal 
pathogen Pneumocystiscarinii f. sp. hominis, isa 
frequent causeof illness in persons with Al DS. I n 
recent years, the number of new cases of PCP has 
greatly diminished (1), largely as a result of 
highly active anti retroviral therapy and effective 
anti-Pneumocystis chemoprophylaxis, primarily 
with tri methopri m-sulfamethoxazole (TM P-SM Z). 
However, PCP is still the most common life- 
threatening AIDS indicator condition in the 
United States in patients whose CD4-Fcel I count 
has deci i ned to <200 cel Is/pL for thefi rst ti me (2). 

Address for correspondence: C. B. Beard, Division of Parasitic 
Diseases, Centers for Disease Control and Prevention, 4770 
Buford Hwy, Mail Stop F-22, Atlanta, GA 30341-3724, USA; 
fax: 770-488-4258; e-mal I : cbbOCgCdc.gov. 


The number of persons at risk, as measured by 
the number of HIV-infected patients with CD4-F 
cell counts <200cells/pL, was recently estimated 
at around 182,000 (3). 

Despite the frequency of PCP, tremendous 
gaps exist in understandingthebasicbiology and 
epidemiology of the causal agent. These gaps, 
which have been reviewed extensively (4-9), 
includelackof knowledge of population structure 
of P. carinii strains circulating among patients or 
in potential environmental reservoirs. The 
implications are important for understanding 
patterns of transmission and developing methods 
of intervention. 

I n the last 5 years, a substantial number of 
genes and gene fragments have been identified 
for potential use in analysis and characterization 
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of P. carinii strains. These loci include the 
mitochondrial largesubunit ribosomal RNAgene 
(mtlsurRNA) (10-12) and the internal tran- 
scribed spacers of the nuclear ribosomal RNA 
array (12-14). The dihydropteroate synthase 
(DH PS) gene locus, which encodes a target for the 
anti-Pneumocystis drugs TMP-SMZ and dap- 
sone, has also been cloned and sequenced (15). 
Substantial variation at this locus (16,17) 
suggests that the widespread use of antimicro- 
bial chemoprophylaxis may be exerting selective 
pressure on P. carinii strains circulating in 
humans. Its polymorphism makes the locus not 
only potent! al I y useful as a marker for changes i n 
antimicrobial susceptibility levels, but also 
valuablefor strain characterization and typing. 

We examined polymorphism at two genetic 
loci of P. carinii isolates from persons with Al DS 
diagnosed in five U.S. cities. One locus, 
mtlsurRNA, is involved in basic metabolic 
functions, and the other, DH PS, is the target of 
sulfoneandsulfonamideantimicrobial drugs. We 
examined the population structure of P. carinii 
strai nsfor i nformation that would el uddatethei r 
patterns of transmission. 

Materials and Methods 

Patient Samples 

Specimens used in the study were obtained 
from March 1995 to J une 1998 during routine 
diagnostic procedures for HIV-infected patients 
hospitalized with POP in Atlanta,Cincinnati, Los 
Angeles, San Francisco, and Seattle. A portion of 
the diagnostic specimen, either induced sputum 
or bronchoalveolar lavage, was preserved 
directly with an equal volume of absoluteethanol 
and stored at 4°C for DNA extraction. 

DNA Purification 

Specimens were divided into approximately 
1-mL aliquots and centrifuged at 14,000 x g for 5 
to 7 minutes. The resulting cell pellet was 
resuspended in 1 mL of phosphate-buffered saline 
(O.OIM, pH 7.2) containing 1 mM EDTA (PBS- 
EDTA), washed twice in PBS-EDTA, centrifuged, 
and stored at -80°C for later DNA extraction. 
DNA was prepared by a commercial purification 
procedure (Wizard Genomic DNA Purification 
Kit, Promega, Madison, Wl) in accordance with 
the product recommendations for DNA purifica- 
tion from blood. Final pellets were resuspended 
inSOpLofTE (10 mM Tris, 1 mM EDTA, pH 7.2). 


Polymerase Chain Reaction (PCR) 
and DNA Sequencing 

PCR amplification was performed at two 
independent genetic loci. A 360-bp fragment was 
amplified from the mtlsurRNA locus by using the 
published primers PAZ102E (5' -GAT GGCTGT 
TTC CAA GCC CA - 3') and PAZ102H (5' - GTG 
TAC GTT GCA AAG TAC TC - 3') (10). PCR 
conditions included a 94°C hot start for 5 
minutes; followed by 35 cycles of a program 
consisting of 92°C for 30 seconds, 55°C for 30 
seconds, and 72°C for 60 seconds; followed by a 
termination step at 72°C for 5 minutes. The 
DHPS locus was amplified by a modification of a 
nested PCR procedure (16,17). I n the first round 
of this PCR, the pri mers DH PS F 1 (5' - CCT GGT 
ATT AAA CCA GTT TTG CC - 3') (S.R. Meshnick, 
pers. comm.) and DHPS B 45 (5' - CAA TTT AAT 
AAA TTT CTT TCC AAA TAG CAT C - 3') (16) 
were used. In the second round, the primers 
DH PS Ahum (S' - GCG CCT ACA CAT ATT ATG 
GCC ATT TTA AAT C - 3') and DHPS BN (5' - 
GGA ACT TTC AAC TTG GCA ACC AC - 3') (16) 
were used. Theconditions for thefirst round were 
94°C for 5 minutes, followed by 35 cycles of 92°C 
for 30 seconds, 52°C for 30 seconds, 72°C for 60 
seconds, and a termination step at 72°C for 5 
minutes. The conditions for the second round of 
PCR were 94°C for 5 minutes, then 35 cycles of 
92°C for 30 seconds, 55°C for 30 seconds, and 
72°C for 60 seconds, followed by a termination 
step at 72°C for 5 minutes. 

For analysis of the PCR-generated frag- 
ments, 10 pL of each 50-pL PCR amplification 
product was examined by horizontal gel 
electrophoresis on 1% agarose gels. The 
remaining 40 pL of successfully amplified 
reactions was purified by a commercial purifica- 
tion procedure (Wizard PCR Purification Kit, 
Promega, Madison, Wl) and suspended in 50 pL 
of TE for DNA sequencing; 5 pL to 10 pL of each 
purified product was sequenced directly by using 
dye terminator chemistry (ABI Prism BigDye 
Terminator Cycle Sequencing Ready Reaction 
Kit, PE Applied Biosystems, Foster City, CA) 
according to the manufacturer's protocol, with 
the DNA oligonuceotide primers used for PCR 
amplification. The PCR fragments were se- 
quenced on an ABI 377 automated DNA 
sequencer (PE Applied Biosystems, Foster City, 
CA) according to manufacturer's recommenda- 
tions. The sequenced DNA fragments were 
analyzed by using Sequence Navigator v.1.0.1 
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(PE Applied Biosystems, Foster City, CA)andthe 
GCG Wisconsin Package Version 9.1. computer 
program (18). 

Statistical Analysis 

Statistical analyses were performed with 
SAS software version 6.12. Logistic regression 
was used to examine the association between 
pi ace of di agnosi s and pi ace of bi rt h and genotype 
for mtlsurRNA and DHPS gene loci. 

Results 

Amplification with Specific Primers 

PCR amplification with the selected primer 
sets gave consistent results in most P. carinii- 
positive samples. The mtlsurRNA primers 
amplified a 360-bp fragment in 223 samples from 
Atlanta, Los Angeles, San Francisco, and Seattle. 
The DFIPS primer sets amplified a 300-bp 
fragment in 220 of these samples. Both genetic 
loci were successfully amplified for 191 samples. 
In the Cincinnati dataset, 101 samples were 
amplified for the mtlsurRNA locus. The DFIPS 
site was not examined for this group of samples. 
All the amplified fragments were sequenced, 
al i gned, and exami ned for genetic polymorphi sm. 

General Observations on Genotype Frequency 

Four unique genotypes were observed for 
each of the two genetic loci examined. At the 
mtlsurRNA locus, all the genotypes were 
distinguished on the basis of polymorphism at 
nucleotide positions 85 and 248. Genotypes 1 
(85:C; 248:C) and 2 (85:A; 248:C) were the most 
common (Table 1), occurring at similar frequen- 
cies and together accounting for 74.7% of the 324 
samples analyzed. Genotypes 3 (85:T ; 248:C) and 
4 (85:C; 248:T) accounted for 9.3% and 5.9%, 
respectively. At the DFI PS locus, genotypes 1 and 
4 accounted for 80.4% of all samples (Table 2). 
Genotypes 2 and 3 were both relatively 
uncommon, seen in 5.9% and 2.3% of the samples 


Table 1. Genotype analysis of Pneumocystis carinii 
isolates at the mtlsurRNA locus 


Genotype 

Nucleotide 

position/identity 

F requency 
n=324 (%) 

1 

85/C; 248/C 

123 (38.0) 

2 

85/A; 248/C 

119 (36.7) 

3 

85/T; 248/C 

30 (9.3) 

4 

mixed 

85/C; 248/T 

19 (5.9) 
33 (10.2) 


Table 2. Genotype analysis of Pneumocystis carinii 
isolates at the DHPS locus 


Genotype 

Nucleotide (ami no acid) 
position/identity 

F requency 
n=220 (%) 

1 

165 (55)/A (Thr) 

68 (30.9) 

2 

171 (57)/C (Pro) 
165 (55)/G (Ala) 

13 (5.9) 

3 

171 (57)/C (Pro) 
165 (55)/A (Thr) 

5(2.3) 

4 

171 (57)/T (Ser) 
165 (55)/G (Ala) 

109 (49.5) 

mixed 

171 (57)/T (Ser) 

25 (11.4) 


analyzed, respectively. All the mutations seen at 
this locus were nonsynonymous changes result- 
ing in amino acid substitutions; no other 
polymorphism was observed. Genotype 1 was the 
designation used to refer to the sequence defined 
by a threonine at position 55 and a proline at 
position 57. Genotype 2 referred to an alanine at 
55 and a proline at 57, Genotype 3 to a threonine 
at position 55 and a serine at 57, and genotype 4 
to an alanine at position 55 and a serine at 
position 57. 

When the results at both genetic loci were 
combined, 14 unique multi locus genotypes of 16 
possible combinations were observed in 191 
samples for which both genes could be amplified 
and sequenced. Thefour most common multi locus 
genotypes accounted for 61.7% of all genotypes. 
The most common multi locus genotypes con- 
sisted of combinations of the most common 
genotypes at each individual locus. No genetic 
linkage of specific genotypes from the two loci 
was observed. 

Coinfection with multiple P. carinii strains 
could be detected in 33 (10.2%) of 324 samples 
typed at the mtlsurRNA locus and 25 (11.4%) of 
220 typed at the DFIPS locus. When the 
genotypes were considered together, 34 (17.8%) 
of 191 samples represented coinfections with 
multiple genotypes. When these samples were 
analyzed according to patient history, 21 (17.6%) 
of 119 samples from patients with no history of 
previous PCP were coinfected with multiple 
strains, by typing at the mtlsurRNA locus. In 
contrast, none of the samples from 41 patients 
with PCP history had multiple genotypes (p = 
0.002, Fisher's exact test). A similar trend was 
observed at the DFIPS locus but was not 
statistically significant. 
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Genotype Frequency Distribution 
Patterns by Geographic Location 

Genotype frequencies in PCP isolates 
differed by city (Figures 1, 2). At the mtlsurRN A 
locus (Figure 1), the genotype distribution was 
significantly different (chi-square test; p =0.001), 
with the rati 0 of genotype 1 to genotype 2 rangi ng 
from 0.7 in San Francisco to 1.8 in Cincinnati. 
Anal ysi s of t he mtl su rR N A genotype f requenci es 
at each city showed that genotype distributions 
also differed significantly or borderline signifi- 
cantly, with p values from 0.001 to 0.08, with the 
exception of Los Angeles. Comparisons involving 
Los Angeles were nonsignificant, perhaps as a 
result of small sample size (n =15). 


Figurel. Distributionof Pnajmocystiscarinii mtlsurRNA 
genotypes, by city (chi -square test, p = .001) 

Atlanta, n =76; Cincinnati, n =107; LosAngeles, n =15; San 
Francisco, n =68; Seattle, n =57. 



Figure 2. Distribution of Pneumocystis carinii DHPS 
genotypes by city (Fisher's exact test, p=.049). 
Atlanta, n=80; Los Angeles, n=20; San Francisco, n=66; 
Seattle, n=53. 


We used logistic regression to investigate 
whether the geographic distribution was associ- 
ated more with place of birth or with place of 
diagnosis. We grouped places of birth and 
diagnosis into either "east" or "west," using the 
M ississippi River as the dividing line. Data from 
Cincinnati were not included, sinceplace-of-birth 
information was not aval I able there. None of the 
four mtlsurRNA genotypes was significantly 
associ ated with pi ace of bi rth when the data were 
adjusted for pi ace of diagnosis. Flowever, genotypes 
2 and4weresignificantly associated with pi ace of 
diagnosis, when the data were adjusted for place 
of birth (p =0.002, p =0.05, respectively); the 
association of genotype 1 with place of diagnosis 
was borderline significant (p =0.08). 

The overall genotype distribution at the 
DFIPS gene locus (Figure 2) also differed 
significantly (Fisher's exact test; p =0.045) for 
the various cities. The ratios of genotype 4 to 
genotype 1 ranged from 0.97: 1 i n Atlanta to 3. 13:1 
in San Francisco. No isolates from Cincinnati 
were analyzed at the DFI PS locus. 

Logistic regression analysis of the DFIPS 
genotypes by pi ace of di agnosi s and pi ace of bi rth 
showed results similar to those for the 
mtlsurRNA genotypes. Genotypes 1 and 4 were 
associated with place of diagnosis (p =0.001, p = 
0.002, respectively) when data were adjusted for 
pi ace of birth. In contrast, pi ace of birth was not 
associated with any DFIPS genotype when data 
were adjusted for place of diagnosis. 

Discussion 

P. carinii, once thought to be a protozoon but 
now regarded an ascomycetel i kefungus (4,19,20), 
has been associated with human disease since 
the 1940s. Despite intense efforts to understand 
this important disease agent, lack of a suitable 
means of propagation has complicated study of 
basic biology and epidemiology. Similarly, 
latency and reactivation (versus recent acquisi- 
tion) have not been resolved. PCR and other 
molecular methods have improved understand- 
ing of genetic variability and host specificity 
(11,12,14,20,21), as well as thecause of recurrent 
infections in persons with Al DS (8,22-24). 

Both an indirect (environmental) source and 
a direct (person-to-person) source have been 
proposed as modes of transmission of P. carinii in 
humans. Three primary observations reported in 
this study address infection sources and 
transmission patterns: geographic variation in 
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genotype frequency distribution, which was 
correlated with the pi ace of diagnosis rather than 
the place of birth; rate of coinfection with 
multiple P. carinii genotypes, which was greater 
in patients with primary rather than secondary 
PCP; and abundance oftheDH PS doublemutant, 
which accounted for 49.5% of all DHPS 
genotypes, strongly suggesting genetic selection. 

Allelic frequency distribution patterns of P. 
carinii isolates were associated with place of 
diagnosis rather than place of birth. Place of 
diagnosis was determined more reliable than 
pi ace of residence as the best i ndicator of the most 
I i kel y pi ace of exposu re to P . ca r i n i i because pi ace 
of residence was frequently recorded as 
permanent or legal residence rather than as 
pi ace of resi dence du ri ng the 3- to 6-month per i od 
before admission. These results suggest that 
infection in adults is acquired later than thefirst 
few months or years of life and that any latency 
has natural limits. Two independent lines of 
study offer a context for these observations and 
suggest that PCP is an actively acquired 
infection. The first line of study comprises 
molecular genetic analyses of P. carinii strains 
from adult AIDS patients with recurrent PCP. 
These analyses have shown that different 
P. carinii genotypes are detected on subsequent 
PCP episodes in a substantial proportion of 
patients (22,23), which suggests that infections 
in adults are actively acquired and that 
subsequent infections do not necessarily repre- 
sent relapses. Thesecond lineof study comprises 
observations of primary PCP in infants with 
perinatally acquired HIV infection (25) who had 
PCP episodes at 3 to 6 months of age; the 
prevalence of PCP among these infants was very 
similar to that among adults with Al DS. These 
observations suggest that P. carinii is common in 
the environment, consistent with the suggestion 
that the organism is easily acquired. 

While the association with geographic 
distribution is consistent with recent acquisition 
of clinical infections, the actual source of the 
infection is not known. The rat P. carinii model 
suggests direct (animal-to-animal) transmission 
(5,26,27); however, if person-to-person transmis- 
sion does occur, it has not been shown to be of 
epidemiologic significance (28,29). PCR amplifi- 
cation of P. carinii DNA from spore traps from 
various sites supports the possibility of an 
environmental source (30-33), but no specific 
plant, animal, or soil source has been identified. 


Identification of a specific environmental reser- 
voir for P. carinii f. sp. hominis, if one exists, 
could elucidate disease transmission and im- 
prove prevention efforts. 

The rate of coinfection with multiple 
P. carinii strains by multi locus typing was 17.8% 
(34 of 191). Of the 191 patients for whom both 
genes were typed, 20 (10.5%) were coinfected at 
the mtlsurRNA locus and 19 (9.9%) at the DH PS 
locus. Five patients tested positive for multiple 
strains concurrently at both loci. While the 
sensitivity to detect coinfections was similar for 
both genes, concordance between the two genes 
was very low in detecting coinfections with 
multiple strains in any particular isolate— only 
four genotypes can be detected at each individual 
locus. When the loci are considered together, 
however, the sensitivity is greatly increased, 
because the number of possible genotypes 
i ncreases from 4 to 16, 14 of which were observed 
in this study. 

The 17.8% coinfection rate is within the 10% 
to 30% range reported by other investigators 
using DNA sequencing-based approaches 
(14,22,34). Although coinfection rates as high as 
69% have been reported with single-strand 
conformation polymorphism analysis (35), these 
differences can probably be resolved by two 
considerations. The first is the sensitivity of the 
selected locus, which is a function of the genetic 
variability or evolutionary rate of the locus (i.e., 
the more variable the locus, the greater its 
sensitivity to detect a different genotype). I n this 
study, mtlsurRNA and DHPS displayed a high 
degree of genetic conservation, with only four 
alleles detected at each locus, all of which have 
been observed (11,16,17). Because these are 
relatively slowly evolving genes, the sensitivity 
to detect coinfections with multiple strains is 
expected to be lessthan in a faster evolving locus. 
The second consideration is PCP prophylaxis. 
Correlation between primaryPCP and coinfection 
with multi pie strains suggests that patients may 
be exposed to multiple P. carinii strains over 
extended periods and harbor short-lived, latent 
infections. Consequently, patients who have not 
been treated are more I i kely to have been exposed 
(andtoharbor) multi pleP. carinii strains. On the 
other hand, patients who have been treated for 
PCP at least once and are taking secondary 
prophylaxis are less likely to become reinfected. 
This hypothesis is consistent with indications 
that P. carinii is ubiquitous in the environment 
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and exposure i n humans i s commonpl ace (22,23). 
A short latency period has been suggested (36) 
and is consistent both with coinfections' 
correlating with primary PCP and allelic 
frequency distributions' correlating with the 
place of diagnosis but not the pi ace of birth. 

The question of latency and reactivation 
versus recent acquisition is of great importance 
as it relates to prevention. I f the preponderance 
of infections in humans results from activated 
latent infections, chemoprophylaxis is the only 
method of preventing disease. If, however, 
infections are actively acquired, identifying 
specific sources of infection would lead to other 
methods of prevention and reduce dependence on 
antimicrobial agents. 

The most important clinical implication of 
the polymorphism observed specifically at the 
DH PS locus relates to the emergence of possi ble 
antimicrobial resistance; however, the observa- 
tion of different DHPS alleles in P. carinii 
populations also has implications for transmis- 
sion routes and patterns and the possibility of 
person-to-person transmission. 

In human P. carinii isolates, four distinct 
genotypes have been reported at the DH PS locus 
(16,17). All fouroftheseDNA base changes result 
in amino acid substitutions at an important 
structural position, in an otherwise very highly 
conserved gene (17,37). Genotype 1 in our study 
corresponds to the "wild-type" genotype, as 
defined by the only allele observed in P. carinii 
from any other mammalian species (16,17). 
Genotype 2 is a point mutation that results in a 
threonine-to-alanine substitution at amino acid 
position 55. Genotype 3 is also a point mutation, 
resulting in a proline-to-serine substitution at 
amino acid position 57. Genotype 4 is a double 
mutant, which has alanine at position 55 and 
serine at position 57. This particular mutation 
has rarely been associated with failure of both 
TMP-SMZ treatment (38) and prophylaxis (17). 
The presence of polymorphism at two different 
positions, both associated with amino acid 
substitutions, at a gene locus that otherwise 
lacks variability, is itself indicative of pro- 
nounced selective pressure. Such a degree of 
selection is not itself an indication of 
resistance, but certainly suggests that resis- 
tance may beemerging. M utations in this same 
region of the DHPS molecule have been 
correlated with specific antimicrobial resis- 
tance to sulfa drugs in a number of other 


microorganisms, including Plasmodium falci- 
parum (39), Streptococcus pneumoniae (40), 
Streptococcus pyogenes (41), Escherichia col i (42), 
and Neisseria meningitidis (43). 

The two single mutations are uncommon by 
themselves (5.9% for genotype 2 and 2.3% for 
genotype 1), yet 50% of all isolates have the 
double mutant, which suggests that, alone, 
neither mutation is highly selected, but together 
they pose a very strong selective advantage. If 
humans are dead-end hosts for P. carinii, how 
can such a pronounced degree of polymorphism 
be explained? Person-to-person transmission 
may be essential to allow genetic selection to 
occur, resulting in this polymorphism. 

Geographic variation in allelic frequency 
detected at the mtlsurRNA locus is not 
unexpected, even though other studies with 
smaller sample sizes failed to observe differences 
(11). Mitochondrial DNA sequence data have 
been highly useful in detecting intraspecific 
differences between populations of diverse 
organisms (44,45). Perhaps more unexpected, 
however, was the detection of geographic 
variation attheDHPS locus. UnlikemtIsurRNA, 
because it encodes a gene product that is the 
target of the primary anti-P. carinii drug TM P- 
SM Z, the DHPS locus is assumed to be subject to 
intense selection pressure. Consequently, this 
selection might be expected to override any 
potential variation fromgeographicseparation or 
geographic patterns to reflect TMP-SMZ expo- 
sure patterns. The relative frequency of the 
double mutant genotype was much higher in 
samples from the West Coast, particularly San 
Francisco, than in samples from Atlanta 
(Figure 2). The reason for this difference is not 
obvious because in preliminary multivariate 
analysis conducted with samples from Atlanta 
and San Francisco, when the data are 
controlled for sulfa exposure, pi ace of diagnosis 
is still the most significant factor influencing 
the frequency of the double mutant genotype 
(data not shown). Explanations for this 
observation are being evaluated. 

An overall frequency of approximately 50% 
for the DHPS double mutant genotype is 
somewhat higher than that reported in other 
recent studies (46,47). The reason for this 
observation is also unclear. One possible 
explanation is that the samples were collected 
more recently (March 1995 to J une 1998) and 
therefore reflect a sulfa-induced increase in 
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frequency of the double mutant. Most patients 
were from San Francisco and Atlanta, the two 
populations that showed the greatest difference 
ingenotypefrequencyattheDHPS locus. Sample 
collection in Atlanta began in 1995, but in San 
Francisco the earliest samples were from 1997. 
When we stratified the patients in Atlanta by 
date, 1995-96 versus 1998-99, we saw no 
significant difference in the frequency of mutant 
genotypes. Similarly, when we stratified data for 
San Francisco by date, 1997 versus 1998-99, we 
saw no differences. Consequently the frequency 
of the mutation did not changeduring this study 
period, nor does the variation observed according 
to geography appear to be confounded by the 
dates when specimens were collected. 

Conclusions 

The pattern of allelic variation differed at 
each of the cities where samples were obtained, 
and this variation correlated with the place of 
diagnosis but not with the place of birth. 
Coinfection with multiple P. carinii genotypes 
was associated with primary rather than 
secondary PCP. The position 55/57 double 
mutant accounted for 50% of al I genotypes of the 
DFIPS locus examined. These observations 
suggest thefollowing possibilities, which contra- 
dict much of the current opinion on the 
epidemiology and transmission of PCP: Most 
cases of PCP are not a result of infections 
acquired very early in life; infections are actively 
acquired from a relatively common source 
(humans or the environment); and humans, 
while not necessarily involved in direct infection 
of other humans, are nevertheless important in 
thetransmission cycleof P. carinii f. sp. hominis. 
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Rhinosporidium seeberi: A Human 
Pathogen from a Novel Group of Aquatic 
Protistan Parasites 
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Rhinosporidium seeberi, a microorganism that can infect the mucosal surfaces of 
humans and animals, has been classified as a fungus on the basis of morphologic and 
histochemical characteristics. Using consensus polymerase chain reaction (PGR), we 
amplified a portion of the R. seeberi 18S rRNA gene directly from infected tissue. 
Analysis of the aligned sequence and inference of phylogenetic relationships showed 
that R. seeberhs a protist from a novel clade of parasites that infect fish and amphibians. 
Fluorescence in situ hybridization and R. seeber/'-specific PGR showed that this unique 
18S rRNA sequence is also present in other tissues infected with R. seeberi. Our data 
support the R. seeber/phylogeny recently suggested by another group. R. seeberhs not 
a classic fungus, but rather the first known human pathogen from the DRIPs clade, a 
novel clade of aquatic protistan parasites (Ichthyosporea). 


Rhinosporidiosis manifests as slow-growing, 
tumorl ike masses, usually of the nasal mucosa or 
ocular conjunctivae of humans and animals. 
Patients with nasal involvement often have 
unilateral nasal obstruction or bleeding due to 
polyp formation. The diagnosis is established by 
observing the characteristic appearance of the 
organism i n tissue biopsies (Figure 1). T reatment 
consists of surgical excision, but relapse occurs in 
approximately 10% of patients (1); antimicrobial 
therapy is not effective (2). Rhinosporidiosis 
occurs in the Americas, Europe, Africa, and Asia 
but is most common in the tropics, with the 
highest prevalence in southern India and Sri 
Lanka. A survey of schoolchildren from Pal lam, 
India, found 11 cases in 781 children examined 
(prevalence 1.4%) (3). Autochthonous cases have 
been reported from the southeastern United 
States (4). Studies have linked infection to 
swimming or bathing in freshwater ponds, lakes, 
or rivers (2,5). 

The etiologic agent of rhinosporidiosis, 
Rhinosporidium seeberi, is an enigmatic microbe 
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Ave, Paio Aito, CA 94304, USA; fax: 650-852-3291; e-maii: 
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that has been difficult to classify. Recently, 
R. seeberi has been considered a fungus, but it was 
originally thought to be a protozoan parasite (2). 
Its morphologic characteristics resemblethoseof 
Cocci dioides immitis: both organisms have 
mature stages that consist of large, thick-walled. 





^ 

r- ■' ' ‘ . • 

\!' ■ , c: .... 

■ . 





Figure 1. Histology of rhinosporidiosis. A formalde- 
hyde-fixed section of human nasal polyp was stained 
with periodic acid-Schiff (PAS) and visualized by 


bright-field microscopy at 400X magnification. The 
thick walls of immature R. seeberi trophocytes stain 
with PAS (pink), and the spherical organisms are 
surrounded by inflammatory cells. 
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spherical structures contai ni ng smal I er daughter 
cells (endospores). In addition, R. seeberi is 
visualized with fungal stains such as methe- 
namine silver and periodic acid-Schiff, as well as 
mud carmine, which stai ns thefungus Cryptococ- 
cus neoformans. R. seeberi has not been detected 
in the environment, and its natural host or 
reser voi r i s u n known . Attempts to propagate this 
organism on artificial media have failed, as has 
continuous cocultivation with human cell lines (6). 

We report a molecular approach for estab- 
lishing the phylogenetic relationships of R. 
seeberi to other eukaryotes. This approach is 
based on amplification of thesmall subunit rRNA 
gene sequence from infected tissue, as in the 
method used to identify the culture-resistant 
bacillus of Whipple disease (7). The sequence of 
the small subunit rRNA gene has proven to be a 
useful gaugeof evolutionary relationships for many 
organisms from diverse taxonomic groups (8). 

Materials and Methods 

Specimens 

We obtained a sample of frozen, minced, 
infected canine nasal polyp in tissue culture 
media that had been used for the limited 
propagation of R. seeberi in cell culture (6). After 
thawing and a 1-minutecentrifugation at 500 xg, 
approximately 0.2 g of the tissue pellet was 
digested by mechanical disruption and the DNA 
was purified by adsorption to glass milk in the 
presence of a chaotropic agent, according to the 
manufacturer's instructions (Fast Prep, Bio 101, 
Vista, CA). The DNA was resuspended in 100 pL 
of 10 mM Tris, 1 mM EDTA buffer at pH 8.5. 

Three blocks of fixed, paraffin-embedded 
nasal polyps from unrelated patients with 
histologically confirmed rhinosporidiosis were 
obtained. One tissue sample came from a patient 
born in southern Asia but living in the United 
States, and two samples came from patients 
living in Spain. Twenty-three blocks of fixed, 
paraffin-embedded nasal polyps from patients 
without rhinosporidiosis were obtained from 12 
consecutive patients at the Palo Alto Veterans 
Affairs (VA) hospital who had undergone nasal 
polypectomy. Two 25-pm sections were cut from 
each block, and the sections were deparaffini zed 
and digested (5). DNA from the digests was then 
purified by the Isoquick method (ORCA 
Research, I nc., Bothell, WA), and the DNA was 
resuspended in 25 pL of Tris/EDTA buffer. 


Sectionsfroma block of fixed, paraffin-embedded 
human lymph node (histologically normal) were 
also us^ in some experiments as negative 
controls. Tissue sections of C. immitis in bone 
were obtained from a patient with disseminated 
coccidioidomycosis at the Palo Alto VA hospital 
and used for fluorescence in situ hybridization 
(FISH). Samples of Rosette agent and 
Dermocystidium salmonis DNA were obtained 
from the Bodega Marine Laboratory of the 
University of California, Davis. 

Fungal specimens in culture were obtained 
from laboratories at Stanford U niversity and the 
Palo Alto VA Health Care System. A cotton swab 
was used to transfer fungal cel Is from agar into a 
1.5-mL microfuge tube containing 0.5 mL 
digestion buffer (9) and 0.1 mL glass beads. 
Samples were incubated at 55°C overnight, and 
the proteinase k was inactivated at 95°C for 10 
minutes and then subjected to two freeze-thaw 
cycles by immersing tubes in a dry ice- 
isopropranol bath followed by vortex mixing. 

Consensus Polymerase Chain Reaction (PCR) 
of the 18S rRNA Gene 

Broad-range fungal PCR primers were 
designed from a database of >4,000 smal I subunit 
rDN A sequences, with the ARB software package 
(Technical University, Munich, Germany). Prim- 
ers were selected that would anneal to most 
fungal and some protist 18S rDN A but not to 18S 
rDNA from the chordata (Fl-fw, F2-rev, F3-rev) 
(Table). When the specificity of the primer pairs 
was tested by using human lymphocyte DNA, no 
amplification was observed (data not shown). The 
broad range of primers Fl-fw/F 2-rev was tested 
by using DNA from several diverse fungi. 
Amplification products of the expected size were 
produced by using DNA from Aspergillusoryzae, 
Alternaria alternata, Candida albicans, Saccha- 
romycescere/isiae, T richyphyton rubrum, Panus 
rudis, Neurospora crassa, Fusarium solani, 
Beauveria bassiana, Flammulina velutipes, 
Gibberellazeae, and Pleurotusostreatus(data not 
shown). 

PCR con si St ed of 40 cycl es of a m pi i f i cat i on on 
a Perkin-Elmer GeneAmp 2400 thermal cycler. 
After an initial activation of Tag gold at 94°C for 
10 mi nutes, each cycle consisted of 30 seconds of 
melt! ng at 94°C, 30 seconds of anneal i ng at 56°C, 
and 30 seconds of extension at 72°C. The last 
cycle was followed by an extension step at 72°C 
for 7 minutes. Amplification products were 


Emerging Infectious Diseases 


274 


Vol. 6, No. 3, May-]une 2000 



Research 


Table. Polymerase chain reaction primers, sequencing primers, and fluorescence in situ hybridization probes 


Primer/Probe 

Nucleotide sequence (5'->3') 

SSU rRNA Position^ 

F 1-fw^ 

F2-reV= 

F 3-rev 

Dermo-fw 

Dermo-rev 

Rhino-fw 

Rhino-rev 

Rhino FISFI probe 

Control FISFI probe 

CAAGTCTGGTGCCAGCAGCC 

GATTTCTCGTAAGGTGCCGA 

AATGCTCTATCCCCAMCACG 

CTGCCAGTAGTCATATGCTTG 

GATCAAGTTTGATCAACTTTTCGGCA 

GGCGTGTGCGCTTAACTGTC 

TGCTGATAGAGTCATTGAAT 

BTGCTGATAGAGTCATTGAATTAACATCTACB 

BACGACTATCTCAGTAACTTAATTGTAGATGB 

554-573 

1068-1087 

1531-1550 

^SSU rRNA pceition based on S. 

cere/isiae 18S rRNA (GenBankJ 01353). 



•Tw = forward primer. 
Tev = reverse pri mer. 


detected by electrophoresis on 2% agarose gels 
stained with ethidium bromide and visualized 
with a UV transilluminator. 

On the basis of the sequences obtained by 
consensus PCR with primers Fl-fw/F 2-rev (-500 
bp) and F 1-fw/F 3-rev (-1000 bp), pri mers Dermo- 
fw and Dermo-rev were designed (Table) and 
used in a PCR to amplify a more complete 185 
rDNA sequence of R. seeberi. 

Rhinosporidium-Specific PCR 

A pair of PCR primers was designed from 
unique regions of the R. seeberi 185 rRNA gene 
sequence (Rhino-fw and Rhino-rev) (Table). 
These primers were used in a 50-pL PCR as 
described, except that AmpliTaq DNA poly- 
merase (PE-ABI ) was used at 1 unit per reaction, 
no dimethyl sulfoxide was added, 50 cycles of 
PCR were run with a 3-mi nute pre-melt at 94°C, 
and the annealing temperature was 55°C. To 
each 50-pL PCR reaction, 1 pL or 5 pL of purified 
DNA were added. 

B-Globin PCR 

15-globin PCR was performed on control 
tissues as described previously (9). 

DNA Cloning, Sequencing, 
and Phylogenetic Analysis 

Amplification products were cloned by using 
theTopo-TA cloning kit (I nvitrogen, Carlsbad, 
CA), andthreeclonesweresequenced. Each clone 
consisted of 1,750 bp of 185 rDNA. Priming 
sequences were removed for further analysis, 
yielding 1,699 bp of meaningful sequence. DNA 
sequencing was performed as described (10). A 
consensus sequence from the three clones was 


made to correct for any Taq polymerase 
incorporation errors. The 185 rDNA primers 
(Table) were used as sequencing primers. 

The R. seeberi 185 rDNA sequence was 
aligned by using the automated aligner of the 
ARB software package. Ambiguously and 
incorrectly aligned positions were manually 
aligned on the basis of the conserved primary 
sequence and secondary structure. The phyloge- 
netic relationship of R. seeberi to other 
eukaryotes was inferred from 1,350 unambigu- 
ously aligned (masked) positions with a 
maximum-likelihood algorithm (11,12), on the 
basis of a previously aligned dataset of the DRI Ps 
clade (named after the organisms Dermocystidi- 
um, the Rosette agent, Ichthyophonus, and 
Psorospermium) (13). The dataset was used to 
empirically determine nucleotide frequencies 
and instantaneous substitution rates with the 
restriction of a 2:1 transition to transversion 
ratio. The organisms used in our tree and the 
accession numbers for their small subunit rRNA 
sequences include Artemi a salina (X01723), 
Xenopus I ae/i s (X04025), M yti I us edu I i s (L 24489), 
Tripedalia cystophora (L 10829), Microciona 
prolifera(L 10825), Diaphanoecagrandis(L10824), 
Rosette agent (L29455), R. seeberi (AF 158369), 
Dermocystidium species (U21336), Dermocystid- 
i um sal moni s (U 21337), Psorospermi um haeckel i i 
(U 33180), Ichthyophonus hoferi (D14358), As- 
pergillus fumagatus (M60300), Chytridium 
confervae (M 59758), M ucor racemosus (X54863), 
Acanthamoeba castellanii (U 07413), Zamia 
pumila (M20017), Porphyra spiralis (L26177), 
Lagenidium giganteum (X54266), Labyrinthu- 
loides minuta (L27634), Perkinsus marinus 
(X75762), 5arcocystis muris (M64244). The tree 
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topology was confirmed by using a neighbor- 
joining algorithm with J ukes-Cantor corrected 
distance values and a maximum-parsimony 
algorithm (ARB).Thenucleotidesequenceforthe 
partial 185 rRNA gene of R. seeberi has been 
deposited in GenBank (accession number 
AF 158369). 

Fluorescence in Situ Hybridization (FISH) 

Tissue sections on slides were dewaxed by 
immersion in 99% octane (Sigma, St. Louis, MO). 
Samples subjected to FISH included R. seeberi- 
infected human nasal polyps, C. immitis-infected 
bone, a Rosette agent- infected cell line, and 
smea rs of C . a I bi ca ns. T he R h i nospor i d i u m probe 
was based on the Rhino-rev 18S rDNA primer 
and was biotinylated at both the 5' and 3' ends 
(Table). Thecontrol probe, which consisted of the 
complement of the Rhinosporidium probe, was 
also biotinylated at both ends. To each slide, 50 
ng of biotinylated probe in 30 pL of hybridization 
buffer was added. Cover slips were placed, and 
the slides were incubated at 40°C overnight in a 
humid chamber. The hybridization buffer 
consisted of 10% dextran, 0.2% bovine serum 
albumin, and 0.01% polyadenosine, in 5X SET 
buffer; the 25X SET buffer consisted of 3.75M 
sodium chloride, 25 mM EDTA, and 0.5M T ris at 
pH 7.8. Cover slips were removed by immersion 
in 5X SET buffer at 4°C, and the slides were 
washed for 10 minutes per cycle, twice in 0.2X 
SET buffer at 25°C and once at 40°C. The slides 
werethen subjected totyramide signal amplifica- 
tion accordingtothe manufacturer's instructions 
(TSA indirect, NEN Life Sciences, Boston, MA). 
Cy5-streptavidin (Amersham, Piscataway, Nj ) 
at 1 mg/mL was diluted 1:500 and added to the 
slides for fluorescence signal detection. Tissue 
sections were visualized on a Bio-Rad confocal 
microscope at 200X magnification after the 
application of 15-20 pL of Vectashield mountant 
(Sigma) and a cover slip. 

Electron Microscopy 

A portion of formalin-fixed, paraffin-embed- 
ded nasal polyp from a patient with rhinosporidi- 
osis was removed from the block, dewaxed with 
xylene, rehydrated with ethanol, post-stained 
with 1.5% osmium tetroxide, then dehydrated 
with ethanol, transferred to propylene oxide 
followed by Epon 12 resin, heat-catalyzed at 
65°C, and ultrasectioned at 50 nm. The grid- 
mou nted secti ons were then seri al I y stai ned with 


lead hydroxide and uranyl acetate and examined 
with a Phillips 201 electron microscope at 75KV. 

Results 

Phylogenetic Classification of R. seeberi 
Inferred from the 18S rRNA Gene 

Consensus PCR of the 185 rRNA gene with 
DNA from a digest of an R. seeberi -infected 
caninenasal polyp produced amplification products 
of the expected size visible on gel electrophoresis 
(pri mer pai rs F 1-fw/F 2-rev = -500 bp and F 1-fw/ 
F 3-rev = -1000 bp) (data not shown). No 
amplification product was detected by using 
control tissue and reagents. On the basis of the 
initial phylogenetic assessment of these se- 
quences, primers Dermo-fw and Dermo-rev were 
designed for amplification of a more complete 
portion of the R. seeberi 18S rRNA gene. Our 
phylogenetic analysis of this gene suggests that 
R. seeberi is a member of the DRIPs clade of 
aquatic protistan parasites (Figure 2). The 
nearest evolutionary neighbors of R. seeberi for 
which a sequence is available are members of 
the Dermocystidium genus, which infect 
salmon and trout. 

Development and Use of a 
Rhinosporidium-Specific PCR Assay 

A PCR assay specific for R. seeberi was 
developed. Primers were created by aligning 18S 
rDNA sequences from R. seeberi, members of the 
DRIPs clade, Saccharomyces cero/isiae, and 
humans. The Rhinosporidium primers (Rhino- 
fw. Rhino-rev) each have three nucleotide 
mi smatches with the sequences from the nearest 
phylogenetic relatives in the Dermocystidium 
genus and multiple other mismatches with 
fungal and human 185 rDNA sequences. An 
assay sensitivity of 1-10 gene copies was 
demonstrated by using a dilution series of cloned 
R. seeberi 185 rDNA. The specificity of the assay 
was assessed with DNA from human lympho- 
cytes, S. cere/isiae D. salmonis, and the Rosette 
agent. No amplification was detected when these 
DNA samples were used in the specific PCR 
assay, although product was amplified from 
these samples with either (5-globin primers 
(lymphocytes) or broad-range 185 rDNA primers 
(F 1-fw/F 2-rev) (Table). 

The original DNA sample from the infected 
cani ne nasal polyp yielded a visi ble product of the 
expected size with the Rhinosporidium-specific 
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Figure 2. Phylogeny of Rhinospo- 
ridium seeberi and the DRI Ps clade 
of protists (I chthyosporea). A. Phy- 
logenetic tree inferred from the 18S 
rDNA sequences of R. seeberi and 
other selected eukaryotes by using a 
maximum likelihood algorithm; 
1,350 masked positions were used 
for analysis. Bootstrap values were 
generated from 100 resamplings. 
The bar, which represents 0.1 base 
changes per nucleotide position, is a 
measu re of evol uti onary di stance. B . 
Phylogenetic tree using the data 
from A, but with pruning and 
grouping to show the broader 
evolutionary position of the DRIPs 
clade. 
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PCR assay (Figure 3). Purified DNA from the 
tissue blocks of human nasal polyps resected 
from three patients with rhinosporidiosis also 
yielded positive results in this PCR assay, with 
visible bands seen on gel electrophoresis (RS 1-3) 
(Figure 3). Direct sequencing of these PCR 
products demonstrated complete identity over 
the 377 bp with the cloned sequence from the 
canine polyp. DNA samples from 23 nasal polyp 
specimens resected from 12 patients without 
rhinosporidiosis were subjected to both Rhino- 
sporidium-specific and 15-globin PCR assays. 
These uninfected polyps failed to yield visible 
ampi icons after gel electrophoresis of the 
Rhi nosporidium-specific PCR reactions (data not 
shown). Flowever, all these samples yielded 
amplification products after (5-globin PCR, 
demonstrating that amplifiable DNA was 
present, without substantial PCR inhibition. 
These results confirmed that the presence of the 
putative R. seeberi 185 rDNA sequence 
correlated with the presence of disease 
(rhinosporidiosis). 



Figure 3. Agarose gel electrophoresis of Rhi nospo- 
ridium-specific PCR products. The specific amplifica- 
tion product is 377 bp. No amplification product is 
seen in the negative control samples consisting of 
water (reagent-only control), digestion buffer (DB), or 
lymph node tissue control (Tis Cnt). The human 
rhinosporidiosis samples (RST3) and the original 
canine nasal polyp show visible amplification 
products. 


Fluorescence in Situ Hybridization 

We sought to determi ne by usi ng F I SFI if the 
R. seeberi 185 rDNA sequence was linked to 
visible pathology in tissue. The Rhinosporidium 
185 rRN A probe but not the control probe bound 
to R. seeberi organisms in tissue, providing 
further evidencethat the putative R. seeberi 185 
rRN A sequence is present in R. seeberi organisms 
(Figure 4). To test the specificity of the 
hybridization, smears of the fungus C. albicans, 
cells infected with the Rosette agent, and tissue 
sections ofC. immitisweresubjectedtoFI5FI. No 
specific hybridization to these organisms was 
detected by usingtheRhinosporidium 185 rRNA 
probe compared with the control probe (data not 
shown). 

Electron Microscopy 

Transmission electron micrographs were 
taken of R. seebai trophocytes in a human nasal 
polyp. Multiple mitochondria were visualized, 
showing that the cristae had tubulovesicular 
morphology (Figure 5), unlike the flat cristae 
found in fungi. 

Discussion 

In the 1890s, first Mai bran and then 
5eeber (14) described an apparent sporozoan 
parasite in nasal polyps from patients living in 
Argentina. 5eeber's teacher, Wernicke, named 
the organism Cocci dium seeberi a after the 
protozoal subdivision Coccidia and his pupil, 
Guillermo 5eeber (2). In 1923, Ashworth 
described the life cycle of the organism, argued 
that it is a fungus, and proposed the name 
R. seeberi (15). 5incethen, the microbe has been 
considered a fungus by most microbiologists, 
although its taxonomy has been debated (1,2,16). 
Using a consensus PCR approach, we amplified a 
unique 185 rDNA sequence from a canine nasal 
polyp infected with R. seeberi. To prove that this 
unique 185 rDNA sequence came from R. seeberi, 
we sought to fulfi 1 1 sequence-based guidel i nes for 
microbial disease causation, since Koch's postu- 
lates cannot be fulfilled for uncultivated 
microbes (17). Using a Rhi nosporidium-specific 
PCR assay, we showed that the unique rDNA 
sequence present in the canine polyp was also 
present in three human polyp samples from 
patients with rhinosporidiosis. We showed the 
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Figure 4. An oligonucleotide probe complementary to a unique region of the Rhinosporidium seeberi 18S rRNA 
sequence localizes to visible organisms in human nasal tissue by using fluorescence in situ hybridization. 
Confocal micrographs at 200X magnification. The Cy-5-labeled Rhinosporidium probe (blue) hybridizes to 
spherical R. seeberi trophocytes (A) and a sporangium with endospores (B). A control probe consisting of the 
complement of the Rhinosporidium probe labeled with Cy-5 does not hybridize to the trophocytes (C) or a 
sporangium with endospores (D). I mages were collected in two wave-length channels: thefirst for Texas Red (red) 
or FITC (green) displays tissue architecture through autofluorescence, and the second channel for Cy-5 (blue 
pseudocolor) displays the probe signal. 
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Figure 5. Transmission electron micrograph of a 
mitochondrion from Rhinosporidium seeberi. The 
cristae of this mitochondrion (arrows) have 
tubulovesicular morphology. Magnification 195,000X. 

specificity of this association by demonstrating 
that this rDNA sequence was not present in 
control nasal polyps from 12 patients without 
rhinosporidiosis. We also used FI SH to link the 
putative R. seeberi 18S rRNA sequence to 
organisms visible in tissue. Although the 
Rhinosporidium probe localized to R. seeberi 
organismsintissueusingFI SH,theprobedid not 
localize to two fungal organisms or another 
member of the DRIPs clade (Rosette agent), 
suggesting some specificity to the hybridization. 

Sequence-based data also provide support for 
a causal relationship when a microbial genotype 
(e.g., phylogenetic placement) correctly predicts 
microbial phenotype and host response (17). In 
other words, the nature of the pathogen i nferred 
from phylogenetic analysis of its nucleic acid 
sequence should be consistent with the known 
biologic characteristics of closely related mi- 
crobes and with the nature of the disease. 
Therefore, the nearest phylogenetic neighbors to 
R. seeberi could be predicted to have similarities 
in morphology, tissue histology, and pathogen- 
esis. D. salmonis, the closest known relative to 
R. seeberi, has a large spherical structure 
containing endospore-like daughter cells (18). 
Flistology of infected hosts shows gill inflamma- 
tion and epithelial hyperplasia. The resemblance 
between R. seeberi and fish pathogens has been 
noted before. In 1960, Satyanarayana wrote in 
his review of 255 cases of rhinosporidiosis (1) that 
sinceR. seeberi hasa morphology similar tothose 


of some fish parasites, it "..may also be a parasite 
or saprophyte of fish and that man, equines, and 
cattleobtain the infection through water in which 
fish harbour! ngthe parasite live." 

A recent independent report based on 
amplification of 185 rDNA from two human 
rhinosporidiosis tissue samples also concludes 
that R. seeberi is a member of the DRI Ps clade of 
microbes (19). Our data support this conclusion 
and provide more evidence for a causal 
relationship. In addition, we describe a Rhino- 
sporidium-specific PCR assay that can be used 
for detecting this organism in clinical and 
environmental samples. Our 185 rDNA sequence 
differs from the sequence determined by these 
investigators (GenBank AF 118851) at a single 
position out of 1,699 common bases. We excluded 
from analysis sequence derived from our PCR 
primers, as ampi icons will contain the primer 
sequences regardless of the target sequence as 
long as there is partial annealing during PCR, 
leading to potent! ally spurious conclusions about 
sequence in these regions. The 185 rDNA 
sequence of R. seeberi determined by this group 
includes the primer sequences. 

Although we describe early trophocytes with 
mitochondrial cristae having vesicular ultra- 
structure, other investigators have found sporan- 
gia with mitochondrial cristae having a flat 
ultrastructure (19). The different mitochondrial 
morphologies observed may be due to differences 
in developmental stage of the organism or in 
methods of tissue preparation for electron 
microscopy. Nevertheless, another member of 
the DRIPs clade, Ichthyophonus hoferi, has 
vesicular mitochondrial cristae. Classic fungi 
(Eumycota) have flat mitochondrial cristae. 

Knowledge of the molecular phylogeny of 
R. seeberi is more than an exercise in taxonomy. 
For organisms such as R. seeberi, which are 
difficult to grow in the laboratory, phylogenetic 
analysis provides some insight into the charac- 
teristics of the organism that can be used to 
further our understanding of disease pathogen- 
esis and epidemiology, as well as to improve 
diagnosis and treatment. Knowingthat R. seebai 
is a member of the DRI Ps clade of microbes allows 
hypotheses to be generated about how it causes 
human disease by analogy, drawing on the 
knowledge and experience of the veterinary 
sciences. The separate but linked observations 
that rhinosporidiosis in humans is associated 
with exposure to water and that R. seeberi 
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belongs to a clade of aquatic parasites lead to a 
testable hypothesis: the natural hosts of 
R. seeberi are fish or other aquatic animals, and 
humans acquire infection when they come into 
contact with water containing these fish and 
their parasites. Investigators should therefore 
look for evidence of infection in fish in ponds and 
rivers in disease-endemic areas. From a public 
health perspective, the R. seeberi -specific PCR 
assay can be used to study environmental sources 
of infection (e.g., specific bodies of water) and 
may provide a means of preventing disease 
through identification of infected water. 

Conversely, knowing that R. seeberi is a 
member of the DRIPs clade may help us 
understand this important, distinct group of 
microbes that appear to form the deepest branch 
in the animal lineage. R. seeberi is a member of a 
newly recognized group of human and animal 
pathogens; the name I chthyosporea has been 
proposed for this expanding taxon of microbes 
(20). Little information is available about these 
organisms and how they cause disease. We hope 
that collaborations between researchers in 
human and animal medicine will correct this 
deficiency. 

Multiple antimicrobials, including antifun- 
gal agents, have been used in the treatment of 
rhinosporidiosis, based on the belief that 
R. seeberi is a fungus. However, no anti microbial 
agent is clearly effective. The medical treatment 
of rhinosporidiosis may be improved through 
screening anti parasitic drugs for an effect on 
disease in Dermocystidium-infected fish or 
infected cell lines. 

In conclusion, phylogenetic analysis of the 
R. seeberi 18S rRNA gene suggests that this 
culture-resistant organism is not a member of the 
Eumycota, but rather is the first known human 
pathogen from a novel clade of aquatic protistan 
parasites that form a branch in the evolutionary 
treeneartheanimal-fungal divergence. R.seeberi- 
specificPCR and FISH confirmthe association of 
this unique 18S rDNA sequence with the 
presence of rhinosporidiosis. This knowledge can 
be used to further our understanding of the 
natural reservoir of this organism and the risk 
factors, pathogenesis, and treatment of this 
disease. This discovery also expands our 
appreciation of the diversity among eukaryotic 
organisms that are pathogenic to humans and 
highlights the limitations of basing phylogenetic 
classification on morphology alone. 
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During 1997, Oklahoma City’s Hospital A reported penicillin-nonsusceptible 
Streptococcus pneumoniae in almost 67% of isolates. To confirm this finding, all 
Hospital A S. pneumoniae isolates from October 23, 1997, through February 19, 1998, 
were tested for antibiotic susceptibility and repeat-tested at two other hospital 
laboratories. Medical records of Hospital A patients with invasive S. pneumoniae 
infections during 1994 through 1997 were also reviewed. These data were compared 
with 1998 statewide sentinel hospital surveillance data for invasive S. pneumoniae. Of 
48 S. pneumoniae isolates from Hospital A during October 23, 1997, through February 
19, 1998, 31 (65%) were penicillin-nonsusceptible S. pneumoniae, and 23 (48%) were 
highly penicillin resistant. Similar prevalences were confirmed at the other hospital 
laboratories; however, significant interlaboratory differences were noted in the 
determination of third-generation cephalosporin susceptibility. During 1994 through 
1 997, a trend toward increasing penicillin nonsusceptibility (p <0.05) was noted among 
S. pneumoniae isolates from nursing home patients. During 1998, 85 (30%) of 282 
invasive isolates reported to the state surveillance system were penicillin- 
nonsusceptible S. pneumoniae; 33 (12%) were highly resistant. The increase in 
resistance observed is notable; the interlaboratory discrepancies are unexplained. To 
respond, a vaccination program was implemented at Hospital A, and vaccination efforts 
were initiated at nursing homes. 


Streptococcus pneumoniaeis a major cause of 
bacterial pneumonia and meningitis in the 
United States. The spread of penicillin- 
nonsusceptible S. pneumoniae (PNSP) has been 
well documented (1-8). Within a large city, 
resistance patterns can vary by region (5,8). 
Oklahoma City has historically had a relatively 
high but stable prevalence of PNSP among 
invasive isolates; penicillin nonsusceptibility 
was reported in 12.2% of invasive isolates during 
1984 (9). From J uly 1989 through J une 1990, 
7.6% of invasive isolates from central Oklahoma 
werepenicillin nonsusceptible, but high penicillin 

Address for correspondence: Ronald L. Moolenaar, Centers for 
Disease Control and Prevention, Mallstop D18, 1600 Clifton 
Rd„ Atlanta, GA 30333, USA; fax: 404-639-4504; e-mall: 
rl mSCgcdc.gov. 


resistance (1.4%) was beginning to emerge (10). 
In late 1997, higher than expected levels of PNSP 
were reported in northwest Oklahoma City at 
Hospital A, a 392-bed community hospital, 
primarily providing adult care. The hospital's 
inpatient census and the proportion of patients 
receiving Medicare had not increased since 1994; 
however, the microbiology laboratory had 
recently begun using the antimicrobial gradient 
strip method for measuring antibiotic 
susceptibility to penicillin and cephalosporins. 
Nonsusceptibility to penicillin among 
pneumococcal isolates exceeded 60%; 
approximately half were highly penicillin 
resistant. We initiated an investigation to 
determine whether the reported prevalence was 
accurate and if so, to explain it, determine a 
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possible trend, and comparethe prevalence with 
that of local hospitals. Findings were used to 
guide local treatment and prevention measures. 

The Study 

Prospective Laboratory Survey 

To confirm the accuracy of the preliminary 
antibiotic susceptibility results from Hospital A, 
we prospectively collected all S. pneumoniae 
isolates identified by Hospital A's microbiology 
laboratory from October 23, 1997, through 
February 19, 1998. An invasive isolate was 
defined as any positiveculturefor S. pneumoniae 
obtained from blood; cerebrospinal fluid; joint, 
pleural, or peritoneal fluid; or other normally 
sterile site. We confirmed the antibiotic 
suscepti bi I ity profi les of each of these isolates by 
retesting them at another community facility. 
Hospital B. To determine antibiotic susceptibil- 
ity, both hospitals used oxacillin disk screening 
followed by the antimicrobial gradient disk 
(E-test). The penicillin susceptibility of the 
invasive isolates was also confirmed using broth 
dilution at the laboratory of Children's Hospital 
of Oklahoma, which serves as the regional 
pediatric referral hospital in Oklahoma. 

For this report, the term penicillin 
nonsusceptible (PNSP)is used to describe any 
S. pneumoniae organism with reduced suscepti- 
bility to penicillin. We used the National 
Committee for Clinical Laboratory Standards' 
cut pointstoidentify anti biotic susceptibility and 
determine whether nonsusceptible organisms 
had intermediate or high resistance (for 
penicillin, < 0.06 pg/mL was considered 
susceptible; 0.10-1.00 pg/mL was intermediate, 
and >2.00 pg/mL was resistant) (11). 

Hospitals A and B used the E -test to evaluate 
the susceptibility of each isolate to third- 
generation cephalosporins, but Hospital A used 
cefotaxime as the test antibiotic whereas 
Hospital B used ceftriaxone. Both hospitals also 
tested isolates for susceptibility to trimethoprim 
sulfamethoxazole (TMP-S). 

The serotypes of 16 of the S. pneumoniae 
isolates (the first 13 invasive isolates collected 
and 3 additional noninvasive isolates) from 
Hospital A were determined by the Division of 
Bacterial and Mycotic Diseases Laboratory, 
Centers for Disease Control and Prevention 
(CDC). The susceptibility profiles and serotypes 
of these isolates were compared with those 


obtained from a 1996 nursing home outbreak of 
invasive multidrug-resistant S. pneumoniaein a 
geographically distant region of Oklahoma (12). 

Hospital A Retrospective Cohort Study 

TodescribetheepidemiologyofS. pneumoniae 
at Hospital A during 1994 through 1997, we 
retrospectively identified all patients with 
invasive isolates by reviewing records of the 
hospital's microbiology laboratory. After linking 
these with patient medical records, we extracted 
information on demographics, medical history, 
clinical course, and results of anti biotic suscepti- 
bility testing of the isolates for each patient. 
Trends over time were determined by the chi- 
square for trend test, and potential predictor 
factors for acqui ri ng nonsuscepti ble organi sms- 
such as hospitalization within the past year at 
Hospital A, nursing home residence, and various 
demographic factors— were evaluated (Epi-Info 
version 6.04b, CDC). 

1998 Oklahoma Sentinel Hospital Surveillance 

During this investigation, a sentinel hospital 
surveillance system was started for invasive S. 
pneumoniae in Oklahoma. We compared the 
prevalence of PNSP at Hospital A with data from 
the 26 sentinel hospitalsthat participated in this 
survei 1 1 ance system i n 1998. T en were among the 
acute-care hospitalsstudied in central Oklahoma 
during 1989 to 1990 (10); the others were 
scattered throughout the state. Some hospitals 
used bacterial broth dilution or antimicrobial 
gradient strips for penicillin susceptibility 
testing; others used bacterial disk diffusion or an 
antimicrobial panel. The sentinel hospitals 
accounted for approximately 58% of all medical 
and surgical hospital beds in the state. 

Results 

Prospective Laboratory Survey 

F rom October 23, 1997, through February 19, 
1998, Hospital A's microbiology laboratory 
identified S. pneumoniae isolates from 48 
patients. Of these, 17 (35%) were invasive: 2 were 
from cerebrospinal fluid, and 15 were from blood. 
The remaining 31 (65%) were noninvasive; 22 
isolates were from sputum and 9 from nose, 
sinus, tonsils, or bronchial washings. The median 
patient age was 60 years; only 5 (10%) were <18 
years of age. Of patients, 26 (54%) were male. 
T wenty (80%) of the 25 for whom race was known 
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were white, 4 (16%) were black, and I (4%) was 
Asian. Twenty-nine (60%) of the patients resided 
in Oklahoma County; the rest were from seven 
other Oklahoma counties. The isolates were 
tested at hospitals A and B for susceptibility to 
several different antibiotics. 

Hospital A reported that 31 (65%) of the 48 5. 
pneumoniae isolates were not susceptible to 
penicillin; 8 (17%)had intermediate resistance 
and 23 (48%) had high penicillin resistance. 
T wenty-three (48%) isolates were nonsusceptible 
to cefotaxime, and 13 (27%) were highly 
cefotaxime resistant. Of 47 isolates tested, 29 
(62%) were nonsusceptible to TMP-S; 23 (49%) 
were highly resistant (Table 1). All isolates were 
suscepti ble to van corny ci n. 

When these isolates were tested for penicillin 
susceptibility at Hospital B, a similar prevalence 
of PNSP was reported (62%), though fewer were 
highly resistant (29%). However, when testing 
for susceptibility to a third-generation cepha- 
losporin was performed at Hospital B, only 11 
(23%) of 48 isolates were reported as 
nonsusceptibletoceftriaxone, compared with the 
23 (48%) nonsusceptible to cefotaxime at 
Hospital B (Table 1). For all 15 isolates for which 


Table 1. Antibiotic-susceptibility test results® for 48 
isolates of Streptococcus pneumoniae from Hospital A, 
Oklahoma, October 23, 1997, through February 19, 
1998, and retest results from Hospital B 



Hospital A 
No. (%) 

Hospital B 
No. (%) 

Oxacillin disk 

Sensitive 

16/47 (34) 

17/46 (37) 

Resistant 

31/47 (66) 

29/46 (63) 

Penicillin 

Sensitive 

17/48 (35) 

18/48 (38) 

1 nter mediate 

8/48 (17) 

16/48 (33) 

Resistant 

23/48 (48) 

14/48 (29) 

Ceftriaxone 

Sensitive 

not performed 

37/48 (77) 

1 ntermediate 


7/48 (15) 

Resistant 


4/48 (8) 

Cefotaxime 

Sensitive 

25/48 (52) 

not performed 

1 ntermediate 

10/48 (21) 


Resistant 

13/48 (27) 


Trimethoprim- 

sulfamethoxazole 


Sensitive 

18/47 (38) 

16/47 (34) 

1 ntermediate 

6/47 (13) 

2/47 (4) 

Resistant 

23/47 (49) 

29/47 (62) 


®By E-test and oxacillin disc. 


cefotaxime susceptibility results at Hospital A 
differed from ceftri axone suscepti bi I ity results at 
Hospital B, the difference was in the direction of 
higherresistanceto cefotaxi me t h an to ceft r i axone 
(sign test, p <0.001). 

When penicillin susceptibilities of the 17 
invasive isolates tested with the E-test by both 
hospitalswerecompared with the susceptibilities 
determined at Children's Hospital of Oklahoma 
using broth dilution, the results were again 
similar. For 11 (65%) of the invasive isolates, all 
three test results were interpreted as the same. 
When discordance was noted, the differences 
were minor (i.e., susceptible was interpreted as 
intermediate [or vice versa], or intermediate was 
interpreted as resistant [or vice versa]). No 
major differences occurred in interpretation of 
penicillin susceptibility (susceptible to resistant 
or resistant to suscepti ble) (T able 2). 

Among the 16 serotyped isolates, 10 
serotypes were identified: 4, 6A, 9V (3 isolates), 
12F, 18C, 19A (2 isolates), 22F, 23F (4 isolates), 
33F, and 35B. Only two (12.5%)(6A and 35B) 
would not have been covered by pneumococcal 
vaccine. Among the 23F isolates, two were 
invasive, and two were not. One of the invasive 

Table 2. Penicillin-susceptibility test results from three 
laboratories for 17 invasive isolates of Streptococcus 
pneumoniae, Hospital A, Oklahoma, November 1997 
through February 1998 


Source 

E-test 

Broth dilution 

Hospital 

A 

Hospital 

B 

Children's Hos- 
pital of Oklahoma 

Blood 

1 

S 

S 

Blood 

R 

1 

R 

Blood 

R 

R 

R 

Blood 

R 

1 

1 

Blood 

S 

S 

S 

Blood 

1 

1 

1 

Blood 

R 

R 

1 

Blood 

S 

S 

s 

Blood 

R 

1 

1 

Blood 

S 

S 

s 

CSF 

S 

s 

s 

Blood 

S 

s 

s 

Blood 

1 

1 

1 

CSF 

R 

R 

1 

Blood 

R 

R 

R 

Blood 

R 

R 

R 

Blood 

S 

S 

S 

1 or R % 

65% 

59% 

59% 


I, having intermediate resistance: R, highly resistant; S, 
suscepti ble to the anti bi oti c. 

CSF = cerebrospinal fluid 
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23F isolates was highly sensitive to everything 
tested; the other three were highly resistant to 
several antibiotics and had similar susceptibility 
profiles to the 23F isolates identified in the 
nursing home outbreak months earlier (12). 

Hospital A Retrospective Cohort Study 

Review of the epidemiology of S. pneumoniae 
infection at Flospital A during 1994 through 1997 
revealed 71 case-patients with invasive infec- 
tions. Forty-three (61%) were female, 62 (87%) 
were white, and 7 (10%) were black. The median 
age was 64 (range: lto94), and 21 (30%) patients 
died. Twenty-three (32%) had been hospitalized 
at Flospital A during the previous 6 months. 
Nineteen (27%) were residents of a nursing 
home, 51 (72%) acquired infection in the 
community, and one (1%) acquired infection 
while hospitalized for another illness. Five (7%) 
had asplenia. Sixty-five (92%) isolates were from 
blood; four (6%) from the joint, pleural, or 
peritoneal fluid; and one (1%) each from 
cerebrospinal fluid and an aortic graft. Fifty-six 
(79%) i nfected patients resided in Oklahoma City 
or its vicinity. 

Fifty-eight (82%) patients with invasive 
S. pneumoniae infections during this period fit 
the Advisory Committee on Immunization 
Practices' (ACIP) criteria for pneumococcal 
vaccine (13). Of the 21 deaths, 20 (95%) would 
have met the criteria. Only one patient's record, 
however, showed receipt of pneumococcal 
vaccine. This patient had had a splenectomy and 
did not survive the infection. 

For the 71 case-patients, 27 (38%) isolates 
were PNSP, and 13 (18%) were highly resistant. 
Fifteen (21%) were nonsusceptible to cefotaxime, 
and one (1%) was highly resistant. All isolates 
tested were sensitive to vancomycin and 
clindamycin. Forty-two (59%) patients with 
invasive pneumococcal disease initially received 
ceftriaxone, and 3 (4%) received vancomycin; 20 
(28%) received vancomycin at some point during 
their hospital stay. Case-patients with PNSP 
infections were not more likely to die than 
patients with penicillin-susceptible infections. 

Over the 4-year period of the study, a trend 
toward increased nonsusceptibility to penicillin 
was noted at FI ospital A (1 of 9, 8 of 24, 7 of 20, 11 
of 18; p = 0.02). The same trend was noted for 
cefotaxime nonsusceptibility (0 of 9, 3 of 24, 5 of 
20, 9 of 18; p =0.008). A significant trend toward 
increased penicillin nonsusceptibility was noted 


in nursing home patients (Oof 2, 2 of 6, 4 of 7, 4 of 
4; p <0.05) but not in non-nursing home residents 
(2 of 7, 7 of 18, 4 of 13, 7 of 14; p >0.05). N ursi ng 
home residents were more likely than non- 
nursing home residents to have a nonsusceptible 
strain of S. pneumoniae, but this finding was not 
significant (10 of 18 vs. 20 of 53, risk ratio [RR] = 
1.5; 95% confidence interval [Cl] = 0.9 - 2.5; 
p > 0.05). Similarly, patients who had been 
hospitalized at Flospital A in the previous 6 months 
were more likely to have PNSP infections than 
those who had not been hospitalized there, but 
againthe association was not significant (18 of 33 
vs. 12 of 38, RR =1.7; 95%CI =1.0-3.0; p =0.09). 

1998 Oklahoma Sentinel Surveillance 

Surveillance data were available from 26 
sentinel hospital laboratories in Oklahoma 
during 1998, including Flospital A's. Seventeen 
laboratories used the E-test, four used broth 
dilution, three used disk diffusion (Kirby Bauer), 
and two used an antimicrobial panel (Microscan) 
for susceptibility testing. Of 282 invasive isolates 
tested, 197 (70%) were penicillin susceptible; 52 
(18%) had intermediate resistance, and 33 (12%) 
were highly penicillin resistant. Of the 26 
sentinel hospitals, 13 had at least 10 invasive 
isolates of S. pneumoniae during 1998. The 
prevalence of PNSP ranged from 10% to 45% in 
these 13 hospitals. 

Discussion 

The high prevalence of PNSP invasive 
isolates observed among vaccine-eligible, elderly 
adults from Flospital A in Oklahoma City is 
similar tothat recently reported asthehighest of 
a range of proportions from CDC's Emerging 
Infections Program's Active Bacterial Core 
Surveillance for 1997 (8). Penicillin 

nonsusceptibility varied amongOklahoma senti- 
nel hospitals, but overall it has increased 
markedly since studies were conducted in 
Oklahoma in the 1980s (9,10). 

The increase in penicillin nonsusceptibility 
at FI ospital A appears to be the result of a trend 
over the past 4 years. While we found no evidence 
of an increase in the number of patients with 
invasive S. pneumoniae coming from nursing 
homes, we did note an increasing prevalence of 
PNSP among such patients. 

The community in which Flospital A is 
located is fairly affluent. Resistance levels can 
vary within acity, and more affluent or suburban 
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communities may have higher rates of antibiotic 
resistance among S. pneumoniae isolates (5), 
perhaps attributable to more frequent use of 
antibiotics. Some resistant isolates were similar 
to those detected in a recent nursing home 
outbreak involving a single clone of a highly 
resistant organism (12). Contribution of an 
outbreak clone to the high prevalence observed 
here is possible; however, our evidence suggests 
that several resistant serotypes contributed to 
the observed increase. Although the reason for 
thePNSP increase at Hospital A is unclear, local 
antibiotic prescribing practices and one or more 
highly resistant clones circulating in the 
community or nursing homes may have been 
involved. 

This investigation had at least three 
limitations. First, thenumber of invasive isolates 
for the retrospective study was small, and 
specific risk factors for acquiring PNSP at 
Hospital A were not identified. Thus, no single 
factor explains the high rates observed. Other 
studies have suggested that age (< 6, or > 65), 
race, recent antibiotic use, socioeconomic status, 
and geographic factors may be associated with 
resistance (2,5). Second, because of the small 
number of isolates collected prospectively 
beginning in October 1997, we could not confirm 
or explain the significant difference in third- 
generation cephalosporin susceptibility reported 
by hospitals A and B. Although these facilities 
used two different test antibiotics (cefotaxime 
and ceftri axone) to determi ne suscepti bi I ity, the 
gen et i c mech a n i sm respon si bl e f or ceph a I ospor i n 
resistance would have been expected to have 
made the same organisms equally susceptible to 
either antibiotic. If confirmed, this finding is 
relevant in clinical settings where cefotaxime is 
used to test for susceptibility and ceftriaxone is 
used for treatment because S. pneumoniae may 
not have equal susceptibility to both antibiotics. 
Although we found reliability among the three 
laboratories for determining penicillin suscepti- 
bility, consistent with reports in the microbiology 
literature (14,15), the lack of interlaboratory 
reliability in differentiating between high and 
intermediate penicillin resistance was also 
unexpl ai ned. F i nal ly, wedid not assess whether 
treatment failure contributed to illness and 
death. Although case-patients with PNSP 
infections were not more likely to die than 
patients with penicillin-susceptible infections, 
the impact of underlying illnesses and the 


virulence of the infecting strain on the outcome is 
not known. 

For Hospital A and the state of Oklahoma, 
our findings have implications for disease 
treatment and prevention. Surveillance, in- 
creased use of the pneumococcal vaccine, and 
judicious use of antimicrobial drugs are 
important components of the effort to limit the 
spread of resistant S. pneumoniae (3, 13, 16). For 
determining community-specific prevalence of 
PNSP infections, use of hospital antibiograms 
has been shown to be comparable to active 
surveillance with centralized testing (17). 
However, our investigation suggested that when 
different third-generation cephalosporin test 
antibiotics are used and when determining high 
versus intermediate penicillin nonsusceptibility, 
significant differences can occur between hospi- 
tal laboratories in the determi nation of suscepti- 
bility. Such differences, whatever their explana- 
tion, have important clinical implications for the 
optimal use of antibiotics (e.g., vancomycin, 
ceftriaxone, and fluoroquinolones). Accurate 
laboratory information is needed at the local level 
to optimize use of antibiotics and minimize the 
development of antibiotic resistance. Accurate 
surveillance information is needed at the state 
level to compare regional prevalences. Further 
study with more isolates is needed to identify 
ways to i mprove the accuracy of thi s survei 1 1 ance 
system. 

This investigation documented that many 
cases of invasive pneumococcal disease could 
have been prevented by improved immunization 
practices: 82% of Hospital A case-patients during 
1994 to 1997 were eligible for vaccine by AC I P 
criteria, but only one had record of receiving it. 
Although hospital records may not have reflected 
true vaccination status, 1997 data from the 
Behavioral Risk Factor Surveillance System 
showed that 40% of Oklahoma residents > 65 
reported havi ng recei ved the pneumoni a vacci ne, 
which is less than the national goal of >60% (18). 
Targeting vaccination programs at nursing 
homes may be particularly effective. Oklahoma 
has recently required that all nursing homes 
offer pneumococcal vaccination totheir residents 
and provide documentation of vaccination status. 

I n addition. Hospital A instituted a hospital- 
based pneumococcal immunization program with 
standing orders to vaccinate inpatients aged 65 
years who are eligible for the vaccine based on 
ACIP criteria (13). Because nearly one out of 
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three of the patients with i nvasi ve pneumococcal 
disease admitted to Hospital A during 1994 to 
1997 had been hospital ized there i n the previous 
6months, an inpatient screeningand vaccination 
program could have prevented a substantial 
number of cases. Hospital -based programs have 
been shown to be effective in vaccinating high- 
risk adults in other settings (19-21), and the use 
of standing orders to vaccinate eligible adults 
with pneumococcal vaccine has been recom- 
mended by the Task Force on Community 
Preventive Services (22). Adult vaccination 
programs in nontraditional settings (e.g., pharma- 
cies, churches, and the workplace) might further 
raise vaccination coverage in this area (23). 

Finally, efforts to promote judicious use of 
antibiotics are needed to minimize the spread of 
antibiotic-resistant S. pneumoniae Several 
studies have demonstrated that antimicrobial 
drug use is associated with resistance in S. 
pneumoniae (24-26). A combination of interven- 
tions involving education of both physicians and 
patients has been successful in reducing 
antibiotic use (27) and would likely reduce the 
trend of increasing resistance in Oklahoma. 
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Rickettsia mongoiotimonae: 

A Rare Pathogen in France 
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andtHopital de la Conception, Marseille, France 


We report a second case of laboratory-confirmed infection caused by Rickettsia 
mongoiotimonae in Marseille, France. This rickettsiosis may represent a new clinical 
entity; moreover, its geographic distribution may be broader than previously 
documented. This pathogen should be systematically considered in the differential 
diagnosis of atypical rickettsioses, especially rashless fevers with lymphangitis and 
lymphadenopathy, in southern France and perhaps elsewhere. 


Rickettsia mongoiotimonae (1) was first 
obtained as isolate HA-91 from Hyalomma 
asiaticum, a tick collected in Inner Mongolia in 
1991 (2). The fi rst evi dence of R . mongol oti monae 
pathogenicity in humans was obtained in our 
laboratory in 1996 in a woman who had an 
atypical tick-transmitted disease. This second 
case of R. mongoiotimonae infection was 
unexpected, especially because the symptoms in 
the second case differed from those of the first 
case and of other rickettsioses cases. 

On May 19, 1998, a 49-year-old HIV-positive 
man was admitted to a hospital in Marseille, 
F ranee, with fever and myalgias for 5 days. The 
patient, who had been HIV positive since 1988, 
was being treated with a combination of 
zidovudine, lamivudine, and nevirapine. With 
this therapy, he had been asymptomatic for 5 
years. His CD4-I- lymphocyte count was 503/pL, 
and his plasma viral level was 5,916 copies per ml 
(3.77 log). He was a gardener who lived in a rural 
area 40 km from Marseille and had contact with 
animals, including cats, rabbits, squirrels, and 
birds. On May 8, 1998, he worked in his garden 
but did not notice any tick bites. On examination 
the day of admission, he had an eschar (tache 
noi re) surrounded by an inflammatory halo on his 
left leg. He had a fever of 39°C, chills, headache, 
and myalgias. A lymphangitis expanded from the 
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inoculation eschar to the left groin, where a 
painful satellite lymphadenopathy was observed 
(Figure). Clinical examination was otherwise 
normal. Laboratory results were unremarkable. 



Figure. Lymphangitis expanding from the inocula- 
tion eschar on the left leg to an enlarged, painful 
lymph node on the left groin of a patient with 
Rickettsia mongoiotimonae infection. 
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An acute-phase serum sample drawn that day 
was sent to our laboratory, along with a biopsy of 
the eschar. Although the patient had no rash, 
Mediterranean spotted fever was suspected. 
After 15 days of treatment with 200 mg per day of 
doxycycline, the symptoms resolved. A second 
serum sample was collected on May 26. 

Microimmunofluorescence testing was per- 
formed (3). Although the first serologic test was 
negative, the second serum sample had low titers 
to R. mongolotimonae and R. conorii (1:16 for 
immunoglobulin [Ig] G and 1:16 for IgM). The 
eschar sample was spread onto human embryon- 
iclungfi broblasts bythecentrifugation shell-vial 
technique (4). After 7 days incubation at 32°C, a 
Gimenez staining of methanol -fixed human 
embryonic lung fibroblasts cells showed Rickett- 
sia-like bacilli. DNA was extracted from the 
ground eschar biopsy and from 200 pL of shell- 
vial supernatant, by using a QIAmp Tissue kit 
(QIAGEN GmbH, H ilden, Germany) accordingto 
the manufacturer's instructions. These extracts 
were used as templates with primers complemen- 
tary to portions of the coding sequences of the 
rOmpA, citrate synthase, and 165 rRNA encoding 
genes in three polymerase chain reaction (PCR) 
assays (5). The basesequences of the obtained PCR 
products were determined (5). All sequences 
shared 100% similarity with R. mongolotimonae 
when compared with GenBank data. 

Eight new rickettsioses, three of which were 
described in Europe, have been identified in the 
past 13 years, including] apanesespotted fever due 
toR.japonica; Astrakhan fever due to Astrakhan 
fever rickettsia; Flinder's Island spotted fever 
due to R. honei; cat flea typhus due to R. felis; 
African tick-bite fever due to R. africae; and the 
proposed rickettsioses d ue to R. mongolotimonae 
(1), R. slovaca (6), and R. helvetica (7). R. conorii, 
responsible for Mediterranean spotted fever, has 
long been considered the only tick-transmitted 
European rickettsia in the Mediterranean Sea 
area. However, the results of this and previous 
studies (1,6) indicatethat both R. mongoloti monae 
and R. slovaca are agents of human rickettsioses 
in France. 

I ntheonly reported caseofR. mongolotimonae 
infection in France, migratory birds carrying 
ticks were hypothesized to be responsible for 
transmittingthe infection, asR. mongolotimonae 
had been isolated only from ticks from China (1). 
On the basis of this hypothesis, the initial case 


was thought to be a unique imported infection. 
However, the occurrence of a second case in the 
same area suggests that the disease may be more 
common than expected around Marseille. R. 
mongolotimonae was initially isolated from 
Hyalomma asiaticum ticks in China, but its 
reservoir in our area is unknown. However, the 
fact that the patient worked in a garden near 
cypress trees where numerous birds nested 
suggests that migratory birds should be 
considered a possible source. We were unable to 
detect this bacterium by using ompA PCR- 
amplification and culture in 540 Rhipicephalus 
ticks from the Marseille area; these ticks are 
prevalent in southern France (8,9). An active 
search for the vector of R. mongolotimonae, 
which may be a rare tick species in Southern 
France, should be conducted. 

The symptoms of both patients occurred in 
spring, which is unusual for Mediterranean 
spotted fever (1). The symptoms of the first 
patient were I ike those of Mediterranean spotted 
fever, except that the rash was discrete. The 
second patient had norash but had lymphangitis, 
expanding from the inoculation eschar to the 
groin, and a painful lymphadenopathy. These 
differences in clinical symptoms may be 
explained by the HIV-positive status of this 
patient. Two cases are inadequate to describe a 
disease; however, our preliminary data suggest 
that R. mongolotimonae infections may have a 
discrete rash or no rash, along with lymphangitis 
and satellite lymphadenopathy. 

This second case of R. mongolotimonae 
infection is important for several reasons. First, 
this rickettsiosis may represent a new clinical 
entity. Second, thegeographicdistribution of the 
illness may be broader than previously docu- 
mented. R. mongoloti monaeshould be systemati- 
cally considered in the differential diagnosis of 
atypical rickettsioses in southern France, 
especiallyfevers without rash but with lymphan- 
gitis and lymphadenopathy, in southern France 
and perhaps elsewhere. 

Dr. Fournier is a physician in the French reference 
center for the diagnosis and study of rickettsial diseases. 
Hisresearch interests include the clinical and epidemio- 
logic features of the rickettsioses. 
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We assessed the societal costs and benefits of a subtype-specific surveillance 
system for identifying outbreak-associated Escherichia co// 01 57;H7 infections. Using 
data from Colorado, we estimated that if it averted five cases annually, the system would 
recover all its costs. 


Escherichia coli 0157:H7 infections pose a 
serious public health threat (1-4). Surveillance, 
rapid reporting of cases, and prompt epidemio- 
logic investigations are essential elements of 
timely public health response (2,5). Surveillance 
that uses molecular subtyping methods has at 
least two advantages over traditional surveil- 
lance systems (6). First, it is sensitive enough to 
identify outbreaks not detected by traditional 
surveillance or can detect them earlier. Second, it 
is specific enough to differentiate sporadic cases 
from outbreak-related cases and distinguish 
between single and multiple outbreaks. 

A subtype-specific surveillance system con- 
sists of 1) mandatory submission of E. coli 
0157:H7 isolates for subtyping; 2) a centrally 
located I aboratory equi pped to perform subty pi ng 
by pulsed-fieldgel electrophoresis; 3)activelinks 
between local and state health officials; and 4) 
epi demi ol ogi c capacity to i nvesti gate the possi bi I - 
ity of an outbreak once identical strains are 
identified. 

I n August 1997, the Colorado Department of 
Public Health and Environment, using subtype- 
specific surveillance, identified an outbreak 
associated with eating hamburgers from beef 
processed in a plant in Nebraska and distributed 
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nationally. After the outbreak was traced to the 
contaminated beef, the company recalled 25 
million pounds of ground beef, the largest meat 
recall recorded (7). 

We used cost-benefit analysis to assess the 
economic feasibility (from a societal perspective) 
of using a system similar to the one in Colorado 
for identifying E. coli 0157:H7 outbreaks. 

The Study 

A system is cost-beneficial if the discounted 
benefits it generates are at least as great as the 
discounted costs of installing and operating the 
system. The life span of the subtype-specific 
surveillance system in Colorado is 5 years, 
yielding benefits over the lifetime of people 
affected by it. Data on the costs of the system 
were obtained from the Colorado Department of 
Public Health and Environment (Table 1). The 
system was not set up only to subtype for E . col i 
0157:H7 but also to identify outbreaks of other 
organisms (e.g.. Salmonella typhi); only E. coli- 
related costs were considered here. The 
sensitivity of the results was examined with all 
the costs of the system attributed to E. coli 
0157:H 7 subtyping. 

I n estimating the costs of outbreak investiga- 
tions, we assumed that, as a result of the system, 
two epidemiologic investigations would be 
carried out each year, with an average cost of 
$9,600 per outbreak (Table 1). The costs 
associated with recalling any outbreak-related 
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Table 1. Costs of installing and operating the subtype-specific surveillance system, Colorado, 1996 


Labor and equipment costs 

Total costs 

Escherichia coli-related costs® 

Equipment 

$40,000 

$16,000 

Laboratory scientist (per year)'’ 

$10,000 

$4,000 

Analyzing the isolates (per year)'’ 

$12,000 

$12,000 

Investigating an outbreak''® 

$9,600 

$9,600 

Present value of outbreak costs (in 5 years)' 

$90,568 

$90,568 

Annual operating costs^ 

$41,200 

$35,200 


^Frorn the proportion of E. coll isolates among the total number of isolates expected to be subtyped each year, we extrapolated 
that 40% of the equipment and labor costs were E . coli-related. 

'The salary and fringe benefits of a full-time laboratory analyst. 

"^Analyzing 300 isolates at a cost of $40 per isolate. 

This cost included, but was not limited to, the value of time (15 days) spent investigating an outbreak, answering telephone 
calls, conducting meetings, improving and transferring pulsed-field gel electrophoresis image files to various groups, creating 
databases, requesting information, responding to media calls, and handling legal issues. Weassumed that, as a result of the 
system, two outbreaks would be investigated each year (6). 

^he costs of additional labor and the epidemiologic investigation of an outbreak were estimated at $5,000 and $4,600, 
respectively. 

^At a discount rate of 3%. 

9Laboratory scientist ($10,000) -F analyzing the isolates ($12,000) -F invest! gating two outbreaks (2 x $9,600 =$19,200) 


prcxJuct were not included. Data (e.g., percentage 
of contaminated beef) that would allow us to 
attribute economic value to the amount of the 
product recalled were not available. In the 
sensitivity analysis, costs were increased by 
100% to account for such missing data. 

The benefits of the surveillance system are 
the economic savings accrued from E. coli 
0157:H 7 cases averted. Determi ni ng the number 
of cases averted as a result of using the system is 
difficult. One way of determining this number is 
by estimating the attack rate and multiplying 
that number by the amount of beef recalled (8). 

However, for the outbreak in Colorado, data 
for estimating specific attack rates were lacking. 
Instead, we estimated two threshold numbers of 
cases that must be averted for the costs to be 
equal to the benefits of the system. The first 
threshold number was calculated by assuming 
the system averts a constant number of cases 
every year. The second number was calculated 
under the assumption that the system averts 
only a given number of cases in the first year and 
no cases in subsequent years. If the estimated 
threshold is below a reasonable number, the 
system is cost beneficial. A reasonable number is 
calculated by consulting the literature and expert 
opinion. 

The average cost of an E. coli 0157:H7 
infection was estimated by using an infection 
outcome tree (4) (Figure). A person infected with 
E. coli 0157:H7 can be in only one disease 
severity category (Figure; Table 2). The far-left 


branch of the tree is designated as severity 
category no. 9. The probability is 0.2% (10% 
hospitalization x 50% hemolytic uremic syn- 
drome [H US] x 4% death) that an infected person 
will be hospitalized for hemorrhagiccolitis, come 
down with HUS, and die after 1 year. F rom the 
time of infection until the time of death, the 
societal costs for this patient are $991,221 
(medical costs $39,204 + productivity losses 
$3,041 -Flost lifetime earnings $948,976). 

Data on the probability of being in any of 
these categories were obtained from Roberts et 
al. (3). The economic costs associated with each 
category were based on the methods and 
assumptions of Buzby et al. (4), with modifica- 
tions (Table 2). Productivity losses were 
estimated by multiplying the average wage in 
1996 by the number of days missed from work. 
The average wage rate was estimated by using 
the average daily earnings of a nonagricultural 
nonsupervisory employee, assuming that fringe 
benefits are 39% of total wages or salaries and a 
labor participation rate of 84% (4). Weestimated 
the costs of a death by using onlythelost lifetime 
earnings, as estimated by Haddix et al. (9) and 
updated by the rates of change in wages. Because 
we did not assess pain and suffering from the 
disease or loss of human life, our estimates of 
benefits should be considered conservatively low. 

All costs and benefits were adjusted to 1996 
dollars according to the consumer price index or 
its components (various issues of the Statistical 
Abstract). All future costs and benefits were 
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Figure. Escherichia coli 0157:H7 infection outcome tree. Severity categories (l)-(9) are described in Table 2. 


Table 2. Assumptions about disease severity following an Escherichia coli 0157:1-17 infection^ 

Severity 

category Assumptions 

No. 1 Patient does not seek medical care, recovers, and misses 2 days of work 

No. 2 Patient seeks medical care for hemorrhagic colitis, has one laboratory test, recovers, and 

misses 4 days of work 

No. 3 Patient is hospitalized for hemorrhagic colitis for 6.5 days and recovers after missing 14 days 

of work 

No. 4 Patient is hospitalized for hemorrhagic colitis for 6.5 days, misses 14 days of work, and dies 

i n the fi rst year 

No. 5 Patient is hospitalized for acute HUS'^ for 5 days in ICU'^ and 10 days in a regular room, and 

recovers after missing 32 days of work 

No. 6 Patient is hospitalized for acute HUS'^ for 5 days in ICU'^ and 10 days in a regular room, 

requires dialysis for 12 days, and recovers after missing 32 days of work 
No. 7 Patient is hospitalized for hemorrhagic colitis; comes down with chronic HUS'’; may require 

dialysis, transplants, or drug therapy; cannot work for an extended period; and recovers 
No. 8 Patient is hospitalized for hemorrhagic colitis; comes down with chronic HUS'’; may require 

dialysis, transplants, or drug therapy; cannot work for an extended period; and dies 
No. 9 Patient is hospitalized for acute HUS'’ for 5 days in ICU'’ and 10 days in a regular room and 

dies after missing 32 days of work 

^Adapted fromBuzby et al. (4). A patient is defined as a person infected with E. coli 0157:H7 whohasat least a gastrointestinal 
1 1 1 ness for more than 1 day. 

‘’HUS, hemolytic uremic syndrome. 

SCU, intensive care unit. 
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discounted at 3%. Other rates were used in the 
sensitivity analysis. 

The discounted average cost of an E. coli 
0157:H7 infection was $7,788 (Table 3). The 
main component of the cost of a case was the 
expected cost of sequelae and death. The 
undiscounted cost of a case was $15,927. The 
discounted cost of installing and operating the 
surveillance system over a period of 5 years was 
$182,042 (Table3). I ncluded in this category were 
the costs of i nvestigati ng outbreaks ($90,568 i n 5 
years), of thesubtypingequipment ($16,000), and 
of analyzing the isolates ($60,000 in 5 years) 
(T able 1). At a 3% discount, five cases per year (or 
14 cases i n thefi rst year only) must be averted for 
the costs of the system to be equal to its benefits 
(Table 3). Without discounting, the threshold 
number dropped to 2.4 cases per year (T able 3). 

Sensitivity Analysis 

T he est i mated costs of a case were sensi t i ve to t he 
estimates of the probability of death after 
infection. If the probability of death is raised to 
2.3% (4), the cost of a case increases to $25,997, 
andthethreshold number of cases averted for the 
system to be economically feasible decreases to 
1.5 per year for 5 years, or 4.3 cases in the first 
year and none in the following years (Table 3). 

If all the costs of subtyping (including 
subtyping for other organisms) were included in 
the analysis, the system would recover its costs 


after averting 6.4 infections annually, or 20.9 
casesin thefi rst year only, with no cases detected 
in subsequent years. If the costs of the system 
doubled or the benefits of a case averted 
decreased by 50%, the threshold number would 
i ncrease to 9.9 cases per year, or 28.4 cases i n the 
first year only (Table 3). Doubling the number of 
outbreaks or considering only direct medical 
costs would raise the threshold numbers to 7.4 
and 11.4 cases per year, respectively. 

Conclusions 

I f 15 cases were averted bytherecall ofthe25 
million pounds of potentially contaminated beef, 
the Colorado system would have recovered all 
costs for the 5 years of start-up and operation by 
detecting a single outbreak (Table 3). In 
comparison, the outbreak-related 1993 recall of 
255,000 regular (0.1-lb) hamburgers in Washing- 
ton State was estimated to have prevented 800 
cases (8). 

The discounted average cost of an E. coli 
0157:H7 infection of $7,788 (Table 3) was a 
relatively conservative estimate compared 
with that of $38,000 (in 1995 dollars) by Buzby 
et al. (4). The major differences are the 
probability of death and theeconomic value of life 
used in the estimation (11). 

If other benefits of the system (e.g., obviating 
the need to investigate sporadic cases) are 
included, the system becomes even more cost 


Table 3. Discounted costs of an Escherichia co// 01 57:1-17 infection, discounted costs of the surveillance system, and 

threshold number of cases, 1996 

Discount rate® 


Components affecting costs 

3% 

0% 

5% 

Discounted average cost of an E. coli 0157:1-17 infection 

$7,788 

$15,927 

$5,847 

Discounted costs of installing and operating the system 
Baseline (best estimate) 

$182,042 

$192,000 

$176,018 

Cases that need to be averted every year for 5 years'’ 

5.0 

2.4 

6.6 

Cases that need to be averted in the first year alone'’ 
One-way sensitivity analysis 

1 ncreasing labor and equipment costs from $20,000 to $50,000 

14.2 

7.2 

18.5 

Cases that need to be averted every year for 5 years'’ 

6.4 

3.1 

8.6 

Cases that need to be averted in the first year alone” 
Decreasing the costs of an infection from $7,788 to $3,894 

20.9 

10.6 

27.2 

Cases that need to be averted every year for 5 years'’ 

9.9 

4.8 

13.2 

Cases that need to be averted in the first year alone” 
1 ncreasing the probability of death from 0.4% to 2.3% 

28.4 

14.5 

36.9 

Cases that need to be averted every year for 5 years'’ 

1.5 

0.5 

2.6 

Cases that need to be averted in the first year alone” 

4.3 

1.5 

7.3 


®The most frequently assumed discount rate is 5%. However, 3% is the recommended rate. No discounting is suggested for 
testing the sensitivity of the results (10). 

•Threshold number of cases averted every year for 5 years above which the system is cost-beneficial . 

Threshold number of cases averted in the first year above which the system is cost-beneficial, assuming the system does not 
avert any cases in subsequent years and continues to incur costs. 
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beneficial. Unproductive extensive traceback 
investigations of sporadic E. coli 0157:H7 
infections have been conducted (12). Investi- 
gating such sporadic cases can be very costly 
(Table 1), and a subtype-specific system can 
reduce such costs. 

According to the National Electronic Tele- 
communications System for Surveillance, 90 
cases of E. coli 0157:H7 were reported in 
Colorado in 1998(13), an annual incidence rate 
of 2.3 per 100,000 population. In comparison, 
the national incidence rate calculated from 
these data was 1.2 per 100,000 population 
(3,161 cases). 

This study was limited by lack of data that 
would have enabled us to estimate attack rates 
from the outbreak, cases averted by the meat 
recall, and the benefit tosociety (money saved) by 
establishing the system. Despite its limitations, 
this study has important implications for public 
health policy. From a societal perspective, a 
surveillance system does not need to prevent a 
large number of cases to yield return on the 
resources invested in it. 
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Dengue Epidemic in Beiem, 

Para, Brazii, 1996-97 
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I nstitutoEvandro Chagas, Belem, Para, Brazil 


We describe clinical and epidemiologic findings during the first epidemic of dengue 
fever in Belem, Para State, Brazil, in 1 996-97. Of 40,237 serum samples, 1 7,440 (43%) 
were positive for dengue by virus isolation or serologic testing. No hemorrhagic cases 
or deaths were reported. 


Dengue fever (DF) and dengue hemorrhagic 
fever (DHF) are caused by infection with one of 
the four serotypes of dengue virus (DEN-1, 
DEN -2, DEN -3, and DEN -4), transmitted by 
Aedes aegypti mosquitoes. In Brazil, DF 
epidemics reported i n the 1980s and 1990s i nvolved 
more than a million cases. However, only 671 
DH F cases were diagnosed, with 26 deaths (1). 

Ae aegypti was reintroduced in Para State in 
1992, and the first dengue cases were reported in 
1995 in the southeast region (Redengao and 
Rondon do Para). I n October 1996, eight cases of 
febriledenguelike illness were reported in Belem 
(population 1,300,000), a city in the Brazilian 
Amazon region at the confluence of the Amazon 
River andtheAtlanticOcean. I n early November, 
DEN-1 virus was isolated and identified (2-4). 
DEN-2 virus was identified in October 1997, and 
si nee then, both serotypes have been responsible 
for illness in Belem. This was the first time 
dengue virus transmission occurred in Belem 
during the last 70 years and the third time the 
disease occurred in theBrazilian Amazon region. 
Previous outbreaks had been reported in 1981-82 
in Boa Vista, Roraima (5), and in 1991 in 
Araguaina, Tocantins State (6). 

We describe cases of denguelike illness 
diagnosed at I nstituto Evandro Chagas. A case of 
dengue was defined as illness with the following 
symptoms: acute onset of high fever, headache. 


Address for correspondence: Pedro Fernando da Costa 
Vasconcelos, I nstituto Evandro Chagas, Av. Almirante 
Barroso492, 66090-000, Belem, Para, Brazil; fax: 55-91-226- 
1284; e-mail: pedrovasconceloscgiec.pa. gov.br. 


myalgia, arthralgia, dizziness, and other symp- 
toms and signs suggestive of denguel i ke i 1 1 ness i n 
a patient with a positive IgM by I gM -capture 
enzyme-linked immunosorbent assay (MAC 
ELI SA), virus isolation or serologic conversion in 
paired serum samples, and an increase of at least 
fourfold in titer in the convalescent-phase serum 
sample (7-9). 

The Study 

From January to December 1997, 40,237 
serum samples were drawn from febri le patients 
in Belem, 20,038(49.8%)of whom were male. The 
patients were all residents of the municipalities 
of Belem (31,506 samples) and Ananindeua 
(8,731 samples). Most were ambulatory patients 
seen in the Arbovirus Department, Evandro 
Chagas I nstitute; some patients were referred by 
public health centers, private physicians, and 
clinics. At the Institute, all patients were 
examined clinically and had blood samples 
drawn. A questionnaire was administered that 
included information about clinical symptoms 
and signs and demographic data. Patients who 
had been ill <10 days and whose serologic tests 
were negative for dengue were requested to 
provide a second blood sample 7 to 14 days later. 
For patients who had at least one hemorrhage 
and either dehydration or hemoconcentration, a 
leukogram, platelet count, and hematocrit were 
performed. 

Serum samples were tested for dengue 
anti bodi es by M AC E L I SA (4) and hemaggi uti na- 
tion-inhibition test (2). The antigens for both 
tests were prepared by using infected mouse 
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brain extracted by the sucrose acetone method. 
The criteria used for establishing primary and 
secondary infection were those recommended by 
the World Health Organization (7,8). 

To isolate dengue virus, 0.1-ml aliquots of 
whole blood from patients with clinical 
symptoms lasting <5 days and negative 
serologic results were injected into cultures of 
C6/36 cells (10). The cultures were visually 
examined daily, and cells were tested on days 7 
and 14 by immunofluorescence (3). 

I sol ated virus strains were initially screened 
by di rect i mmunofi uorescenceagai nst a fl avi vi rus 
standard hyperimmune fluid prepared at the 
I nstitute. Strai ns that reacted were identified to 
serotype by using an indirect i mmunofi uores- 
cencetest with monoclonal anti bodies against the 
four dengue viruses provided by the Centers for 
Disease Control and Prevention. 

The epidemic distribution was accompanied 
early in 1997 by a seasonal increase in rainfall 
typical of the Brazilian Amazon. However, the 
highest dengue positivity rates were reported in 
the dry months of September to December. 
Coincidentally, DEN -2 virus was isolated in 
October . At f i rst , cases were reported on I y i n B el em 
and Ananindeua, but in December 1997, at least 15 
other municipalities reported transmission 
cycles involving either DEN-1 or DEN -2 or both. 

Of all sera collected, 17,525 (43.5%) were 
positive by serologic testing or DEN virus 
isolation: 13,805 (78.8%) from Belem and 3,720 
(21.2%) cases from Ananindeua (Figure). Both 
primary and secondary serologic responses were 
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Figure. Monthly distribution of rain precipitation 
(line) and dengue fever cases as seen at Institute 
Evandro Chagas, Para, Brazil. Source for rain 
precipitation. Institute Nacional de Meteorologia, 
Brazil. 

found. Among the positive samples, 9,469 
(54.25%) were among female and 7,971 (45.75%) 
among male patients (p <0.0001); 45.95% of the 
patients were 25 to 44 years of age (Table 1). 
Paired serum samples were obtained from 3,558 
patients, 2,997 (84.2%) of whom had serologic 
conversions. 

No DH F cases were reported, in spite of the 
fact that since October 1997 DEN -2 had been 
isolated from 24 patients with previous DEN-1 
infection (Table2). Clinical symptoms were more 
severe in older children and adults than in young 
children, but no differences were observed 
between the clinical symptoms of patients with 
DEN-1 and DEN-2 infection. 


Table 1. Distribution of patients with serologic tests positive for dengue 


Age 


Male 



Female 


Total 

No. 

Pos 

% 

No. 

Pos 

% 

Pos 

% 

<4 

186 

60 

0.76 

325 

76 

0.8 

136 

0.78 

5-9 

308 

89 

1.12 

425 

186 

1.96 

275 

1.58 

10-14 

1,033 

444 

5.56 

816 

356 

3.76 

800 

4.59 

15-24 

4,080 

1,596 

20 

3,765 

1,643 

17.35 

3,239 

18.56 

25-34 

5,083 

2,184 

27.37 

4,375 

2,127 

22.47 

4,311 

24.7 

35-44 

3,919 

1,684 

21.1 

4,638 

2,023 

21.37 

3,707 

21.25 

45-54 

3,319 

1,016 

12.73 

3,572 

1,558 

16.46 

2,574 

14.75 

>55 

2,110 

898 

11.36 

2,283 

1,500 

15.83 

2,398 

13.79 

Total 

20,038 

7,971 

100 

20,199 

9,469 

100 

17,440 

100 
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Table 2. Strains of dengue virus isolated from patients 
from Belem and Ananindeua, Institute Evandro Chagas, 
October 1996 to August 1998 



DEN-1 

DEN -2 

Total 

1996 

31 

- 

31 

1997 

751 

48 

799 

1998^ 

58 

60 

118 

Total 

840 

108 

948 


^Through August. 


Conclusions 

Duri ng 1953 and 1954, Causey andThei ler (11) 
used seroneutralization in mice to survey several 
municipalities of the Amazon region; they found 
9.8% and 2.2% of serum samples positive for 
dengue virus serotypes DEN-1 and DEN -2, 
respectively, among residents >50 years of age. 
These results suggest that dengue viruses were 
circulating in Belem and other municipalities of 
the Amazon Valley early in the 20th century. 

The dengue epidemic in Belem was unusual 
in that a "lag phase" (12) lasted for at least 4 
months before extensive transmission began. I n 
spite of laboratory diagnosis of cases from the 
beginning of the outbreak and the reporting of 
these results to health authorities, control 
measures were unsuccessful. Consequently, 
explosive transmission began in March 1997 and 
is still occurring (1998-99). 

I n the last 4 years, the Brazilian M inistry of 
Health (13) has reported an increase in dengue 
cases from approximately 56,000 in 1994 to 
>530,000 in 1998. High indexes of Ae aegypti 
infestation exist in all important urban centers of 
Brazil, reflecting a poor health education 
program. The fact that more cases occurred in 
femalethan in malepatients(p <0.0001) suggests 
that women are at increased risk for dengue 
infection because of high peridomestic exposure. 

Denguetransmission has been reported i n 22 
of the 27 Brazilian states, and the mosquito 
vector is present in all states; therefore, the 
situation in Brazil may be rapidly approaching 
hyperendemicity, with the cocirculation of two 
serotypes (DEN-1 since 1986 and DEN -2 since 
1990). The risk for DHF will increase if a new 
serotype (DEN-3 or DEN-4) is introduced. 
E ndemicity is the most constant factor associated 
with the evolution of epidemic DHF in a 


geographic area (14). Although few cases have 
been reported until recently, DHF may become 
an important cause of hospitalization and death 
in the Americas, including Brazil. 
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To assess whether the Mycobacterium fubercu/os/s Beijing genotype is emerging in 
Vietnam, we analyzed 563 isolates from new cases by spoligotyping and examined the 
association between the genotype and age, resistance, and BCG vaccination status. 
Three hundred one (54%) patients were infected with Beijing genotype strains. The 
genotype was associated with younger age (and hence with active transmission) and 
with isoniazid and streptomycin resistance, but not with BCG vaccination. 


A high degree of diversity of Mycobacterium 
tuberculosis has been shown with restriction 
fragment length polymorphism (RFLP) typing 
using 1 56110 as a probe, particularly in countries 
like the Netherlands, where many tuberculosis 
(TB) cases occur among immigrants (1). 
However, in the Beijing region of China, a 
particular genotype was found in >80%oftheTB 
patients and was thus designated the Beijing 
genotype (2). I n other parts of China and in Asian 
countries such asMongolia, Thailand, and Korea, 
40% to 50% of the tested M . tu bercu I osi s i sol ates 
represented this genotype (2). Although Beijing 
genotype strai ns carry a large number of 1 56110 
insertion elements, the 156110 RFLP patterns 
are highly similar (2). Moreover, the 
spoligopatterns of Beijing genotype strains are 
identical and distinct from those of other 
M. tuberculosis strains (1), which suggests that 
strains of the Beijing genotype have emerged 
recently from a single ancestor. 

Reasons for the predominance of a narrow 
range of genotypes may include limited contact 
with other populations or a selectiveadvantageof 
certain strains due to reduced sensitivity to 


Address for correspondence: D. van Soolingen, Diagnostic 
Laboratory for Infectious Diseases and Perinatal Screening, 
National I nstituteof Public Health andtheEnvironment, P.O. 
Box 1, 3720 BA Bilthoven, The Netherlands: fax: 31-30- 
2744418; e-mail: d.van.soolingen(3rivm.nl. 


vaccine-induced immunity. For instance, wide- 
scale application of vaccines against whooping 
cough (Bordetella pertussis) has led to shifts in 
thepopulationsof dreu I ating pathogens (3). BCG 
vaccination, which has been applied widely in 
China since the early 1950s, may have led to a 
similar shift in the population of M. tuberculosis. 
Alternatively, a selective advantage may be 
provided by reduced sensitivity toanti-TB drugs. 
The Beijing genotype was associated with recent 
transmission of drug-resistant strains in Cuba, 
Germany, Russia, and Estonia (4-7). The largest 
known epidemic of multidrug-resistant TB in 
North America was caused by the "W" strain, a 
variant of the Beijing genotype (8). 

Vietnam is one of 22 countries in which 80% 
of the world's new TB cases occur (9). Among 
these 22 countries, Vietnam has one of the most 
successful directly observed therapy short- 
course programs, with a cure rate of approxi- 
mately 90% and a case-detection rate estimated 
at 67% in 1996 and at >70% si nee then (National 
Tuberculosis Programme, unpub. data). BCG 
coverage has been high (>80%) during the past 
decade, and the level of primary multidrug 
resistance was recently estimated at 2.3% (9,10). 
We investigated the spread of M. tuberculosis 
Bei j i ng genotype strai ns i n Vietnam and whether 
the spread is associated with BCG vaccination 
status or drug resistance. 
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I n total, 822 isolates were obtained from TB 
patients whose age and BCG status were known. 
Of these, 563 had newly diagnosed disease (Table 
1). The isolates were collected in 1998 and the 
first quarter of 1999 at the Tuberculosis and 
Lung Diseases Centre in Ho Chi M inh City and 
the National Institute of Tuberculosis and 
Respiratory Diseases in Hanoi, respectively. The 
isolates were analyzed by spoligotyping (11). In 
spoligotyping, the genomic direct repeat (DR) 
region of M. tuberculosis complex bacteria is 
amplified by polymerase chain reaction (PCR) 
and the presence of 43 spacer sequences between 
theDRs is examined in a reversed line-blot assay 
(2,11). In previous studies, this method has 
proven highly reliable for distinguishing Beijing 
genotype strains (2). Among the 563 
spoligopatterns analyzed, two predominant 
genotypes were recognized (Figure), of which the 
Beijing type was the most frequent (n = 301; 
53%). Thesecond most frequent genotypewas not 
found in the database of spoligopatterns of 2,500 
M. tuberculosis isolates from countries all over 
theworld, held at the National I nstituteof Public 
Health and the Environment in Bilthoven and 
designated Vietnam genotype (n = 152; 27%) 
(Figure). The remaining 110 isolates exhibited 18 
different spoligopatterns. 

The Bei j i ng genotype was strongl y associ ated 
with younger age (X ^trend’ P <0-001), but not with 


BCG status, after the data were adjusted for age 
(Table 1). As levels of drug resistance varied 
between Hanoi and Ho Chi Minh City and the 
number of samples from Hanoi was small, the 
association of the Beijing genotype and drug 
resistance was restricted to Ho Chi Minh City 
(Table 2). Drug resistance was found more 
commonly in theBeijingthan in other genotypes 
(Table 2). The association with drug resistance 
was significant for isoniazid (OR 1.7; 95% Cl 1.1- 
2.6) and streptomycin (odds ratio [OR] 3.1; 95% 
confidence interval [Cl ] 2. 0-4.6) resistance. 

TB occurs partly as primary disease 
(typically defined as occurring within 5 years of 
infection) and partly as endogenous reactivation 
or exogenous reinfection (occurring >5 years 


If ■ * mmitie 

jo JO 30 iO 


Figure. Representative spoligotype patterns of the 
Beijing (A) and the Vietnam (B) genotypes. Numbers 
indicate the spacer oligonucleotide sequences, present 
on the reversed line blot, which are derived from 
reference Mycobacterium tuberculosis strain H37Rv 
and M. bovis BCG vaccine strain P3. 


Table 1. Beijing genetype of Mycobacterium tuberculosis in 563 new tuberculosis cases, by age and BCG status, in 
Hanoi and Ho Chi Minh City, Vietnam, 1998 



No. 

Beijing genotype 

OR3 (95% C 

1^) 

No. 

(%) 

Crude 

Adjusted‘S 

Province 






Hanoi 

64 

37 

(58) 

1 (p>0.05) 


HoChi Minh City 

499 

264 

(53) 

0.8 (0.5-1.4) 


Age group 






<25 

76 

54 

(71) 

1 (p<0.001) 

1 

25-34 

173 

102 

(59) 

0.6 (0.3-1.0) 

0.6 (0.3-1.1) 

35-44 

164 

83 

(51) 

0.4 (0.2-0.7) 

0.4 (0.2-0.8) 

45-54 

74 

31 

(42) 

0.3 (0.1-0.6) 

0.3 (0.2-0.6) 

55-64 

39 

16 

(41) 

0.3 (0.1-0.6) 

0.3 (0.1-0.7) 

65-F 

37 

15 

(41) 

0.3 (0.1-0.6) 

0.3 (0.1-0.7) 

BCG scar 






Yes 

285 

167 

(59) 

l(p<0.05) 

1 

No 

278 

134 

(48) 

0.7 (0.5-0.9) 

0.9 (0.6-1.3) 

Total 

563 

301 

(54) 




®OR, odds ratio. 

‘^Cl, confidence interval. 

‘^Adjusted for age and BCG vaccination. 
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Table 2. Risk factors for drug resistance in 499 new tuberculosis cases in Ho Chi Minh City 



No. 


Resistance (%) 


0RM95%CI‘’) 

INH 

SM 

RIF 

EMB 

MDR 

INH 

SM 

Genotype 









Beijing 

264 

28 

42 

3 

3 

3 

1.7 (1.1-2.6) 

3.1 (2.0-4.6) 

Other 

235 

19 

19 

2 

1 

2 

1 (p<0.05) 

1 (p<0.001) 

Age group 









<25 

66 

17 

36 

3 

0 

3 

1 (p>0.05) 

1 (p>0.05) 

25-34 

165 

27 

32 

3 

2 

3 

1.8 (0.9-3.8) 

0.8 (0.4-1.5) 

35-44 

147 

25 

33 

1 

1 

1 

1.7 (0.8-3.5) 

0.8 (0.5-1.6) 

45-54 

59 

19 

25 

3 

5 

3 

1.1 (0.5-2.9) 

0.6 (0.3-1.3) 

55-64 

30 

23 

23 

3 

3 

3 

1.5 (0.5-4.4) 

0.5 (0.2-1.4) 

65-r 

32 

31 

31 

0 

0 

0 

2.3 (0.8-6.1) 

0.8 (0.3-2.0) 

BCG scar 









Yes 

259 

22 

31 

3 

2 

3 

1 (p>0.05) 

1 (p>0.05) 

No 

402 

73 

22 

2 

2 

2 

1.3 (0.9-2.0) 

1.1 (0.7-1.6) 

Total 

499 

24 

31 

2 

2 

2 




®OR, cxdds ratio. 

‘’Cl, confidence interval. 

INH, isoniazide; SM, streptomycin: RIF, rifampin; EMB, ethambutol; MDR, multidrug resistant 


after infection). With increasing age, a decreas- 
ing proportion of cases is due to primary TB. 
Thus, the association of the Beijing genotype and 
young age suggests a recent spread of the Beijing 
genotype in Vietnam. The study of Qian et al., in 
which spoligotyping was performed on paraffin- 
embedded material, indicated that the Beijing 
genotype was presumably already prevalent in 
the Beijing region 30 to 40 years ago (12). 

I n Vietnam, the Beijing genotype occurs more 
commonly in those with a BCG scar than in those 
without it. However, this is likely to represent a 
cohort effect of BCG vaccination, rather than 
reduced sensitivity to vaccine-induced immunity 
of Bei j i ng genotype strai ns. Because of i ncreasi ng 
BCG vaccination coverage in Vietnam over the 
past 2 decades, young people are more I i kel y to be 
vaccinated than older people. Within age groups, 
occurrence of the Beijing genotype is not 
associated with BCG vaccination status. 

Its striking association with anti-TB drug 
resistance may explain the Beijing genotype's 
predominance in recently infected patients. Anti- 
TB drugs are widely used in Vietnam, and anti- 
TB drug resistancewould thus provide a selective 
advantage. In vitro experiments should deter- 
mine whether Beijing genotype strains have an 
increased intrinsic resistance to anti-TB drugs or 
an enhanced capacity to gain resistance against 
these drugs. 


Dr. Anh is head of the Department of Bacteriology, 
National Institute of Hygiene and Epidemiology, Hanoi, 
Vietnam. Areas of interest include molecular epidemiol- 
ogy of mycobacteri al i nfecti on and epi demi ol ogi c su rvei I - 
lance of bacterial meningitis in children. 
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Bartonella species were isolated from 49% of 128 cattle from California and 
Oklahoma, 90% of 42 mule deer from California, and 15% of 100 elk from California and 
Oregon. Isolates from all 63 cattle, 14 deer, and 1 elk had the same polymerase chain 
reaction/restriction fragment length polymorphism profiles. Our findings indicate 
potential for inter- and intraspecies transmission among ruminants, as well as risk that 
these Bartonella spp. could act as zoonotic agents. 


Bartonella species have been identified as 
important zoonotic agents (1,2). Cats are the 
main reservoir of Bartonella henselae, the agent 
that causes cat scratch disease in humans (1). 
Long-term bacteremia in cats and flea transmis- 
sion from cat tocat, as confirmed by experimental 
infection, support a vectorbornetransmission (3). 
Some human cases of cat scratch disease were not 
associated with any known exposure to cats (4), 
suggesting that other animal species may serve 
as reservoirs of Bartonella. Recently, new 
B a rtonel I a speci es h ave been isolatedfromawide 
range of mammals, including rodents (5-10), 
lagomorphs (11), carnivores (12-14), and cervids 
(14,15). Similarly, 90% of 42 mule deer 
(Odocoileus hemionus) from California were 
bacteremic with Bartonella isolates that were 
similartoisolatesfrom roedeer in France(15) by 
polymerase chain reaction/restriction fragment 
length polymorphism (PCR/RFLP) of the 16S 

Address for correspondence: BrunoB. Chomel, Department of 
Population Health and Reproduction, School of Veterinary 
Medicine, University of California, Davis, CA 95616, USA; 
fax: 530-752-2377; e-mail: bbchomelcgucdavis.edu. 


rRNA and citrate synthase genes (14). Modes of 
transmission in these ruminants need to be 
established. Tick transmission has been suspect- 
ed but not yet proven for dogs infected with 
B. vinsonii subsp. berkhoffii (16). Since fleas are 
less likely than ticks to infest cattle (17), ticks 
may playan important roleinthetransmission of 
Bartonella species from wild ruminants. 

Our objectives were to determine if elk 
(Cervuselaphus), bighorn sheep (Ovis canadensis), 
and domesti c cattle (Bos taurus) are i nfected with 
Bartonella and to determine the molecular 
relationships between Bartonella isolated from 
cattle and wi Id rumi nants. We performed a cross- 
sectional study to compare the prevalence of 
Bartonella infection in a beef cattle herd in the 
California Sierra Nevada foothills and a dairy 
herd from the California Central Valley. 

The Study 

In February 1997, 42 samples from free- 
ranging mule deer were obtained from the Round 
Valley population. Mono and Inyo counties, 
California. I n November 1997, 84 samples were 


^An earlier version of this paper was presented at the Second International Conference on Emerging Zoonoses, Strasbourg, 
France, November 5-9, 1998. 
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collected from bighorn sheep herds in Califor- 
nia and New Mexico. During J anuary and 
February 1998, 100 blood samples were 
collected from elk in California and Oregon. 
One hundred twenty-eight cattle samples were 
collected: 12 from Oklahoma beef cattle in April 
1998 and 116 from two California herds from 
May to J uly 1998. Fifty-three samples were 
collected from a >4, 000-head beef cattle herd in 
the Sierra Nevada foothills and 63 samples 
from a >1, 500-head dairy herd in the Central 
Valley. These 116 cattle were all > 2 years of 
age. Blood samples collected into lysis- 
centrifugation tubes were plated within 48 
hours. Blood samples collected into EDTA 
tubes were frozen at -70° until plated. Wildlife 
and domestic herds were selected on the basis 
of ongoing surveys by the California and 
Oregon Departments of Fish and Game and 
researchers at the Universities of California 
and Oklahoma.^ 

Blood samples were cultured on heart 
infusion agar containing 5% rabbit blood and 
incubated in 5%C02 at 35°C for at least 4 weeks 
(18). Gram staining and biochemical tests were 
performed on representative isolates, which were 
defined as isolates with a unique PCR/RFLP 
profile for each of the three ruminant species. 
Nine representative isolates were identified, 
including one cattle strain (pattern I), five deer 
strains (patterns I, II, IV, V, and VI ), and three 
elk strains (patterns I, II, and III). Standard 
methods were used to test for various preformed 
enzymes and carbohydrate use. Preformed 
bacterial enzyme activity was tested by M icroscan 
Rapid Anaerobe Panel (Dade I nternational Inc., 
West Sacramento, CA) (19). 

An approximately 400-bp fragment of the 
citrate synthase gene was amplified as described 
(20). The amplified product was digested with 
TaqI and Flhal and Msel restriction endonu- 
cleases and visualized by gel electrophoresis. 
Banding patterns werecompared with B. henselae 
(strain U-4; University of California, Davis, CA). 

Cellular fatty acid composition was analyzed 
for representative cattle, deer, and elk isolates. 
Isolates were grown on rabbit blood agar at 35°C 
for 5 days. Fatty acid methyl ester derivatives 
were separated on a Flewlett-Packard series II 
5890 gas chromatograph. 


The PCR products used for DNA sequencing 
were purified with Microcon centrifugal filter 
devices (Millipore Corp., Bedford, MA) and 
sequenced with a fluorescent-based automated 
sequencing system. Primer BhCS.1137n (5'- 
AATGCAAAAAGAACAGTAAACA-3') (20) was 
used for partial sequencing of the 400-bp product 
of the citrate synthase gene. N i ne representative 
strains from ruminants and one B. henselae 
strain (strain U-4, University of California, 
Davis) were sequenced. The GAP program of 
GCG software (Wisconsin Sequence Analysis 
Package, Genetics Computer Group, version 10) 
was used for alignments and comparisons of 
sequences, based on the 276 bp of the citrate 
synthase gene. 

Using E pi I nfo version 6.03, we performed a 
chi-square test to assess association between 
prevalence of bacteremia of Bartonella infection 
and herd location. The Bartonella infection 
prevalence ratio (PR) was calculated to show the 
proportionate increase of infection prevalence 
due to difference in herd location. 

Results 

Bartond la spp. were isolated from 5 (42%) of 
12 Oklahoma cattle, 58 (50%) of 116 California 
cattle, 38 (90%) of 42 California mule deer, 15 
(15%) of 100 elk, and noneof 84 bighorn sheep. I n 
the California beef cattle herd, 25 (96%) of 26 
bulls and 22 (81%) of 27 cows were Bartond la 
bacteremic; in the dairy herd, 11 (17%) of 63 cows 
were bacteremic. Bartond I a bacteremia preva- 
lence in the Sierra Nevada foothills beef cattle 
herd was therefore significantly higher than in 
the Central Valley dairy cattle herd (PR = 5.1; 
95% confidence interval [Cl] = 2. 9-8.8). Preva- 
lence of Bartondia bacteremic cows in the 
foothills herd was also significantly higher (81% 
vs. 17%) than in the Central Valley dairy cattle 
herd (PR = 4.7; 95% Cl = 2.7-8.2). For elk, 
bacteremia prevalence differed significantly (p = 
0.0002) between California (0 of 47) and Oregon 
(15 [28%] of 53). No Bartond I a-bacteremic elk 
were found in the two California herds, but 11 
(38%) of 29 elk from southwestern Oregon and 4 
(17%) of 24 elk from northwestern Oregon were 
bacteremic. 

The organisms isolated were short, slender 
gram-negative rods. By measuring preformed 


^Collalion sites for bighorn sheep werethe Peninsuiar Rangesin Caiifornia and the San Francisco River, Turkey Creek, and 
Red Rock in New Mexico. For eik, coiiection sites werethe San Luis Nationai Wiidiife Refuge in MercedCounty and theT upman 
T uie Eik State Reserve in Kern County (Caiifornia): the Roseburg, Drain, and Demet herds, Dougias County (southwestern 
Oregon): and thej eweii Wiidiife Area, Ciatsop County (northwestern Oregon). 
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enzymes (Rapid Anaerobe Panel), the tested 
strains were found to be biochemically inert 
except for the production of peptidases, charac- 
teristicof the Bartonella profile (10077640). 

Several strain profiles were observed by 
PCR/RFLP of the citrate synthase gene, using 
TaqI and Hhal and Msel endonucleases for deer 
(five profiles) and elk (three profiles) isolates 
(Figure). Conversely, all 63 cattle isolates had the 
same PCR/RFLP profile (Figure) with the same 
restriction enzymes. Overall, six different PCR/ 
RFLP profiles were obtained from Bartonella 
isolated from cattle, deer, and elk. Bartonella 
isolated from cattle (63 of 63 tested; lanes 2, 12, 
and 22), mule deer (14 of 38 tested; lanes 3, 13, 
and 23), and an elk from southwestern Oregon (1 
of 11 tested; lanes 10, 20, and 30) yielded the 
samePCR/RFLP profile (pattern I ) withthethree 
enzymes used. A second profile (pattern 1 1 ) was 
obtained for Bartonella isolated from elk 
captured in northwestern Oregon (4 of 4 tested; 
lanes 8, 18, and 28) and from mule deer (5 of 38 
tested; lanes 4, 14, and 24). A third profile 
(pattern III) was obtained for 10 of the 11 
Bartonella isolated from elk captured in 
southwestern Oregon (lanes 9, 19, and 29). The 
other three profiles (patterns IV, V, and VI ) were 
obtained for Bartonella isolated from mule deer 
([pattern IV: 12 of 38 tested; lanes 6, 16, and 26]; 


[pattern V: 5 of 38 tested; lanesS, 15, and 25]; and 
[pattern VI : 2 of 38 tested; lanes 7, 17, and 27]). 

The cellular fatty acid composition was 
characteristic of the Bartonella genus for all 
isolates. The main fatty acids observed for the 
cattle, deer, and elk strains were octadecanoic acid 
(Ci8;i, 45%-66%), octadecanoic acid (Ci8;o, 12%- 
23%), and hexadecanoic acid (Ci6;o, 13%-20%). 

After pairwise comparisons, the partial 
sequencing analysis (276 bp) of the citrate 
synthase gene for the nine representative 
ruminant strains showed a high percentage of 
DNA similarity, from 93.12% to 100% (Table 1). 
The strains cattle- 1, deer-1, and elk-1 belonging 
tothePCR/RFLP pattern I had 95.65% to 99.64% 
DNA similarity. The strains deer-2 and elk-2 
with PCR/RFLP pattern II had 100% DNA 
similarity. The strain deer-1 with PCR/RFLP 
pattern I was closely related (98.91% DNA identity) 
tothestrain deer-2 with PCR/RFLP pattern II. For 
strains deer-4 and deer-5, corresponding to PCR/ 
RFLP patterns IV and V (similar digestion profiles 
with Flhal and Msel endonucleases and different 
profiles from TaqI endonuclease), a 98.55% DNA 
similarity was observed. Partial sequence 
analysis (276 bp) of the citrate synthase gene 
showed that all strains from ruminants were 
closely related to B. weissii, a Bartonella species 
isolated from domestic cats (Table 2). 



Figure. Polymerase chain reaction/restriction fragment length polymorphism of the citrate synthase gene of 
isolates from cattle, deer, and elk, with TaqI, Hhal, and Msel endonucleases. Lanes 1 and 32, standard 100-bp 
molecular ladder; lanes 2, 12, and 22, cattle isolate; lanes 3 to 7, 13 to 17, and 23 to 27, deer isolates; lanes 8 to 
10, 18 to 20, and 28 to 30, elk isolates; lanes 11, 21, and 31, B. henselae strain. 
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Table 1 . DNA similarity values and GenBank accession numbers based on 276 bp of the citrate synthase gene of the 
nine representative ruminant strains 


Organism/ 
accession no. 




% Similarity bv strain 




Cattle-1 

Deer-1 

Deer-2 

Deer -4 

Deer-5 

Deer-6 

Elk-1 

Elk-2 

Elk-3 

Cattle-1 AF 228768 

100.00 

95.65 

96.01 

94.57 

94.57 

94.57 

99.64 

96.01 

94.57 

Deer-1 A F 228769 

- 

100.00 

98.91 

93.84 

93.84 

93.12 

95.65 

98.91 

93.84 

Deer-2 AF 228771 

- 

- 

100.00 

94.20 

94.20 

93.48 

96.01 

100.00 

94.20 

Deer-4 A F 228774 

- 

- 

- 

100.00 

98.55 

94.93 

94.57 

94.20 

94.57 

Deer-5 A F 228775 

- 

- 

- 

- 

100.00 

94.20 

94.57 

94.20 

94.57 

Deer-6 A F 228776 

- 

- 

- 

- 

- 

100.00 

94.20 

93.48 

96.01 

Elk-1 AF 228770 

- 

- 

- 

- 

- 

- 

100.00 

96.01 

94.57 

Elk-2 AF 228772 

- 

- 

- 

- 

- 

- 

- 

100.00 

94.20 

Elk-3 AF 228773 

- 

- 

- 

- 

- 

- 

- 

- 

100.00 


T able 2. DNA similarity values based on 276 bp of the citrate synthase gene of the nine ruminant strains compared with 
those of the Bartonella strains in GenBank 


% Similarity by strain 


Organism/accession no. 

Cattle-l 

Deer-1 

Deer-2 

Deer-4 

Deer-5 

Deer-6 

Elk-1 

Elk-2 

Elk-3 

B. bacilliformis U28076 

86.59 

87.68 

87.32 

84.78 

85.51 

84.78 

86.59 

87.32 

87.68 

B. grahamii Z70016 

90.22 

90.22 

90.58 

91.67 

90.22 

90.58 

90.58 

90.58 

89.49 

B. taylorii Z70013 

88.41 

87.32 

86.96 

87.68 

87.68 

87.68 

88.04 

86.96 

88.04 

B. tribocorumAJ 005494 

89.86 

89.13 

89.49 

90.58 

89.13 

88.41 

89.49 

89.49 

88.04 

B. doshiaeZ70017 

88.41 

86.59 

86.96 

86.96 

86.23 

87.68 

88.04 

86.96 

85.87 

B. vinsonii subsp. vinsonii U 28074 

88.69 

89.05 

87.96 

88.69 

88.69 

89.42 

88.32 

87.96 

87.96 

B. vinsonii subsp. berkhoffii U28075 

89.86 

89.49 

89.13 

87.68 

87.68 

88.41 

89.49 

89.13 

86.96 

B. vinsonii subsp. arupensis U 77057 

90.22 

89.13 

89.13 

90.94 

90.94 

90.22 

89.86 

89.13 

88.77 

B. weissii AF 071190 

99.64 

95.65 

96.01 

94.57 

94.57 

94.20 

100.00 

96.01 

94.57 

B. clarridgeiae U84386 

90.58 

89.49 

89.86 

88.77 

88.04 

88.77 

90.22 

89.86 

89.49 

B. henselae strain U-4 

90.58 

88.41 

89.49 

88.77 

88.77 

87.32 

90.22 

89.49 

87.68 

B. henselae strain Houston-1 L 38987 

90.58 

88.41 

89.49 

88.77 

88.77 

87.32 

90.22 

89.49 

87.68 

B. koehleraeAF 176091 

89.13 

88.41 

88.77 

89.49 

88.77 

87.32 

88.77 

88.77 

87.32 

B. quintana Z70014 

90.22 

88.04 

88.41 

87.68 

86.96 

88.41 

89.86 

88.41 

88.77 

B. elizabethaeZ70009 

88.41 

88.04 

88.41 

90.22 

88.77 

89.49 

88.04 

88.41 

88.04 

Strain C7-rat Z70020 

88.41 

88.04 

88.41 

90.22 

88.77 

89.49 

88.04 

88.41 

88.04 

Strain C5-rat Z70018 

88.77 

88.77 

89.13 

90.22 

88.77 

87.32 

88.77 

89.13 

87.68 

Strain C4-phyZ70019 

87.32 

86.23 

86.96 

86.96 

86.96 

86.59 

86.96 

86.96 

85.15 

Strain Cl-phyZ70022 

86.59 

85.51 

86.23 

86.23 

86.23 

85.87 

86.23 

86.23 

84.42 

Strain R-phy2 Z70011 

87.32 

86.23 

86.96 

86.96 

86.96 

86.59 

86.96 

86.96 

85.15 

Strain R-phyl Z70010 

88.04 

87.68 

88.04 

87.32 

87.32 

87.32 

87.68 

88.04 

85.87 

Strain N40Z70012 

90.22 

88.77 

89.13 

88.77 

87.32 

87.32 

89.86 

89.13 

86.96 

Strain A1 U 84372 

88.77 

87.68 

88.77 

88.04 

87.32 

88.04 

88.41 

88.77 

86.23 

Strain A2 U 84373 

88.41 

87.32 

88.41 

87.68 

87.68 

87.68 

88.04 

88.41 

86.23 

Strain A3 U 84374 

88.77 

88.04 

88.77 

88.04 

88.04 

88.04 

88.41 

88.77 

86.23 

Strain B1 U 84375 

88.49 

89.13 

88.77 

88.77 

88.77 

89.49 

89.13 

88.77 

88.04 

Strain B2 U 84376 

89.86 

89.49 

89.13 

88.41 

88.41 

89.13 

89.49 

89.13 

87.68 

Strain Cl U 84377 

88.77 

89.13 

88.77 

87.68 

86.96 

88.41 

88.41 

88.77 

86.59 

Strain C2 U 84378 

88.77 

89.13 

88.77 

87.68 

86.96 

88.41 

88.41 

88.77 

86.59 

Strain D1 U84379 

89.86 

88.77 

88.77 

90.58 

90.58 

89.86 

89.49 

88.77 

88.41 

Strain D2 U 84380 

90.22 

89.13 

89.13 

90.94 

90.94 

90.22 

89.86 

89.13 

88.41 

Strain D3 U 84381 

90.58 

89.86 

89.49 

90.58 

90.58 

90.58 

90.22 

89.49 

88.77 

Strain D4 U 84382 

90.22 

89.13 

89.13 

90.94 

90.94 

90.22 

89.86 

89.13 

88.77 

Strain D5 U 84383 

89.49 

89.13 

88.41 

90.22 

90.22 

89.49 

89.13 

88.41 

88.04 

Strain D6 U 84384 

90.58 

89.86 

89.49 

90.58 

90.58 

90.58 

90.22 

89.49 

88.77 

Strain D7 U 84385 

90.22 

89.13 

89.13 

90.94 

90.94 

90.22 

89.86 

89.13 

88.41 
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Conclusion 

This is the first published report of isolation 
of Bartonella spp. from free-ranging wild 
ruminants and domestic ruminants in North 
America. Our results suggest that deer, elk, and 
domestic cattle are possible reservoirs of 
Bartonella spp. Selected bighorn sheep popula- 
tions from California and New Mexico appeared 
to be free of Bartonella. The first report of 
infection of cattle with a Bartonella organism was 
made in 1934 by Donatien and Lestoquard, who 
proposed the name B. bovis or Haemobartonella 
bovis (21). In 1942, Lotze and Yiengst also 
described Bartonel la-1 ike structures in American 
cattle (22); however, their identifications of 
B artonel I a-l i ke structures were based only on the 
morphologic aspects of these organisms in red 
blood cells also infected with Theileria or 
Anaplasma, two well-known tickborne infections. 

Partial sequencing analysis of the citrate 
synthase gene of the ruminant strains showed 
that they were all closely related to each other 
and to a feline strain, B. weissii. Further studies 
by DNA-DNA hybridization may determine if 
these strains are specific to ruminants but 
closely related to B. weissii, or if they are in 
fact B. weissii. If the ruminant strains are 
identical toB. weissii, the high prevalence (89%) 
of Bartonella bacteremia observed in beef cattle 
may indicate that ruminants are the main 
reservoirs of B. weissii, which is not commonly 
isolated from cats. 

The prevalence of Bartonella bacteremia was 
high in beef cattleand muledeer, possi bly related 
to exposure to potential vectors. Since fleas are 
rarely observed on cattle and tick infestation is 
common in both cattle and deer, ticks are a 
possible source of infection for ruminants (17). 
Furthermore, Bartonella DNA has recently been 
demonstrated in a high percentage of ticks 
infesting roedeer in Europe (23,24). The herd of 
beef cattle from the Sierra Nevada foothills, 
wheretick infestation is common, has permanent 
access to open pastures. I n contrast, the dairy 
cattle herd from the Central Valley has little or 
no access to pastures and tick infestations are 
not commonly observed (R. Bon Du rant, pers. 
comm.). Therefore, geographic differences in 
the prevalence of Bartonella infection in 
California cattle herds warrant further investi- 
gation for possible tick transmission of 
Bartonella spp. among these animals. 


PCR/RFLP analysis of the citrate synthase 
gene has been widely used for identification of 
Bartonella organisms tothespecies level (25-27). 
We identified one PCR/RFLP profile for all the 
cattle isolates, but several profiles for deer and 
elk. This diversity by geographic location is of 
epidemiologic interest and warrants further 
investigation. Only one elk from southwestern 
Oregon had a strain with a similar PCR/RFLP 
profiletothat of domestic cattle, suggesting that 
wild ruminants could be infected with Bartonella 
speciesthat arenot commonly shared with cattle. 

Our findings also suggest that transmission 
of Bartonella may occur among cattle and 
wildlife, especially mule deer, which are more 
abundant in the western USA than elk and are 
morelikelytobesympatricwith cattle. Collection 
and analysis of ticks on wild animals and cattle 
and from the environment will be necessary to 
determine if ticks can be infected with Bartonella 
species. Whether Bartonella isolated from these 
ruminants are human pathogens is still unclear. 
The recent report of a cattle rancher who was 
infected with a new B. vinsonii subspecies (28) 
warrants further investigation to establish if 
these Bartonella species could be zoonotic and 
whether humans could potentially be infected by 
tick bites during work or recreation. 

Dr. Chang is pursuing his Ph.D. in epidemiology at 
the University of California, Davis, under the direction 
of Bruno B. Chomel. His research interests include epl- 
demi ol ogy of zoonoses, especi al I y the mol ecul ar epi demi - 
ology of Bartonella infections and potential vectors for 
Bartonella spp. transmission. 

Dr. Chang's research was funded by a grant from 
theCenter for Companion Animal Health, University of 
California, Davis, California, USA. 
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Carbapenem-Resistant Pseudomonas 
aeruginosa with Acquired biavm 
Metaiio-f}-Lactamase Determinants, itaiy 

To the Editor: Acquired metalloB-lactamase 
determinants in Pseudomonas aeruginosa and 
other major bacterial pathogens are of concern 
for development of anti microbial drug resistance. 
The carbapenemase and extended-spectrum 
cephalospori nase activity of metal lo- 15-lactamases, 
as well as their resistance to 15-lactamase 
inhibitors, may severely limit the antimicrobial 
agents active against bacterial strains that 
produce such enzymes (1,2). Antimicrobial 
chemotherapy may become ineffective against 
P. aeruginosa strains with a multidrug-resistant 
phenotype that have acquired a metal lo-l5- 
I actamase determi nant. 

We recently described a new acquired 
metallo-l5-lactamase determinant, biavm, in a 
carbapenem-resistant P. aeruginosa clinical 
isolate (VR-143/97) from the U ni versity H ospital 
of Verona, Italy (3). This isolate was the index 
strain of an outbreak of biavm -positive 
P. aeruginosa, which was caused both by strains 
that were clonal ly related to VR-143/97 and by 
clonally unrelated strains (4). biavm isthesecond 
known metal I o-15-l actamase determinant that 
can spread among P. aeruginosa; the first was 
blaiMP, which was detected in the early 1990s in 
nosocomial isolates of various Enterobacteri- 
aceae P. aeruginosa, and other nonfastidious 
gram-negative nonfermenters from the Far East 
(l,2,5-7) and, recently, in an Ad netobacter 
baumannii clinical isolate from Italy (8). 
Although completely unrelated at the sequence 
level, biavm resembles blamp in being carried on 
an integron-borne mobile gene cassette and in 
encoding an enzyme (VIM-1) with broad 
substrate specificity (3). Because of these 
properties, biavm has the potential to become a 
dangerous resistance determi nant. 

An analysis of carbapenem-resistant 
P. aeruginosa from Italian hospitals since 1998 
showed production of metal lo-l5-l actamase, as- 
sayed as described (3), in five isolates from three 
hospitals in Italy. Two isolates (PPV-97 and PPV- 
108) were from the University Hospital of Pavia 
(PPV-97 was isolated in September 1998 from the 
urine of an inpatient in the neurosurgery 
department, and PPV-108 was isolated in 
November 1998 from a decubitus ulcer of an 
inpatient in the vascular surgery department); 


two (TS-832035 and TS-832347) were isolated in 
February 1999 from the University Hospital of 
T rieste (both from the blood of inpatients, in the 
intensive care unit and in the internal medicine 
department, respectively); and one (SAP-01/99) 
was isolated in September 1999 from the Rome 
University Hospital "PoliclinicoUmbertol "(from 
the blood of an inpatient in the vascular surgery 
department). Except for those from Pavia 
Hospital, where the two departments share the 
same surgical unit, noepidemiologic relationship 
could be established among any of them or with 
those previously isolated in Verona (3,4). 

The five isolates were highly resistant to 
carbapenems (M I Csfor i mi penem and meropenem 
were >64 pg/mL) and tocarbenicillin, ticarcillin, 
ticarcillin/clavulanate, piperacillin, piperacillin/ 
tazobactam, mezlocillin, ciprofloxacin, gentami- 
cin, tobramycin, and netilmicin. All fivewerealso 
resistant to ceftazidi me and cefepi me, except for 
SAP-01/99, which had intermediate resistanceto 
the above drugs. Some isolates retained 
suscepti bi I i ty to aztreonam (PPV-97, TS-832035, 
and SAP-01/99) or amikacin (PPV-108, TS- 
832035, and TS-832347). 

In a colony-blot hybridization assay (3), all 
the above isolates were recognized by a blavm- 
specific probe consisting of an amplicon that 
contained the entire biavm coding sequence (3). 
None were recognized by a probe specific for 
bl aiMP and consisti ng of a 0.5-kb H i ndl 1 1 fragment 
from the blamp gene (8). 

Our results indicate that circulation of 
carbapenem-resistant P. aeruginosa carrying 
biavm metal lo-15-l actamase determinants, origi- 
nally detected in one Italian hospital (3-4), could 
soon become widespread. The detection in 
different hospitals of blavm-positive isolates that 
apparently were epidemiological I y unrelated 
suggests that the environmental reservoir of 
blavm-containing strains is relatively broad and 
that this novel determinant has potential 
relevance for the emerging phenomenon of 
carbapenem resistance in P. aeruginosa. Since 
we did not sequence the blavm-related genes 
carried by the various isolates, we do not yet 
know whether they differ from that cloned from 
P. aeruginosa VR-143/97 (3). The five isolates 
descr i bed i n t h i s report a re bei ng ch a racter i zed to 
ascertain their clonal relatedness and identify 
the sequences of their blavm-related determi- 
nants. The recent appearance of this and other 
acquired metallo-15-lactamasesamongP. aerugi nosa 
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and other gram-negative pathogens i n E urope (8- 
10) underlines the need for systematic surveil- 
lance to monitor the spread of similar resistance 
determinants. 

This work was supported by grants no. FMRX-CT98- 
0232 from the EuropeanTMR Research Network on metallo- 
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Malaria and Global Warming 
in Perspective? 

To the Editor: I read with great interest the 
article "From Shakespeare to Defoe: malaria in 
England in the Little Ice Age" (1). Unfortunately, 
the article is not as balanced as a presentation 
last year by Paul Reiter, which clearly illustrated 
that, although climate is important in the 
transmission of malaria, the influence of other 
factors (e.g., access to medical care and improved 
housing) is likely to be of more importance in 
Europe. 

Malaria indeed was quitecommon in Europe, 
even in the Roman Empire and in Medieval 
Europe, and until a few decades ago, it was still 
present in parts of Europe, Australia, and North 
America. I n fact, the failure of the 1806 British 
invasion of Zeeland in the Netherlands may be 
attributableto infection of theBritish forces with 
malaria. Flowever, the authors referenced by 
Reiter have never made the claim that in the 
coming years warmer "temperatures will result 
in malaria transmission in Europe and North 
America." On the contrary, the reports of the 
Intergovernmental Panel on Climate Change 
Reiter quotes conclude that "Although climate 
change could increase the potential transmission 
of malaria [in Europe and North America], 
existing public health resources— disease sur- 
veillance, surface water management, and 
treatment of cases— would make reemergent 
malaria unlikely" (2,3). 

Reiter's argument that some scientists 
attribute the recent observed increase in malaria 
risk to climate trends is also not accurate. While 
acknowledging the sensitivity of the malaria 
mosquito and parasite to cl i mate, these research- 
ers examine insect and incidence data to explore 
multiple factors underlying malaria emergence. 
Another group of scientists uses mathematical 
simulation models to estimate changes in 
malaria risk over the next few decades. These 
models, which are heuristic tools not meant to 
predict future worlds, assess how potential risk 
for malaria may by affected by changes in 
climate (4). The goals of both types of research 
are to i mprove knowl edge of the compi ex mal ari a 
transmission cycle, define epidemic-prone areas, 
identify the reasons for increased malaria risk, 
and develop solutions to protect vulnerable 
communities. 
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Dr. Reiter acknowledges the sensitivity of 
malariatoclimaticinfluences, and I amsurethat 
he agrees that change in climate will affect risk 
for transmission— he may be skeptical as to 
whether global warming will ever become a fact, 
but t h at i s a not her quest i on . Wh i I e Rei ter 's paper 
offers an interesting perspectiveon thehistoryof 
malaria in Europe, it provides no illuminating 
i nformati on on the i nf I uence of cl i mate change on 
human health. 

Pirn Martens 

Maastricht University, Maastricht, the Netherlands 
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For P. Reiter’s response, please see 
http;//www.cdc.gov/ncidocl/E I D/vol6no4/reiter.htm 


Serologic Evidence of Human Monocytic 
and Granulocytic Ehrlichiosis in Israel 

To the Editor. We read with great attention the 
article by Dr. Keysary et al., who reported the 
first evidence of human monocytic and granulo- 
cytic ehrlichiosis in Israel (1); however, we 
disagree with their conclusions. 

Ehrlichiae comprise a large group of 
intracellular organisms pathogenic for animals 
and occasionally for humans. Because these 
organisms are closely related, serologic cross- 
reactions occur within and between groups, 
leading to mistakes in identification. For 
example. Ehrlichia chaffeensis was misdiag- 
nosed as E. canis in humans (2) and human 
granulocytic ehrlichiosis as human monocytic 


ehrlichiosis in areas where the vector was not 
present (3). Because of such cross-reactions, 
serology alone is not sufficient to establish the 
existence of a new ehrlichial disease. 

With the exception of Rhipicephalus 
sanguineus, the brown dog tick, which is 
distributed worldwide, tick species of medical 
importance are very geographically specific. For 
example, the I xodes and Dermacentor spp. found 
in Europe are not those found in the United 
States. Consequently, tick-transmitted organ- 
isms and diseases are also very specific 
geographically. For example, Borrdia spp. found 
i n the Old World are not found i n America (except 
for B. burgdorferi strictosensu, which isfound in 
both Europe and America). R. rickettsii, 
transmitted by Dermacentor andersoni and 

D. vari abl ills, is reported i n the U nited States but 
not in Europe, where the vectors are not present. 

American monocyticehrlichiosis iscaused by 

E. chaffeensis, which is transmitted by the tick 
Amblyomma americanum, found only in America. 
The main reservoir is the deer Odocoileus 
virgi nanus (4). 

It is very unlikely that a tick-borne disease 
occurred in a country where neither the vector 
nor the reservoir of the bacterium exists. All 
attempts to demonstrate the presence of 
E. chaffeensis in theOld World, including Africa, 
have failed. Indeed, there is no convincing 
evi dence of the existence of E. chaffeensis outside 
America. 

Philippe Brouqui* and J . Steven Dumlerf 

Unite des Rickettsies, Faculte de Medecine, 
Marseille, France; and J ohns Hopkins University 
School of Hygiene&i Public Health, Baltimore, 
Maryland, USA 

References 

1. Keysary A, Amram L, Keren G, Sthoeger Z, Potasman 
l,J acobA.etal. Serologic evidenceof human monocytic 
and granulocytic ehrlichiosis in Israel. Emerg Infect 
Dis 1999;5;775-8. 

2. Maeda K, Markowitz N, Hawley RC, Ristic M, Cox D, 
McDadeJ E. Human infection with Ehrlichia canis, a 
leukocytic rickettsia. N EnglJ M ed 1987;316;853-6. 

3. BrouquiP.RaoultD.Humanehrlichiosis.N EnglJ Med 
1994;330; 1760-1. 

4. DumlerJ S, BakkenJ S. Ehrlichial diseases of humans; 
emerging tick-borne infections. Clin Infect Dis 
1995;20;1102-10. 


Emerging Infectious Diseases 


314 


Vol. 6, No. 3, May-]une 2000 



Letters 


Serologic Evidence of Human Monocytic 
and Granulocytic Ehrlichiosis in Israel 

To the Editor. Our articleon ehrlichiosis in Israel 
is one of several reporting serologic evidence of 
ehrlichiosis outside North America (1-4). Sero- 
logic tests were performed for immunoglobulin 
(IgG) antibodies to Ehrlichia chaffeensis, E. 
canis, and human granulocytic ehrlichiosis 
(HGE) agent to prevent misinterpretation due to 
cross-reaction. The conclusions meet the criteria 
published by the American Society for 
Rickettsiology, in which a confirmed diagnosisof 
human monocyticehrlichiosis (H M E) is based on 
a "single serum titer of 256" in a patient with 
cl i n i ca 1 1 y com pat i bl e d i sease ( 5) . 

The argument that the same epidemiologic 
circumstances that exist in the United States 
have to be prevalent in Israel for the disease to be 
present is not compelling. Reporting serologic 
evidence of ehrlichiosis in Israel alerts physi- 
cians to the possibility of H M E and HGE when 
they see patients with symptoms compatible 
with these diseases. We agree that isolation and 


molecular identification of both agents are 
essential to confirming the presence of these 
diseases in Israel. 

Avi Keysary and Trevor Waner 

Israel Institute for Biological Research, 

Ness Ziona, I srael 
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Upcoming Events 


Parasite Genomes — Thematic issue 
Now Avaiiabie 

The International J ournal for Parasitology, 
Volume 30/4 (April 2000), announces the 
publication of a thematic issue on Parasite 
Genomes, guest-edited by Paul Brindley at 
Tulane University, New Orleans, USA. 

Experts in the field have contributed articles 
on medically and economically important para- 
sites, both helminths and protozoans. Topics 
include 'The African trypanosome genome," by 
Najib M. El-Sayed, Priti Hegde, John 
Quackenbush, Sara E. Melville, and J ohn E. 
Donelson; "Genes and genomes of Necator 
americanusand related hookworms," by Mark L. 
Blaxter; 'The conserved genome organisation of 
non-falciparum malaria species: the need to 
know more," by Leonard H .M . van Linn, Chris] . 
J anse, and Andrew P. Waters; "Mitochondrial 
genome diversity in parasites," by Jean E. 
Feagin; and 'The Giardia lamblia genome," by 
Rodney D. Adam. 

F ree access to the full text of the 20 articles 
can be found at http://www.elsevier.nl/locate/ 
ijpara. 


5th internationai Conference on Legioneiia 
Ulm, Germany 
September 26-29, 2000 

The 5th International Conference on 
Legionella will be held in Ulm, Germany, 
September 26-29, 2000. Topics include pathogen- 
esis, immunology, ecology, clinical microbiology, 
epidemiology, surveillance, and prevention. 
Proceedings will be published by the American 
Soci ety for M i crobi ol ogy. 

For more information, contact] ensThomsen, 
Department of Medicine, Microbiology and 
Flygiene, University of Ulm, 89081 Ulm, 
Germany; telephone: 49-731-50-24603; fax: 49- 
731-50-24619; e-mail :]ens.thomsen(§medizin. uni- 
ulm.de; http://www.uni-ulm.de/legionella/. 


Third internet Conference on Saiivarian 
Trypanosomes and Trypanosomatids 
October 2-1 8, 2000 

The third I nternet Conference on Saiivarian 
Trypanosomes and Trypanosomatids will take 
place October 2-18, 2000, at http://www. 
dbbm.fiocruz.br/trypnews/events/ticstt.html or 
http://pubweb.nwu.edu/~kmt564/TI CSTT.html. 
Through featured presentations, posters, and 
discussion boards, the conference will cover all 
aspects of trypanosome research, including 
biology and ultrastructure, biochemistry and 
drug development, immunology and pathology, 
molecular biology, epidemiology, vectors, and 
control. Proceedings will be published in the 
I nternational J ournal for Parasitology. 

For presentations (6 megabyte limit), the 
deadlines are September Ifor FITML format and 
August 15 for other formats. For posters (1.5 
megabyte limit), the deadline is August 15; only 
FITML format will be accepted. Instructions for 
authors are available on-line at the above URL. 
For more information, contact Alberto Davila, 

I nstitutoOswaldoCruz, Riode] aneiro, Brazil; e- 
mail: davila(§gene.dbbm.fiocruz.br; or Kevin 
Tyler, Northwestern University, Chicago; e-mail: 
k-tyler (gnwu.edu. 


Nosocomiai infections Conference 
Institut Pasteur, Paris, France 
November 15-18, 2000 

The rapid proliferation of nosocomial infec- 
tions, mostly in intensive care and surgical units, 
is related to increases in surgical interventions, 
multidrug-resistant microorganisms, and the 
number of immunocompromised patients at high 
risk. Prophylaxis and treatment may be delayed 
because of poor diagnosis, biofilm formation on 
implants, and the emergence of new pathogens 
from bacterial genetic exchange and in vivo 
mutagenesis. 

T opi cs to be discussed during the N osocomi al 
I nfections Conference include surveillance, new 
methods for rapid microorganism identification 
and typing, identification of genes linked to 
antibiotic resistance, molecular studies of 
pathogenicity, and biofilm formation. 
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The I nstitut Pasteur Eur<xonferences are 
intended to facilitate communication between 
academicians and industry. Committee members 
for the Nosocomial Infections Conference are 
Patrice Courvalin, J ean-Christophe Lucet, 
Philippe Sansonetti, and Cecile Wandersman. 

For more information, contact Ludovic Drye, 
Conference Coordinator, I nstitut Pasteur 
EuroConferences, 28 ruedu Docteur Roux, 75724 
Paris Cedex 15, France; telephone: 331-406- 
13374; fax: 331-406-13405; e-mai I : ldrye(gpasteur.fr; 
conference web site: http://www.pasteur.fr/appli- 
cations/euroconf/. 


Erratum Vol. 5, No. 5 

In the article "Using a Spatial Filter and a 
Geographic Information System to Improve Rabies 
Surveillance Data," by Andrew Curtis, there is an error 
in Figure 1 on page 604. The topmost and center maps 
have been transposed. The on-line article contains the 
correct figure, at URL: http://www.cdc.gov/ncidod/eid/ 
vol 5no5/cu rtis.htm 

We regret any confusion this error may have 
caused. 
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Editorial Policy and Call for Articles 


Emerging I nfectious Diseases is a peer-reviewed journai estabiished expressiy to promote the recognition of new and reemerging infectious 
diseases around the worid and improve the understanding of factors invoived in disease emergence, prevention, and eii mi nation. 

Thejournai has an internationai scopeand is intended for professionais in infectious diseases and reiated sciences. We weicome contributions 
from infectious disease speciaiists in academia, industry, dinicai practice, and pubiic heaith, as weii as from speciaiists in economics, demography, 
socioiogy, and other disci pi ines. Inquiries about thesuitabiiity of proposed articies may be directed to the Editor at 404-371-5329 (tei), 404-371-5449 
(fax), or eideditor@cdc.gov (e-maii). 

E merging I nfectious Diseases is pubiished in Engiish and features thefoliowing types of articies: Perspectives, Synopses, Research Studies, Policy 
Reviews, and Dispatches. The purpose and requirements of each type of articie are described in detaii below. To expedite pu hi i cation of information, 
we post journai articies on the Internet as soon as they are cleared and edited. 

Spanish and French translations of some articles can be accessed through thejournai 's homepage at www.cdc.gov/eid. Articles by authors from 
non-English-speaking countries can be made si multaneously aval lablein English and in theauthor'snativelanguage(electronic version of thejournai 
only). 


Instructions to Authors 

Manuscript Preparation 

Follow "Uniform Requirements for Manuscripts Submitted to 
Biomedical Journals" (Ann Intern Med 1997:126[l]36-47) (http:// 
www.acponline.org/journals/annals/01jan97/unifreqr.htm). 

Begin each of thefoliowing sections on a new pageand in this order: 
title page, abstract, text, acknowledgments, references, tables, figure 
legends, and figures. 

Title page. Give complete information about each author (i.e, full 
name, graduate degree(s), affiliation, andthenameoftheinstitution in 
which the work was done). Also provide address for correspondence 
(include fax number and e-mail address). 

Abstract and key words. Avoid citing references in the abstract. 
Include up to 10 key words; use terms listed in the Medical Subject 
Fleadings from Index Medicus (http://www.nlm.nih.gov/mesh/ 
meshhome.html). 

Text. Double-space everything, including the title page, abstract, 
references, tables, andfigurelegends. Typeonly on oneside of the paper 
and number all pages, beginning with the title page. Indent paragraphs 
5 spaces: leave no extra space between paragraphs. After a period, leave 
only onespace before beginning thenext sentence. Use Courier font size 
10 and ragged right margins. Italicize (rather than underline) scientific 
names when needed. 

Electronic formats. For word processing, use WordPerfect or MS 
Word. Send graphics in either .TIF (Tagged Image File), or .EPS 
(Encapsulated Postscript) formats. The preferred font for graphics files 
is Flelvetica. Convert Macintosh files into one of the suggested formats. 
Submit slides or photographs in glossy, camera-ready photographic 
prints. 

References. F ol I ow the U ni form Requi rements styl e. P I ace reference 
numbers in parentheses, not in superscripts. Number citations in order of 
appearance (including in text, figures, and tables). Cite personal 
communications, unpublished data, and manuscripts in preparation or 
submitted for publication in parentheses in text. Consult List ofj ournals 
Indexed in Index Medicus for accepted journal abbreviations: if ajournal 
is not listed, spell out thejournai title in full. List the first six authors 
followed by "et al." 

Tables and figures. Create tables within the word processing 
program's table feature (not columns and tabs within the word 
processing program). For figures, use col or as needed; send files, slides, 
photographs, or prints. Figures, symbols, lettering, and numbering 
should be clear and large enough to remain legible when reduced. Place 
figure keys within the figure. 

Manuscript Submission 

Include a cover letter verifying that the final manuscript has been 
seen and approved by all authors. 

Submit three copies of the original manuscript with three sets of 
original figures and an electronic copy (on diskette or by e-mail) to the 
Editor, Emerging Infectious Diseases, Centers for Disease Control and 
Prevention, IBOOClifton Rd., MS D 61, Atlanta, GA 30333, USA; e-mail 
eideditor@cdc.gov. 


Types of Articles 

Perspectives, Synopses, Research Studies, and Poiicy Reviews: 

Articles should be approximately 3,500 words and should include refer- 
ences, not to exceed 40. Useof subheadings in the main body of the text 
is recommended. Photographs and illustrations are encouraged. Pro- 
vide a short abstract (150 words) and a brief biographical sketch. 
Perspectives: Articles in this section should provide insightful 
analysis and commentary about new and reemergi ng i nfectious diseases 
or related issues. Perspectives may also address factors known to 
influencetheemergenceof diseases, including microbial adaptation and 
change: human demographics and behavior; technology and industry: 
economic development and land use; international travel and 
commerce: and the breakdown of public health measures. If detailed 
methods are included, a separate section on experimental procedures 
should immediately follow the body of the text. 

Synopses: This section comprises concise reviews of infectious diseases 
or closely related topics. Preference is given to reviews of new and 
emerging diseases; however, timely updates of other diseases or topics 
are also welcome. Use of subheadings in the main body of the text is 
recommended. If detailed methods are included, a separate section on 
experimental procedures should immediatelyfollowthebody ofthetext. 
Photographs and illustrations are encouraged. 

Research Studies: These articles report laboratory and epidemiologic 
results within a public health perspective. Although these reports may 
be written in the style of traditional research articles, they should 
explain the value of the research in public health terms and place the 
findings in a larger perspective (eg., "FI ere is what we found, and here 
is what thefindings mean"). 

Policy Reviews: Articles in this section report public health policies 
that are based on research and analysis of emerging disease issues. 
Dispatches: These brief articles are updates on infectious disease 
trends and research. The articles include descriptions of new methcds 
for detecting, characterizing, or subtyping new or reemerging 
pathogens. Developments in antimicrobial drugs, vacci nes, or infectious 
disease prevention or elimination programs are appropriate. Case 
reports are also welcome. Dispatches (1,000 to 1,500 words) need not be 
divided intosections. Providea short abstract (50 words); references, not 
to exceed 10; figures or illustrations, not to exceed two; and a brief 
biographical sketch. 

Book Reviews: Short reviews (250 to 500 words) of recently published 
books on emerging disease issues are welcome. 

Letters: This section includes letters that give preliminary data or 
comment on published articles. Letters (500 to 1,000 words) shouldnot 
be divided into sections, nor should they contain figures or tables. 
References (not more than 10) may be included. 

News and Notes: We welcome brief announcements (50 to 150 words) 
oftimely events of interest toour readers. (Announcements can be posted 
on thejournai web page only, depending on the event date.) In this sec- 
tion, we also include summaries (500 to 1,500 words) of conferences 
focusing on emerging infectious diseases. Summaries may provide ref- 
erences to a full report of conference activities and should focus on the 
meeting's content. 



